New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3944978	NA (NA)	12 (0.200)	9/186	6.0	42.9%	noncoding (1275/2904 nt)	rrlC	23S ribosomal RNA of rrnC operon
	?	NC_000913	= 4038811	17 (0.260)					noncoding (1293/2905 nt)	rrlA	23S ribosomal RNA of rrnA operon
Rejected: Coverage evenness skew score above cutoff.

GCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAG  >  NC_000913/4038712‑4038826                                                                                                    |                gCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCt                 >  1:826160/1‑100 (MQ=33)   gTTCTGTAAGCCTGTGAATGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCg               <  2:876813/100‑1 (MQ=17)   gTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCg               >  1:1351207/1‑100 (MQ=33)   gTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCg               >  1:396712/1‑100 (MQ=33)    ttCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGc              >  2:978699/1‑100 (MQ=32)       tGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCgg           <  2:996156/100‑1 (MQ=32)        gTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGa          <  1:1102482/100‑1 (MQ=33)        gTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGa          >  1:620788/1‑100 (MQ=32)         tAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGaa         <  2:360723/100‑1 (MQ=32)          aaGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAg        >  1:1146096/1‑100 (MQ=33)           aGCCTGTGAAGTTCTTCTCGCAGGCATGCTGGAGTTAGCAGAAGTGCGAATGCTCACATAAGTAACGATAAAGCGGGTCAAAAGCCCGCTCGCCgg           <  2:1207813/96‑1 (MQ=11)           aGCCTGTGAAGGTGTGCTGTGAAGGATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTCACGAGAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGa       >  2:920201/1‑100 (MQ=11)            gCCGGTGAAGGTGTGATGTAAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTAAAAAGCCCGCTCGCCGGaaa        <  2:1250120/98‑2 (MQ=11)              cTGTAAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCa    <  2:372072/100‑1 (MQ=14)               tgtgAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCaa   >  2:1465508/1‑100 (MQ=32)               tgtgAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCaa   >  2:1468899/1‑100 (MQ=32)                gtgAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAg  >  2:1073200/1‑100 (MQ=31)                gtgAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAg  <  1:1547511/100‑1 (MQ=31)                                                                                                    |                GCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAG  >  NC_000913/4038712‑4038826

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ 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.