New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3944978	NA (NA)	6 (0.090)	4/272	9.2	24.8%	noncoding (1275/2904 nt)	rrlC	23S ribosomal RNA of rrnC operon
	?	NC_000913	= 4038811	19 (0.290)					noncoding (1293/2905 nt)	rrlA	23S ribosomal RNA of rrnA operon
Rejected: Coverage evenness skew score above cutoff.

ACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCC  >  NC_000913/4038669‑4038868                                                                                                                                               |                                                          aCCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCt                                                           >  2:811683/1‑143 (MQ=255) aCCGAAGCTGCGGCAGCGACACTATGTGTTGGTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCt                                                           <  1:375056/143‑1 (MQ=255)  ccGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTc                                                          >  1:347541/1‑143 (MQ=255)    gAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGCTAAAGCGGGTGAGACGCCCGCTCGc                                                        >  1:417384/1‑143 (MQ=255)      aGCTGCGGAGGAGACGCTAGGTGTTGTTGGGTAGGGGGGGGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCg                                                      <  2:133071/143‑1 (MQ=11)         tGCGGTAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGaa                                                   <  1:1055056/143‑1 (MQ=255)             gCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAGAGCCCGCTCGCCGGAAGAcc                                               >  1:741150/1‑143 (MQ=255)             gCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGATGTGCTGGGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATACGTAACGATGAAGCGGGGGAAAAGGCCGCCCGCCCGACGAcc                                               >  2:510933/1‑143 (MQ=14)               aGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCaa                                             <  1:836506/143‑1 (MQ=40)               aGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGGGAAGGTGTGCTGTGCGGCATGCTGGAGGTATCAGAAGTGCGAATGCGGCCATGAGTAACGAAAAAGCGGGCGAACAGCCCGCTCGCCGGAAGACGaa                                             >  2:1035791/1‑143 (MQ=14)               aGCGACACCATGTGGTGTTGGGTAGGGCAGCCTTCTGTAACGCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGTTATCAGAAGGGCGAATGCTGACATAAGTAACGATAAAGCGGGCGAAAAGCCCGCTCGCCGGAAGACCaa                                             <  1:473766/143‑1 (MQ=14)                gagatacTATGTGTTGTTGGGTAGGGGAGCGTTCTGGAAGCCTGTGAAGTTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAg                                            <  2:699011/138‑1 (MQ=255)                 cGACACTATGTGTTGTTGGGTAGGGGCGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGGGGCAGCAGAAGTGAGAATGCTGACCTAAGTAGCGATCAAGCGGGAGAAAAGCGCGCGCGCCTGGAGACCAAgg                                           >  2:630870/1‑143 (MQ=14)                   acacTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGt                                         >  1:693516/1‑143 (MQ=38)                         tgtgTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTACCGATCAAGCGGGTGAGACGCCCGCTCGCCGGAAGACCAAGGGTTCCTGt                                   >  1:745959/1‑143 (MQ=255)                           tgttgtTGGGTAGGGGGGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTcc                                 >  1:658471/1‑143 (MQ=18)                              cgtcgggTAGGGGAGCGATCGGTAAGCCTGTGAAGGTTTGCGGGGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGTCATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTTGCCGGAAGACCAAGGGTGCCTGTCCAAc                              <  2:455292/139‑1 (MQ=11)                                ttgGGTAGGGGAGCGTTCGGTAAGCCGGTGAGGGTGTGCTGTGAGGCATGCTGGGGGGATCAGAAGTGCGAAGGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGt                            <  2:739146/143‑1 (MQ=11)                                                tCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTg            >  1:1035609/1‑143 (MQ=32)                                                    tAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAACGCCCGCTCGCCGGAAGACCAAGGGGTCCTGTCCAACGTTAGTCGGGGCAGGGTGAggc        >  1:579733/1‑141 (MQ=12)                                                          tgtgAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGAcccc  >  1:6311/1‑143 (MQ=32)                                                                                                                                               |                                                          ACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCC  >  NC_000913/4038669‑4038868

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.