New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	W3110S.gb	= 2018583	1 (0.040)	13 (0.600)	12/94	NT	93.1%	coding (579/687 nt)	fliH	flagellar biosynthesis protein
	?	W3110S.gb	= 2521059	NA (NA)					intergenic (+170/+30)	insL/yfeA	predicted transposase/predicted diguanylate cyclase

ACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATC  >  W3110S.gb/2521003‑2521119                                                         |                                                             aCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCg                                                >  1:497269/1‑71 (MQ=39)     ggCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGttt                                            >  1:280446/1‑71 (MQ=255)       cACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGttttt                                          <  1:904994/71‑1 (MQ=255)        aCCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                         <  1:12061/71‑1 (MQ=255)        aCCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                         <  1:853245/71‑1 (MQ=255)                    cATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCt                                  <  1:1095858/66‑1 (MQ=255)                       tcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAAt                            >  1:2117792/1‑69 (MQ=255)                           cAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt                       <  1:571178/70‑1 (MQ=255)                             gATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAAt                    <  1:1108181/71‑1 (MQ=255)                                   tCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTa                                       >  1:2172667/1‑46 (MQ=255)                                            aaGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGAcc      >  1:378041/1‑70 (MQ=255)                                               ttAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATc  <  1:1122269/71‑1 (MQ=255)                                                  gcgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTaa                             >  1:1485499/1‑41 (MQ=255)                                                  gcgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTaa                             >  1:369031/1‑41 (MQ=255)                                                         |                                                             ACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATC  >  W3110S.gb/2521003‑2521119

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