New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	W3110S.gb	2018571 =	4 (0.090)	45 (1.040)	26/90	NT	92.3%	coding (567/687 nt)	fliH	flagellar biosynthesis protein
	?	W3110S.gb	2519723 =	NA (NA)					intergenic (+33/‑54)	nupC/insL	nucleoside (except guanosine) transporter/predicted transposase

TGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTC  >  W3110S.gb/2519654‑2519787                                                                      |                                                                 tgtgCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCa                                                                  <  1:4615848/70‑1 (MQ=255)            cGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTa                                                         <  1:112864/68‑1 (MQ=255)                aaTCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAg                                                  <  1:1107113/71‑1 (MQ=255)                          cTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCAt                                                                 <  1:4379444/46‑1 (MQ=255)                             gTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                       <  1:693611/69‑1 (MQ=255)                             gTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                       <  1:2010920/69‑1 (MQ=255)                             gTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                       <  1:4372219/69‑1 (MQ=255)                                cTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATaa                                                               >  1:1438760/1‑42 (MQ=255)                                cTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCt                                                          <  1:687281/47‑1 (MQ=255)                                cTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACtt                                                     <  1:833690/52‑1 (MQ=255)                                 tGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGt                                         <  1:1524877/63‑1 (MQ=255)                                   tAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                >  1:4318803/1‑70 (MQ=255)                                       aCCATACATAAAAAAGCCGGGGATAATTCCCATAAgcg                                                            <  1:2692592/38‑1 (MQ=255)                                          atacataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                <  1:1587042/63‑1 (MQ=255)                                           tacataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGTTTGTGGTATTa                                     >  1:2174378/1‑57 (MQ=255)                                           tacataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                       <  1:927712/55‑1 (MQ=255)                                             cataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                       <  1:564617/53‑1 (MQ=255)                                             cataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                      >  1:2124802/1‑70 (MQ=255)                                             cataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                      <  1:667074/70‑1 (MQ=255)                                               taAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                      <  1:2571152/68‑1 (MQ=255)                                                aaaaaaGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                  >  1:1254066/1‑71 (MQ=255)                                                aaaaaaGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                  >  1:3936446/1‑71 (MQ=255)                                                 aaaaaGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                >  1:378955/1‑56 (MQ=255)                                                    aaGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGa                              >  1:3022984/1‑55 (MQ=255)                                                    aaGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGcc              <  1:2359131/71‑1 (MQ=255)                                                     aGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                     >  1:2734211/1‑47 (MQ=255)                                                         ggggATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGa         <  1:4041137/71‑1 (MQ=39)                                                           ggATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt       <  1:1764566/71‑1 (MQ=39)                                                             aTAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa     <  1:4463345/71‑1 (MQ=38)                                                              tAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                      <  1:2181386/53‑1 (MQ=38)                                                               aaTTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct   >  1:4227433/1‑71 (MQ=34)                                                                aTTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGa         >  1:2152319/1‑64 (MQ=33)                                                                aTTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt       <  1:1447771/66‑1 (MQ=33)                                                                aTTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActc  >  1:3965466/1‑71 (MQ=32)                                                                      |                                                                 TGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTC  >  W3110S.gb/2519654‑2519787

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.