New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? W3110S.gb = 20185830 (0.000)51 (1.340) 5/118 0.2 100% coding (579/687 nt) fliH flagellar biosynthesis protein
?W3110S.gb = 2521059 NA (NA)intergenic (+170/+30) insL/yfeA predicted transposase/predicted diguanylate cyclase

GCGCTACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCCTCCATCCTGGCGGCTGTAAAGTCTCCGCCGATGAAGGCGATCTCGACGCCAGTGTCGC  >  W3110S.gb/2018537‑2018642
                                              |                                                           
gcgcTACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagg                                           >  1:1889767‑M1/1‑47 (MQ=255)
gcgcTACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagg                                           >  1:1237338‑M1/1‑47 (MQ=255)
gcgcTACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagg                                           >  1:754514‑M1/1‑47 (MQ=255)
gcgcTACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagg                                           >  1:2429716‑M1/1‑47 (MQ=255)
gcgcTACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagg                                           >  1:2691465‑M1/1‑47 (MQ=255)
gcgcTACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagg                                           >  1:1329864‑M1/1‑47 (MQ=255)
    tACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttg                                       <  1:1913916‑M1/65‑23 (MQ=255)
         tAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaacca                                  >  1:3079567‑M1/1‑38 (MQ=255)
         tAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaacca                                  >  1:2218565‑M1/1‑38 (MQ=255)
         tAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaacca                                  >  1:2749437‑M1/1‑38 (MQ=255)
         tAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaacca                                  >  1:2952642‑M1/1‑38 (MQ=255)
         tAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaacca                                  >  1:2956051‑M1/1‑38 (MQ=255)
         tAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaacca                                  >  1:2353349‑M1/1‑38 (MQ=255)
         tAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaacca                                  >  1:400492‑M1/1‑38 (MQ=255)
         tAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaacca                                  >  1:498024‑M1/1‑38 (MQ=255)
         tAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaacca                                  >  1:607786‑M1/1‑38 (MQ=255)
         tAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaacca                                  >  1:820984‑M1/1‑38 (MQ=255)
         tAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaacca                                  >  1:99450‑M1/1‑38 (MQ=255)
         tAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaacca                                  >  1:162294‑M1/1‑38 (MQ=255)
         tAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaacca                                  >  1:1497409‑M1/1‑38 (MQ=255)
         tAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaacca                                  >  1:1252828‑M1/1‑38 (MQ=255)
         tAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaacca                                  >  1:1758033‑M1/1‑38 (MQ=255)
         tAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaacca                                  >  1:1250763‑M1/1‑38 (MQ=255)
         tAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaacca                                  >  1:1803417‑M1/1‑38 (MQ=255)
         tAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaacca                                  >  1:1835827‑M1/1‑38 (MQ=255)
             tttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:705308‑M1/64‑31 (MQ=255)
             tttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:1180329‑M1/64‑31 (MQ=255)
             tttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:603239‑M1/64‑31 (MQ=255)
             tttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:1173615‑M1/64‑31 (MQ=255)
             tttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:126919‑M1/64‑31 (MQ=255)
             tttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:1273994‑M1/64‑31 (MQ=255)
             tttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:1338476‑M1/64‑31 (MQ=255)
             tttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:1395909‑M1/64‑31 (MQ=255)
             tttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:1419828‑M1/64‑31 (MQ=255)
             tttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:2912267‑M1/64‑31 (MQ=255)
             tttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:2879746‑M1/64‑31 (MQ=255)
             tttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:2712168‑M1/64‑31 (MQ=255)
             tttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:2700856‑M1/64‑31 (MQ=255)
             tttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:2699268‑M1/64‑31 (MQ=255)
             tttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:1529411‑M1/64‑31 (MQ=255)
             tttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:2678165‑M1/64‑31 (MQ=255)
             tttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:2586220‑M1/64‑31 (MQ=255)
             tttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:2546116‑M1/64‑31 (MQ=255)
             tttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:2314490‑M1/64‑31 (MQ=255)
             tttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:2165911‑M1/64‑31 (MQ=255)
             tttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:1907644‑M1/64‑31 (MQ=255)
              ttGCATGGCTGGCGCTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:1903172‑M1/63‑31 (MQ=255)
                         gcgcTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatctgaagaatgcgacg                  >  1:1000970‑M1/1‑22 (MQ=255)
                         gcgcTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatctgaagaatgcgacg                  >  1:2145264‑M1/1‑22 (MQ=255)
                         gcgcTTGCGGGGCGATCCCACCataagcgctaacttaagggttgaaccatctgaagaatgcgacg                  >  1:759785‑M1/1‑22 (MQ=255)
                                   ggCGATCCCACCataagcgctaacttaagggttgaaccatct                               <  1:688294‑M1/42‑31 (MQ=255)
                                             ccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgtta  <  1:2874102‑M1/61‑60 (MQ=255)
                                              |                                                           
GCGCTACCTTAAGTTTGCATGGCTGGCGCTTGCGGGGCGATCCCACCCTCCATCCTGGCGGCTGTAAAGTCTCCGCCGATGAAGGCGATCTCGACGCCAGTGTCGC  >  W3110S.gb/2018537‑2018642

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