New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4037450	NA (NA)	56 (0.450)	22/184	10.4	25.3%	intergenic (+115/‑69)	alaT/rrlA	tRNA‑Ala/23S ribosomal RNA of rrnA operon
	?	NC_000913	= 4168568	178 (1.320)					intergenic (+121/‑73)	gltT/rrlB	tRNA‑Glu/23S ribosomal RNA of rrnB operon
Rejected: Coverage evenness skew score above cutoff.

GTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATC  >  NC_000913/4037383‑4037494                                                                    |                                             gTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAaca              <  2:1674774/100‑1 (MQ=11) gTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAaca              <  1:118140/100‑1 (MQ=11)  ttCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAAGACTGAACAACGAAAGTTGTTCGAGAGTCTCTGAAATTTTCGCAaca              >  2:2653243/1‑99 (MQ=2)  ntCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAacac             >  2:62077/2‑100 (MQ=14)    cACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACga           <  1:876438/100‑1 (MQ=33)    cACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACga           <  1:876454/100‑1 (MQ=33)    cACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACga           >  2:497261/1‑100 (MQ=33)     aCTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgat          >  2:2039800/1‑100 (MQ=33)      cTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatg         <  1:1961057/100‑1 (MQ=33)      cTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatg         >  2:1197942/1‑100 (MQ=33)      cTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatg         >  2:2334410/1‑100 (MQ=33)        gCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatgat       >  2:1864559/1‑100 (MQ=32)        gCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatgat       >  2:191275/1‑100 (MQ=33)         cGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGCAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGatgatg      <  2:2237344/100‑1 (MQ=18)         cGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGatgatg      <  2:2211096/100‑1 (MQ=33)         cGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGatgatg      <  1:1641026/100‑1 (MQ=33)             gTTTTGCTCCTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTCAGTCTCTCAAATTTTCGCAACACCATGATGAATc  <  2:126631/100‑1 (MQ=255)                                                                    |                                             GTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATC  >  NC_000913/4037383‑4037494

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