New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3943624	NA (NA)	84 (0.460)	26/182	6.2	77.7%	intergenic (+114/‑80)	gltU/rrlC	tRNA‑Glu/23S ribosomal RNA of rrnC operon
	?	NC_000913	= 4037457	26 (0.130)					intergenic (+122/‑62)	alaT/rrlA	tRNA‑Ala/23S ribosomal RNA of rrnA operon
Rejected: Coverage evenness skew score above cutoff.

TCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATC  >  NC_000913/4037385‑4037494                                                                         |                                      tCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAacac             <  1:2202274/99‑1 (MQ=32) tCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAacac             >  1:3350096/1‑99 (MQ=32)   aCTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACga           >  1:1564106/1‑99 (MQ=34)   aCTGCGAAGTTTTGCTCTTTAAAAAGCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCTTGAGTCTCTCGAATTTTCGCAACACga           <  2:3537797/99‑1 (MQ=255)     tGCGAAGTTTTGCTCTTTGAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatg         <  2:4289289/99‑1 (MQ=18)     tGCGAAGTTTTGCTCTTTAAAAATCTGGGTCAAGCTTAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGCGTCTCTCAAATTTTCGCAACACgatg         >  2:1029445/1‑99 (MQ=255)     tGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCACCACgatg         >  2:3409927/1‑99 (MQ=18)     tGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatg         >  2:3410905/1‑99 (MQ=34)     tGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatg         >  1:1094417/1‑99 (MQ=33)     tGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatg         >  2:1821279/1‑99 (MQ=34)      gCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGCAACACTGAACAACGAAAGTTGTTCGTTAGTCTCTCAAATTTTCGCATCGCgatga        <  2:4261339/99‑1 (MQ=11)      gCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatga        >  2:307176/1‑99 (MQ=34)       cGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACATCGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatgat       <  2:2845769/99‑1 (MQ=18)       cGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAACTTTTCGCAACACgatgat       >  2:4869169/1‑99 (MQ=18)       cGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatgat       >  2:1826554/1‑99 (MQ=34)       cGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGAAACACgatgat       <  1:4824888/99‑1 (MQ=21)        gAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGatgatg      >  1:2227061/1‑99 (MQ=33)        gAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGatgatg      >  1:3804557/1‑99 (MQ=34)        gAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGatgatg      <  2:2030121/99‑1 (MQ=34)         aaGTTTTGCTCTTTGAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGAtgatga     <  2:1547668/99‑1 (MQ=21)         aaGTTTTGCTCTTTAAAAATCTGGATCCAGCTGAAAATTGAAACAATGAACAACGAAAGTTGTGGGTGAGTCTATCAAATTTTCGCAACACGAtgatga     <  1:5018106/99‑1 (MQ=11)          aGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGaa    <  1:4448598/99‑1 (MQ=33)            ttttGCTCTTTAAAAATCTGGATCGAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATc  <  1:5161369/99‑1 (MQ=14)            ttttGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATc  <  1:4992943/99‑1 (MQ=31)            ttttGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATc  >  1:2429289/1‑99 (MQ=31)            ttttGCTCTTTAAAAATCTGCATCAAGCTGAAAATTGAAACACTCAACAACGAAAGTTGTTCCTGAGTCTCTCAAATTTTCGCAACACGATGATGAATc  <  2:847567/99‑1 (MQ=255)                                                                         |                                      TCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATC  >  NC_000913/4037385‑4037494

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