New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	W3110S.gb	= 2018583	18 (0.640)	24 (0.880)	3/104	0.4	58.1%	coding (579/687 nt)	fliH	flagellar biosynthesis protein
	?	W3110S.gb	= 2521059	NA (NA)					intergenic (+170/+30)	insL/yfeA	predicted transposase/predicted diguanylate cyclase

CGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCG  >  W3110S.gb/2521017‑2521112                                           |                                                      cGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttt                                    >  1:529002/1‑62 (MQ=255) cGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttt                                    >  1:2863904/1‑62 (MQ=255) cGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttt                                    >  1:2835322/1‑62 (MQ=255) cGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttt                                    >  1:2285212/1‑62 (MQ=255)            tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                  >  1:1544505/1‑53 (MQ=255)            tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                  >  1:855594/1‑53 (MQ=255)            tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                  >  1:1570459/1‑53 (MQ=255)            tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCt                           >  1:187405/1‑60 (MQ=255)              aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTa                       >  1:1027153/1‑62 (MQ=255)                 tGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                    >  1:861818/1‑62 (MQ=255)                 tGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                    >  1:3092954/1‑62 (MQ=255)                 tGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                    >  1:2322871/1‑62 (MQ=255)                 tGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                    >  1:2282067/1‑62 (MQ=255)                 tGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                    >  1:2145671/1‑62 (MQ=255)                 tGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACACCTGTTAATc                    >  1:15494/1‑62 (MQ=255)                                    gcgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCg  >  1:26893/1‑61 (MQ=255)                                    gcgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCg  >  1:1952583/1‑61 (MQ=255)                                    gcgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCg  >  1:655045/1‑61 (MQ=255)                                    gcgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCg  >  1:812481/1‑61 (MQ=255)                                         tATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt                >  1:1789666/1‑42 (MQ=255)                                         tATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt                >  1:1124455/1‑42 (MQ=255)                                           |                                                      CGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCG  >  W3110S.gb/2521017‑2521112

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