New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? W3110S.gb = 20185835 (0.100)41 (0.860) 27/94 NT 89.5% coding (579/687 nt) fliH flagellar biosynthesis protein
?W3110S.gb = 2521059 NA (NA)intergenic (+170/+30) insL/yfeA predicted transposase/predicted diguanylate cyclase

GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATC  >  W3110S.gb/2520997‑2521119
                                                              |                                                            
gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAAt                                                      >  1:649460/1‑71 (MQ=34)
     aaCGAGGCACCGAGGCTTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTcccc                                                 >  1:2051241/1‑71 (MQ=255)
     aaCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTcccc                                                 >  1:701132/1‑71 (MQ=39)
          ggCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGttt                                            >  1:4324219/1‑71 (MQ=255)
             aCCGAGGCGTCGCATTCTTCATATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                         <  1:2696413/71‑1 (MQ=255)
              ccGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATaa                                                       >  1:277008/1‑56 (MQ=33)
              ccGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttt                                           <  1:1426646/68‑1 (MQ=255)
                gAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTAc                                      >  1:3969586/1‑71 (MQ=255)
                     gTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTccc                                                  <  1:2858225/54‑1 (MQ=255)
                      tCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTcccc                                                 <  1:4076045/54‑1 (MQ=255)
                      tCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGt                                <  1:1584362/71‑1 (MQ=255)
                      tCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGt                                <  1:3543982/71‑1 (MQ=255)
                           ttcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                         <  1:4550229/57‑1 (MQ=255)
                             cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                         >  1:3126752/1‑55 (MQ=255)
                             cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAg                         >  1:95407/1‑71 (MQ=255)
                             cttcAGATGGTTCAAACCTTAAGTTAGCGCTTATGGGATa                                                        >  1:3214148/1‑40 (MQ=255)
                               tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTcc                                                   >  1:115892/1‑43 (MQ=255)
                               tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAg                         >  1:110213/1‑69 (MQ=255)
                                cAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt                       <  1:3727945/70‑1 (MQ=255)
                                  gATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAAt                    <  1:582428/71‑1 (MQ=255)
                                     ggTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                           <  1:631233/61‑1 (MQ=255)
                                              cTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                         <  1:4475096/38‑1 (MQ=255)
                                               ttAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGc              <  1:4473446/64‑1 (MQ=255)
                                                tAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                           <  1:574986/50‑1 (MQ=255)
                                                     tAGCGCTTATGCGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATc  >  1:3866903/1‑70 (MQ=255)
                                                              |                                                            
GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATC  >  W3110S.gb/2520997‑2521119

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