New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? W3110S.gb = 20185831 (0.020)38 (0.830) 30/98 NT 97.5% coding (579/687 nt) fliH flagellar biosynthesis protein
?W3110S.gb = 2521059 NA (NA)intergenic (+170/+30) insL/yfeA predicted transposase/predicted diguanylate cyclase

GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCT  >  W3110S.gb/2520997‑2521128
                                                              |                                                                     
gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATaa                                                                <  1:3614744/70‑1 (MQ=33)
 tCTTAACGAGGGACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAAtt                                                              <  1:3981248/71‑1 (MQ=25)
     aaCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTc                                                             <  1:345568/68‑1 (MQ=39)
             aCCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                                  <  1:250335/71‑1 (MQ=255)
                gAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTAc                                               >  1:2485511/1‑71 (MQ=255)
                        gCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTa                                       <  1:142260/71‑1 (MQ=255)
                        gCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTa                                       <  1:3314778/71‑1 (MQ=255)
                        gCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTa                                       <  1:1635837/71‑1 (MQ=255)
                            tcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                                    <  1:4011606/70‑1 (MQ=255)
                             cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACga                                             >  1:2212609/1‑59 (MQ=255)
                              ttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATaa                                                                >  1:3323293/1‑40 (MQ=37)
                               tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGc                                 >  1:175621/1‑70 (MQ=255)
                                cAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt                                <  1:1922729/70‑1 (MQ=255)
                                cAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt                                <  1:3628176/70‑1 (MQ=255)
                                       ttCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGc                          >  1:417713/1‑69 (MQ=255)
                                       ttCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCt                         >  1:2167245/1‑70 (MQ=255)
                                       ttCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTg                        >  1:354983/1‑71 (MQ=255)
                                        tCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGc                       >  1:1751726/1‑71 (MQ=255)
                                               ttAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGc                       >  1:116021/1‑64 (MQ=255)
                                                  aGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGACTGTTAATCAGCTAATGGCTGCGGGCGACCaa             >  1:1706873/1‑71 (MQ=255)
                                                     tAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATc           <  1:69605/70‑1 (MQ=255)
                                                       gcgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCaaa        <  1:657128/71‑1 (MQ=255)
                                                             aTGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCt  >  1:4121513/1‑71 (MQ=255)
                                                             aTGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCt  >  1:3213965/1‑71 (MQ=255)
                                                              |                                                                     
GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCT  >  W3110S.gb/2520997‑2521128

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.