New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4168341	33 (0.210)	27 (0.200)	18/178	13.9	47.7%	intergenic (+141/‑31)	rrsB/gltT	16S ribosomal RNA of rrnB operon/tRNA‑Glu
	?	NC_000913	= 4168363	NA (NA)					intergenic (+163/‑9)	rrsB/gltT	16S ribosomal RNA of rrnB operon/tRNA‑Glu
Rejected: Coverage evenness skew score above cutoff.

GATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTC  >  NC_000913/4168242‑4168408                                                                                                    |                                                                    gATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAat                                                                     >  2:756945/1‑100 (MQ=40) gATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAat                                                                     >  2:3727014/1‑100 (MQ=40) gATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAat                                                                     >  2:3522774/1‑100 (MQ=40)  aTAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAata                                                                    >  1:609827/1‑100 (MQ=38)  aTAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAat                                                                     <  1:3001865/99‑1 (MQ=38)   tAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAattg                                                                   >  2:3581364/1‑98 (MQ=38)   tAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAatat                                                                   >  2:1772315/1‑100 (MQ=38)   tAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTGCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAatat                                                                   >  1:893475/1‑100 (MQ=25)   tAGAAAGTGAAAAGCAAGGCGTTTACGAGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTGCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAatat                                                                   <  1:745024/100‑1 (MQ=255)     gAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGGGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCa                                                                 <  1:3641826/100‑1 (MQ=25)     gAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAattgct                                                                 >  2:98293/1‑96 (MQ=25)     gAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCa                                                                 <  1:1037134/100‑1 (MQ=35)     gAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCa                                                                 >  1:1635190/1‑100 (MQ=35)     gAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCa                                                                 <  2:3192148/100‑1 (MQ=35)     gAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCa                                                                 >  1:2124549/1‑100 (MQ=35)     gAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCa                                                                 >  1:1983707/1‑100 (MQ=35)     gAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCa                                                                 >  2:1293907/1‑100 (MQ=35)      aaaGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCAc                                                                >  1:1472507/1‑100 (MQ=33)      aaaGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCAc                                                                <  1:2982406/100‑1 (MQ=33)      aaaGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCAc                                                                <  2:1646697/100‑1 (MQ=33)       aaGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACg                                                               <  2:2148375/100‑1 (MQ=33)       aaGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACg                                                               <  1:2435970/100‑1 (MQ=33)        aGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGc                                                              <  2:3845197/100‑1 (MQ=33)        aGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGc                                                              <  1:1656441/100‑1 (MQ=33)        aGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGc                                                              >  2:3384774/1‑100 (MQ=33)        aGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCCATTAATGAAAGCTCACCCTACACGAAAATATCACGc                                                              <  1:2731067/100‑1 (MQ=18)         gTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAACTATCACGCa                                                             >  2:220243/1‑100 (MQ=18)         gTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCa                                                             >  1:1404213/1‑100 (MQ=33)         gTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCa                                                             >  2:131264/1‑100 (MQ=33)          tGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCaa                                                            >  2:745694/1‑100 (MQ=32)           gAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAAc                                                           >  1:3771053/1‑100 (MQ=32)           gAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAAc                                                           >  1:3137651/1‑100 (MQ=32)           gAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAAc                                                           >  1:1144519/1‑100 (MQ=32)            aaaaGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAAcg                                                          >  1:1712053/1‑100 (MQ=32)             aaaGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAAcgc                                                         >  2:3170260/1‑100 (MQ=31)              aaGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGcgttgcg                                                        <  2:893472/100‑7 (MQ=2)              aaGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAAcgcg                                                        <  1:605798/100‑1 (MQ=11)                    ggggTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGa                                                     <  2:3296464/94‑1 (MQ=1)                                                                     tcgaTTAATGAAAGCTCACCCTACACGAGAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGCCACCGCCCTTTc  <  2:1289362/96‑1 (MQ=0)                                                                                                    |                                                                    GATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAATATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTC  >  NC_000913/4168242‑4168408

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