New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	W3110S.gb	= 2018583	0 (0.000)	56 (1.040)	42/98	NT	100%	coding (579/687 nt)	fliH	flagellar biosynthesis protein
	?	W3110S.gb	= 2521059	NA (NA)					intergenic (+170/+30)	insL/yfeA	predicted transposase/predicted diguanylate cyclase

CGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTT  >  W3110S.gb/2520996‑2521129                                                                |                                                                       cGTCTTAACGAGGCACCGAGGCGTCGCATTCTTAAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATaa                                                                 >  1:321201/1‑71 (MQ=17)  gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCATATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAAt                                                                >  1:4657040/1‑71 (MQ=21)  gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTTAACCCTTAAGTTAGCGCTTATGGGATaa                                                                 <  1:1345909/70‑1 (MQ=18)       aaCTAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTc                                                              <  1:2316856/68‑1 (MQ=255)         cGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCgg                                                         <  1:5154431/71‑1 (MQ=40)           aGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtt                                                       >  1:681260/1‑71 (MQ=255)              cACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGttttt                                                    <  1:595847/71‑1 (MQ=255)               aCCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                                   <  1:2204995/71‑1 (MQ=255)               aCCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                                   >  1:3980504/1‑71 (MQ=255)                 cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTa                                                 <  1:1675342/71‑1 (MQ=255)                          gCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGtt                                         >  1:4602363/1‑70 (MQ=255)                          gCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGtt                                         >  1:4106613/1‑70 (MQ=255)                          gCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGt                                          >  1:4119283/1‑69 (MQ=255)                             ttcttcAGATGGTTCAAGCCTTAAGTTAGCGCTTATGGGATAATTgcc                                                            <  1:4915808/48‑1 (MQ=14)                             ttcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAAt                                      >  1:4382221/1‑70 (MQ=255)                             ttcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                                     <  1:2815989/71‑1 (MQ=255)                              tcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAAt                                      >  1:386329/1‑69 (MQ=255)                                 tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGc                                  >  1:977777/1‑70 (MQ=255)                                   aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTa                                <  1:1118504/70‑1 (MQ=255)                                   aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTa                                <  1:1230159/70‑1 (MQ=255)                                     aTGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCg                                                          <  1:4134200/42‑1 (MQ=255)                                     aTGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATg                             <  1:1159990/71‑1 (MQ=255)                                       ggTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATgg                            >  1:2715155/1‑70 (MQ=255)                                         ttCAACCCTTAAGTTAGCGCTTATGGGATAATTccc                                                            <  1:939676/36‑1 (MQ=255)                                         ttCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGc                           >  1:1462602/1‑69 (MQ=255)                                           cAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGc                        <  1:4962501/70‑1 (MQ=255)                                           cAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGc                        <  1:167988/70‑1 (MQ=255)                                            aaCCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttt                                                     <  1:4757555/40‑1 (MQ=255)                                            aaCCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCgg                      <  1:4697110/71‑1 (MQ=255)                                               ccTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                                   <  1:4122608/39‑1 (MQ=255)                                                      ttAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATc            <  1:3055808/71‑1 (MQ=255)                                                          cgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt                                 <  1:1334146/46‑1 (MQ=255)                                                            cTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATAc      >  1:2856760/1‑71 (MQ=255)                                                               aTGGGATAATTCCCCGGTTTTTTTACGCCTGTTAAt                                      >  1:2258185/1‑36 (MQ=255)                                                                tGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCtt  <  1:1948507/71‑1 (MQ=255)                                                                tGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCtt  <  1:2373717/71‑1 (MQ=255)                                                                |                                                                       CGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTT  >  W3110S.gb/2520996‑2521129

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