New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	W3110S.gb	943798 =	94 (1.150)	7 (0.090)	3/96	NT	7.8%	coding (2418/2445 nt)	dmsA	dimethyl sulfoxide reductase, anaerobic, subunit A
	?	W3110S.gb	943824 =	73 (0.920)					coding (2444/2445 nt)	dmsA	dimethyl sulfoxide reductase, anaerobic, subunit A

GGAATCAATAAAAAATCCATACTGGGTTGTCATCGGTTACTCCTTACACCTTTTCAACCTGAACAAGGTTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  W3110S.gb/943868‑943798 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tAAGGAGTAACCGATGACAACCCAGTATGGATTTTTTATTGATTCCAGCCGTTGCACCGGTTGCAAAACCT  >  W3110S.gb/943824‑943893                                                                                                                                               GGAATCAATAAAAAATCCATACTGGGTTGTCATCGGTTACTCCTTACACCTTTTCAACCTGAACAAGGTTT                                                                        >  1:4243023/1‑71                                           CTTACACCTTTTCAACCTGAACAAGGTTTAAGGAGTAACCGATGACAACCCAGTATGGATTTTT                                     <  1:743192/64‑1                                           CTTACACCTTTTCAACCTGAACAAGGTTTAAGGAGTAACCGATGACAACCCAGTAT                                             <  1:1055110/56‑1                                             TACACCTTTTCAACCTGAACAAGGTTTAAGGAGTAACCGATGACAACCCAGTATGGATTTTTTA                                   >  1:1072198/1‑64                                             TACACCTTTTCAACCTGAACAAGGTTTAAGGAGTAACCGATGACAACCCAGTATGGATTTTTTA                                   >  1:1766222/1‑64                                             TACACCTTTTCAACCTGAACAAGGTTTAAGGAGTAACCGATGACAACCCAGTATGGATTTTTTA                                   >  1:3016622/1‑64                                             TACACCTTTTCAACCTGAACAAGGTTTAAGGAGTAACCGATGACAACCCAGTATGGATTTTTTA                                   >  1:4729360/1‑64                                             TACACCTTTTCAACCTGAACAAGGTTTAAGGAGTAACCGATGACAACCCAGTATGG                                           >  1:1248471/1‑56                                                                       TAAGGAGTAACCGATGACAACCCAGTATGGATTTTTTATTGATTCCAGCCGTTGCACCGGTTGCAAAACCT  <  1:1738704/71‑1                                                                       TAAGGAGTAACCGATGACAACCCAGTATGGATTTTTTATTGATTCCAGCCGTTGCACCGGTTGCAAAACC   <  1:3084524/70‑1                                                                                                                                               GGAATCAATAAAAAATCCATACTGGGTTGTCATCGGTTACTCCTTACACCTTTTCAACCTGAACAAGGTTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  W3110S.gb/943868‑943798 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tAAGGAGTAACCGATGACAACCCAGTATGGATTTTTTATTGATTCCAGCCGTTGCACCGGTTGCAAAACCT  >  W3110S.gb/943824‑943893

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.