New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? W3110S.gb = 20185830 (0.000)32 (0.870) 22/96 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

GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACC  >  W3110S.gb/2520997‑2521126
                                                              |                                                                   
gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGCTGGTTCAACCCTTAAGTTAGCGCTTATGGGATAAt                                                             >  1:98450/1‑71 (MQ=21)
     aaCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATa                                                               >  1:1994411/1‑64 (MQ=31)
     aaCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTc                                                           <  1:908363/68‑1 (MQ=39)
        gAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGt                                                     <  1:3363964/71‑1 (MQ=255)
          ggCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCATTAAGTTAGCGCTTATGGGATAATTCCCCGGttt                                                   >  1:1355297/1‑71 (MQ=255)
              ccGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATaa                                                              <  1:1547204/56‑1 (MQ=33)
                gAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATaa                                                              >  1:794717/1‑54 (MQ=33)
                gAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATaa                                                              >  1:3207871/1‑54 (MQ=33)
                   gCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGcc                                          >  1:258245/1‑71 (MQ=255)
                        gCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGtt                                      >  1:3440093/1‑70 (MQ=255)
                        gCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGtt                                      >  1:262097/1‑70 (MQ=255)
                         cATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTa                                     <  1:1429272/70‑1 (MQ=255)
                           ttcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                                  <  1:727357/71‑1 (MQ=255)
                            tcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                                <  1:1869024/56‑1 (MQ=255)
                             cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGc                                           >  1:2968466/1‑60 (MQ=255)
                             cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGc                                           >  1:1458953/1‑60 (MQ=255)
                             cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAg                                >  1:2985846/1‑71 (MQ=255)
                                cAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt                              <  1:2208349/70‑1 (MQ=255)
                                 aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTa                             <  1:1944478/70‑1 (MQ=255)
                                               ttAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                                <  1:2113559/37‑1 (MQ=255)
                                               ttAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGc                 <  1:2539339/68‑1 (MQ=255)
                                               ttAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGAc              >  1:2020781/1‑71 (MQ=255)
                                                tAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                                <  1:2202949/36‑1 (MQ=255)
                                                  aGTTAGCGCTTATGGGATAATTGCCCGGTTTTTTTACGCCTGTTaa                                    >  1:876473/1‑46 (MQ=255)
                                                      aGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCa        >  1:2666443/1‑70 (MQ=255)
                                                      aGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCa        <  1:3163496/70‑1 (MQ=255)
                                                      aGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCa        <  1:3617718/70‑1 (MQ=255)
                                                           ttATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATAcc  >  1:3385964/1‑71 (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.