New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 227052	33 (0.490)	17 (0.260)	8/270	7.5	35.0%	noncoding (1294/2904 nt)	rrlH	23S ribosomal RNA of rrnH operon
	?	NC_000913	2727908 =	NA (NA)					noncoding (1275/2904 nt)	rrlG	23S ribosomal RNA of rrnG operon
Rejected: Coverage evenness skew score above cutoff.

AAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTAAGGCGAGGCCGAAAG  >  NC_000913/226913‑227126                                                                                                                                            |                                                                           aaGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCtcgctcgc                                                                         <  2:1059529/143‑1 (MQ=38)    cggcggcggcgACGCGTATGCGTTGTTGGGGAGGGGAGCGTTCTGTAAGCCTGTTAAGGTGTACTGTGAGGTATGCGGGAGGTATCAGAAGTGCGAGTGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCgg                                                                      <  2:743144/141‑1 (MQ=11)      gctgCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGaa                                                                    <  1:150670/140‑1 (MQ=38)        ggCAGCGACGCTTATGCGTTGTTGGGTAGGAGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGa                                                                  >  2:671422/1‑143 (MQ=37)          cAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGAcc                                                                <  1:832739/143‑1 (MQ=38)             cGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCagg                                                             >  1:817188/1‑141 (MQ=38)                cggctgtGCCTTTTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACCGTCAGTTATGCTGGAGGTATCATAAGTGCGAATGCTGACATCAGTAACGATAAAGCGGGTGAAACGCCCGCTCGCCGGAAGACCAAGGGt                                                          <  1:1009439/137‑1 (MQ=11)                 gCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGtt                                                         <  2:272090/143‑1 (MQ=37)                 gCTTATGCGTTGTTGGGTAGGGGAGCGTCCTGTCACTCTGTGAAGGTGTAATGTGAGGGATGCTGGAGGTATCAGAAGTGCGAACGCTGACATAAGGAAAGAGAAAGCGGGGGACAAGCACGCCCGCCTGAAgggccagccgc                                                         >  2:533454/1‑133 (MQ=2)                    tATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCt                                                      >  1:477794/1‑143 (MQ=38)                      tGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGGGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGGGGGAGAAAAGCCAGCTAGCCGGAAGACCAGGGGCTGCTGt                                                    >  2:410985/1‑143 (MQ=11)                       gCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTc                                                   <  1:318740/143‑1 (MQ=37)                       gCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTc                                                   >  1:1054507/1‑143 (MQ=37)                        cGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACAAAGGGTTCCTGTcc                                                  >  2:421953/1‑143 (MQ=25)                            gttgGGTAGTGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACg                                              <  1:5937/143‑1 (MQ=25)                            gttgGGTAGGGCAGCTTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACg                                              <  1:453348/143‑1 (MQ=21)                             ttgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTAGCCGGAAGACCAAGGGTTCCTGTCCAACGt                                             >  1:438576/1‑143 (MQ=25)                             ttgGGGAAGGGAGAGCTCTGTAAGCCTGTGAAGGTGTACTGTGAGGGATGCTGGAGGGATCAGAAGTGCGACTGCTGACATAAGTGACGATAGAGGGGGTGAGAAGGCCGCTCGCCGGAAGACCAAGGGTACCTGTCCAGCGt                                             >  2:421742/1‑143 (MQ=11)                               gggTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTa                                           <  1:1030897/143‑1 (MQ=35)                                       gAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGCTAGTCGGGGc                                   >  2:43645/1‑143 (MQ=25)                                         gCGTTCTGCAAGCCTGTGAAGGTGTACTGCGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAg                                 <  1:19743/143‑1 (MQ=21)                                           gtgctgggaggcTGGGAAGGTGTACGGTGAGGTATGCGGGAGGTATCGGGAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAggg                               <  1:393461/132‑1 (MQ=11)                                                gcaaGCCCGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGt                          <  1:136512/141‑1 (MQ=21)                                                gTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGt                          >  1:1065524/1‑143 (MQ=35)                                                    gCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACAAAGGATTGCTGTCCAACGTTAATCGGGGCAGGGTGAGTccgc                      >  2:1128319/1‑140 (MQ=255)                                                     ccTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAAGGATGAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTACACCGCTCATCGGTGCAGGGTCAGCCGAca                     >  2:211534/1‑142 (MQ=11)                                                         tgAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAGCGATAAAGCGGGTGAAACGCCCGCTCGCCGGAAGACGAAGGGTTCGTGTCCAACGTTAATCGGGGCAGCGCGCGTCGCCCCCTa                 >  2:1124690/1‑143 (MQ=11)                                                                  aCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTGCTGTCCAACGTGAATCGGGGCAGGGGGAGTCGACCCCTAAGGCgcagc        >  2:937140/1‑139 (MQ=255)                                                                  aCTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTAAGGCGAGGc        >  2:931887/1‑143 (MQ=32)                                                                    tgtgAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACTGAGGGTTGCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTAAGGCGAGGCCg      >  1:381906/1‑143 (MQ=12)                                                                     gtgAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTAAGGCGAGGCCGa     <  1:463578/143‑1 (MQ=32)                                                                      tgAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGGCGACCCCTAAGGCGAGGCCGaa    >  2:1098449/1‑143 (MQ=17)                                                                        aGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTAAGTCGCGGCCGaagg  >  1:13319/1‑141 (MQ=12)                                                                                                                                            |                                                                           AAGCTGCGGCAGCGACGCTTATGCGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTAAGGCGAGGCCGAAAG  >  NC_000913/226913‑227126

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.