New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4037450	NA (NA)	20 (0.250)	10/164	NT	18.9%	intergenic (+115/‑69)	alaT/rrlA	tRNA‑Ala/23S ribosomal RNA of rrnA operon
	?	NC_000913	= 4168568	93 (1.080)					intergenic (+121/‑73)	gltT/rrlB	tRNA‑Glu/23S ribosomal RNA of rrnB operon

GTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCG  >  NC_000913/4037383‑4037495                                                                    |                                              gTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAaca               <  2:1724398/100‑1 (MQ=11)    cACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACga            >  1:2041722/1‑100 (MQ=33)    cACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACga            >  2:867640/1‑100 (MQ=33)      cTGCGAAGTTTTGCTCTTTAAAAACCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAACTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgat           >  2:695283/1‑99 (MQ=12)        gCGAAGTTTTGCTCTTTAAAAATCTGGATCGAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAAAACgatga         >  2:1569528/1‑99 (MQ=12)        gCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatgat        >  1:6854/1‑100 (MQ=33)        gCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatgat        <  2:1318311/100‑1 (MQ=33)         cGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGatgatg       >  1:373867/1‑100 (MQ=33)         cGAAGTTTTGCTCTTGAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGatgatg       <  1:2353794/100‑1 (MQ=18)         cGAAATTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGatgatg       <  2:754530/100‑1 (MQ=17)             gTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATc   >  2:799617/1‑100 (MQ=31)             gTTTTGCTCTGTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGCGTGTCTCAAATTTTCGCAACACGATGATGAATc   <  2:1566906/100‑1 (MQ=255)              ttttGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCg  <  1:1065719/100‑1 (MQ=11)                                                                    |                                              GTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCG  >  NC_000913/4037383‑4037495

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.