New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4037477	NA (NA)	16 (0.180)	5/170	NT	33.0%	intergenic (+142/‑42)	alaT/rrlA	tRNA‑Ala/23S ribosomal RNA of rrnA operon
	?	NC_000913	= 4168637	34 (0.360)					intergenic (+190/‑4)	gltT/rrlB	tRNA‑Glu/23S ribosomal RNA of rrnB operon

GTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATC  >  NC_000913/4037383‑4037494                                                                                               |                  gTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCTAGCTGAAAATTGATACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAaca              <  2:2106346/100‑1 (MQ=2)   tCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACg            <  1:1872767/100‑1 (MQ=32)   tCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACg            >  2:2511393/1‑100 (MQ=32)    cACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTTAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACga           <  1:851573/100‑1 (MQ=17)      cTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatg         >  1:1699475/1‑100 (MQ=33)      cTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatg         >  1:734503/1‑100 (MQ=33)      cTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatg         <  2:1057243/100‑1 (MQ=33)      cTGCGAAGTTTTGCTCTTCAAACATCTGGATCTAGCTGAAAATTGAAACACTGACCACCGCAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatg         <  2:650627/100‑1 (MQ=11)        gCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatgat       >  2:1312589/1‑100 (MQ=32)        gCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatgat       >  2:1873138/1‑100 (MQ=33)        gCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatgat       >  2:712508/1‑100 (MQ=33)        gCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatgat       >  2:712733/1‑100 (MQ=33)        gCGAAGTTTTGCTCTTTAAAAATCTGAAACAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatga        <  2:337743/99‑1 (MQ=16)          gAAGTTTTGCTCTTTAAAAATCTGGATTAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGAtgatga     <  1:1904077/100‑1 (MQ=18)           aaGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGaa    <  2:2143113/100‑1 (MQ=33)             gTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATc  <  1:1111752/100‑1 (MQ=31)                                                                                               |                  GTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATC  >  NC_000913/4037383‑4037494

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