New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	W3110S.gb	2018571 =	0 (0.000)	8 (0.450)	5/96	NT	100%	coding (567/687 nt)	fliH	flagellar biosynthesis protein
	?	W3110S.gb	2519723 =	NA (NA)					intergenic (+33/‑54)	nupC/insL	nucleoside (except guanosine) transporter/predicted transposase

CTGCAACGTGTGGATGATATGCTCGGCGCTACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCCTCCATCCTGGCGGCTGTAAAGTCTCCGCCGATGAAGGCGATCTCGACGCCAGTGTCG  >  W3110S.gb/2018512‑2018641                                                            |                                                                       cTGCAACGTGTGGATGATATGCTCGGCGCTACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCAC                                                             >  1:1563882‑M1/1‑71 (MQ=255) cTGCAACGTGTGGATGATATGCTCGGCGCTACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCAC                                                             >  1:1801896‑M1/1‑71 (MQ=255)   gCAACGTGTGGATGATATGCTCGGCGCTACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCa                                                           >  1:412280‑M1/1‑70 (MQ=255)         cacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatGGGCGATCCCACCCTCCATc                                                     >  1:136083‑M1/52‑71 (MQ=255)           cgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatGGGCGATCCCACCCTCCATCCt                                                   <  1:999411‑M1/22‑1 (MQ=255)           cgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatGGGCGATCCCACCCTCCATCCt                                                   <  1:1734729‑M1/22‑1 (MQ=255)                    tGCTCGGCGCTACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataag                                                       >  1:1633248‑M1/1‑53 (MQ=255)                    tGCTCGGCGCTACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataag                                                       >  1:478380‑M1/1‑53 (MQ=255)                       tcaggcgtaataccacaacccttaagttagcgcttatGGGCGATCCCACCCTCCATCCTGGCGGCTGTaaa                                       >  1:1290126‑M1/38‑71 (MQ=255)                       tcaggcgtaataccacaacccttaagttagcgcttatGGGCGATCCCACCCTCCATCCTGGCGGCTGTaaa                                       >  1:940162‑M1/38‑71 (MQ=255)                       tcaggcgtaataccacaacccttaagttagcgcttatGGGCGATCCCACCCTCCATCCTGGCGGCTGTaaa                                       >  1:322383‑M1/38‑71 (MQ=255)                                 taccacaacccttaagttagcgcttatGGGCGATCCCACCCTCCATCCTGGCGGCTGTAAAGtc                                    >  1:1677473‑M1/28‑64 (MQ=255)                                   tAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatctg                             >  1:1425883‑M1/1‑38 (MQ=255)                                   tAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatctg                             >  1:188248‑M1/1‑38 (MQ=255)                                       aacccttaagttagcgcttatGGGCGATCCCACCCTCCATCCTGGCGGCTGTAAAGTCTCCGCCGATGAAg                       >  1:55743‑M1/22‑71 (MQ=255)                                                   gcgcTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcg            <  1:1757564‑M1/70‑49 (MQ=255)                                                   gcgcTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcg            <  1:358213‑M1/70‑49 (MQ=255)                                                      cTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaacca                                 <  1:1044090‑M1/46‑28 (MQ=255)                                                        tGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcct      >  1:335379‑M1/1‑17 (MQ=255)                                                        tGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcct      >  1:803037‑M1/1‑17 (MQ=255)                                                        tGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatctgaagaatgc                     <  1:141955‑M1/56‑40 (MQ=255)                                                         gCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctc     >  1:1179224‑M1/1‑16 (MQ=255)                                                            GGGCGATCCCACCCTCCATCCTGGCGGCTGTAAAGTCTCCGCCGATGAAGGCGATCTCGACGCCAGTGTCg  <  1:1409028‑M1/71‑1 (MQ=255)                                                            |                                                                       CTGCAACGTGTGGATGATATGCTCGGCGCTACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCCTCCATCCTGGCGGCTGTAAAGTCTCCGCCGATGAAGGCGATCTCGACGCCAGTGTCG  >  W3110S.gb/2018512‑2018641

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