New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4038798	38 (0.720)	11 (0.220)	4/272	7.3	25.2%	noncoding (1280/2905 nt)	rrlA	23S ribosomal RNA of rrnA operon
	?	NC_000913	= 4038805	29 (0.570)					noncoding (1287/2905 nt)	rrlA	23S ribosomal RNA of rrnA operon
Rejected: Coverage evenness skew score above cutoff.

CCATGCACCGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCC  >  NC_000913/4038663‑4038867                                                                                                                                                |                                                               ccATGCACCGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGCTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAg                                                                >  1:403546/1‑143 (MQ=255)   aTGCACCGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCGTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAGGcc                                                              >  1:888681/1‑143 (MQ=255)    ggCACGGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGccc                                                             <  1:250943/142‑1 (MQ=255)     gCACCGAAGCTGCGGCAGCGACACTGTGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCg                                                            <  1:531112/143‑1 (MQ=255)     gCACCGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGACTGCTTACATGAGTAACGATAAAGCGGGTGAAAAGCCCg                                                            >  1:274050/1‑143 (MQ=255)        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)                                                                                                                                                |                                                               CCATGCACCGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCC  >  NC_000913/4038663‑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.