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

GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACC  >  W3110S.gb/2520997‑2521126                                                               |                                                                    gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATa                                                               >  1:765700/1‑69 (MQ=31)    ttAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAAtt                                                            <  1:2268773/69‑1 (MQ=37)    ttAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAAtt                                                            <  1:1684542/69‑1 (MQ=37)       aCGAGGCACCGAGGCGTCGCATTGTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTccc                                                         >  1:2505815/1‑69 (MQ=255)        cGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTccc                                                         >  1:885754/1‑68 (MQ=40)        cGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTcc                                                          >  1:1177141/1‑67 (MQ=40)            gCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtt                                                    >  1:3018222/1‑69 (MQ=255)                cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                                <  1:1982585/69‑1 (MQ=255)                           aTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTa                                     <  1:1421655/69‑1 (MQ=255)                           aTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTa                                     >  1:2496933/1‑69 (MQ=255)                             tcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTa                                     >  1:1718734/1‑67 (MQ=255)                              cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTa                                              >  1:2593497/1‑57 (MQ=255)                                   gATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt                              <  1:2777130/68‑1 (MQ=255)                                   gATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt                              <  1:1031058/68‑1 (MQ=255)                                     tGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAAt                           <  1:1900595/69‑1 (MQ=255)                                     tGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAAt                           <  1:2862998/69‑1 (MQ=255)                                      ggTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAAt                           <  1:3115814/68‑1 (MQ=255)                                         tCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCt                       >  1:1650739/1‑69 (MQ=255)                                         tCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCt                       >  1:2291014/1‑69 (MQ=255)                                          cAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                                >  1:1153701/1‑43 (MQ=255)                                                 tAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGc                 <  1:3615255/67‑1 (MQ=255)                                                 tAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCg                >  1:3686241/1‑68 (MQ=255)                                                  aaGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGa               >  1:2627031/1‑68 (MQ=255)                                                  aaGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGa               >  1:734105/1‑68 (MQ=255)                                                      tAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGcc                                          >  1:152609/1‑37 (MQ=255)                                                       aGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGtt                                      >  1:3092055/1‑40 (MQ=255)                                                       aGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTa                                     >  1:2879026/1‑41 (MQ=255)                                                         cgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCa        <  1:2216675/68‑1 (MQ=255)                                                         cgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCa        <  1:3216276/68‑1 (MQ=255)                                                         cgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCa        <  1:2102541/68‑1 (MQ=255)                                                         cgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCa        <  1:1371990/68‑1 (MQ=255)                                                              aTGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATAcc  >  1:2875426/1‑69 (MQ=255)                                                              aTGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATAcc  >  1:2013661/1‑69 (MQ=255)                                                               tGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt                              >  1:562743/1‑40 (MQ=255)                                                               tGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt                              <  1:867564/40‑1 (MQ=255)                                                               |                                                                    GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACC  >  W3110S.gb/2520997‑2521126

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.