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

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

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.