New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3944978	NA (NA)	13 (0.260)	8/272	5.3	31.9%	noncoding (1275/2904 nt)	rrlC	23S ribosomal RNA of rrnC operon
	?	NC_000913	= 4038811	29 (0.550)					noncoding (1293/2905 nt)	rrlA	23S ribosomal RNA of rrnA operon
Rejected: Coverage evenness skew score above cutoff.

CCGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCC  >  NC_000913/4038670‑4038867                                                                                                                                               |                                                         ccGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTc                                                         >  2:232298/1‑143 (MQ=255)       cttcGGTAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGa                                                   <  2:753486/140‑1 (MQ=255)       cTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGa                                                   <  2:574529/143‑1 (MQ=255)        tGCGGGAGCGACA‑TTTGTGGTGGTTGGGTGGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGaa                                                  <  2:731603/143‑1 (MQ=11)        tGCGGCAGCGAGACGCGGT‑GTTGTTGGGTAGGGGAGCGTGCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGaa                                                  <  1:164127/143‑1 (MQ=11)             cAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATCATGGCACCGCCATATGTGCGCCTCGATACGTAAGTAACGATAAAGCGGGTGGAGAGCCCGCTCGCCGGAAGACCa                                             >  2:423184/1‑143 (MQ=9)              aGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCaa                                            <  2:750844/143‑1 (MQ=255)                cGACGCTATGG‑GGTGTTTGGTAGGGGAGCGTTCTGGAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAgg                                          <  2:930706/143‑1 (MQ=11)                cGACACTATGG‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAgg                                          <  1:328327/143‑1 (MQ=38)                 gACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGCTAAAGCGGGTGAAAAGCCCGCTCGCCGCCAGACCAAggg                                         >  1:897296/1‑143 (MQ=255)                  acacTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGt                                        >  1:925933/1‑143 (MQ=38)                    acTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTCTGAACGTGTGCTGTGAGGCATGCTGGAGGCACCAGAAGCGCGAATCCCTACATAAGTAACGATAAAGCGCGCAAGAAGCCCCCCCCCCGCAACACCAAgcgccc                                      >  2:109310/1‑138 (MQ=2)                      tATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCACGGGTGCCt                                    >  2:206065/1‑143 (MQ=25)                            gttgttGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGCGGTCTCCGAACTCCGAATGCTGACATAAGTAACGATAAAGCGGGTGACAAGCCCGCTCGCCGGAAGCCCAAGGGATCCTGTCCa                               >  1:342071/1‑143 (MQ=11)                              tgttgGTGAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACTATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAAc                             <  1:107293/143‑1 (MQ=12)                              tggtgGGTAGGGGGGCGTTCTGTAGGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAAc                             <  2:605705/140‑1 (MQ=12)                               gtgtggTAGGGGAGCGGTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACg                            <  2:31980/139‑1 (MQ=17)                                ttgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAACGGTTCCTGTCCAACGt                           >  2:784254/1‑143 (MQ=17)                                      aGGGGAGCGTTCTGTAAGCTCGTGAAGGTGTGCTGTGAAGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGACAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCg                     <  1:660471/143‑1 (MQ=255)                                      aGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCg                     >  1:928085/1‑143 (MQ=32)                                        gggAGCGTTCTGTGAGCGTGTGATGGTGTGCTGTGAGTCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCggg                   <  2:567319/143‑1 (MQ=255)                                        gggAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCggg                   <  1:609341/143‑1 (MQ=32)                                             cggtCTGGAAGCTTGTGACGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAgg              <  2:31806/140‑1 (MQ=255)                                             cGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAgg              >  2:318567/1‑143 (MQ=32)                                              gtgcTGTAAGCCTGTGAACGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAggg             <  2:667624/140‑1 (MQ=17)                                                taggTAAGGCTGTGAAGGGGTGCGGTGAGGCATGCGGGAGGGAGCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTg           <  1:775359/140‑1 (MQ=11)                                                       gCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCCACGTGAATCGGGGCAGGGTGAGTCGAc    >  2:915539/1‑143 (MQ=17)                                                       gCCTGTGAAGGTGAGCTCTTTTGCATTCTGGAGGTAGCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGAc    <  2:150335/143‑1 (MQ=11)                                                         cTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGAccc  <  1:318567/143‑1 (MQ=32)                                                                                                                                               |                                                         CCGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCC  >  NC_000913/4038670‑4038867

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.