Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,101,397 Δ36,321 bp [wbbL]wza 34 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2101397 2137717 36321 34 [0] [0] 34 [wbbL]–wza [wbbL],insH‑7,wbbL,wbbK,wbbJ,wbbI,wbbH,glf,rfbX,rfbC,rfbA,rfbD,rfbB,galF,wcaM,wcaL,wcaK,wzxC,wcaJ,cpsG,cpsB,wcaI,gmm,fcl,gmd,wcaF,wcaE,wcaD,wcaC,wcaB,wcaA,wzc,wzb,wza

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 21013960 (0.000)34 (1.270) 29/96 0.0 100% pseudogene (349/349 nt) wbbL interrupted rhamnosyltransferase WbbL
?NC_000913 2137718 = 0 (0.000)intergenic (‑475/‑184) wza/yegH outer membrane polysaccharide export protein Wza/inner membrane protein YegH

ATTAAAGTATAAATAGCTTATCCATGCTTATATGCTTACGGCTTTATATTACGGGTGAAAAACTGATGAAATTCGATCAAAGTTGCGATTTGAT  >  NC_000913/2101348‑2101441
                                                |                                             
attaAAGTATAAATAGCTTATCCATGCTTATATGCTTACGGCTTTATATg                                              >  1:1001627‑M1/1‑49 (MQ=255)
    aaGTATAAATAGCTTATCCATGCTTATATGCTTACGGCTTTATATgtct                                           <  1:2251556‑M1/49‑5 (MQ=255)
    aaGTATAAATAGCTTATCCATGCTTATATGCTTACGGCTTTATATgtct                                           <  1:746497‑M1/49‑5 (MQ=255)
      gTATAAATAGCTTATCCATGCTTATATGCTTACGGCTTTATATgtctcca                                        <  1:2206933‑M1/50‑8 (MQ=255)
       tataAATAGCTTATCCATGCTTATATGCTTACGGCTTTATATgtctcca                                        >  1:1908405‑M1/1‑42 (MQ=255)
            aTAGCTTATCCATGCTTATATGCTTACGGCTTTATATgtctccacatgcg                                  <  1:803404‑M1/50‑14 (MQ=255)
               gCTTATCCATGCTTATATGCTTACGGCTTTATATgtctccacatgcgata                               <  1:1773578‑M1/50‑17 (MQ=255)
                 ttATCCATGCTTATATGCTTACGGCTTTATATgtctccacatgcgatatt                             >  1:811793‑M1/1‑32 (MQ=255)
                  tATCCATGCTTATATGCTTACGGCTTTATATgtctccacatgcgatattt                            >  1:701594‑M1/1‑31 (MQ=255)
                   aTCCATGCTTATATGCTTACGGCTTTATATgtctccacatgcgatatt                             <  1:2175830‑M1/48‑19 (MQ=255)
                     ccATGCTTATATGCTTACGGCTTTATATgtctccacatgcgatatttct                          >  1:1399230‑M1/1‑28 (MQ=255)
                      cATGCTTATATGCTTACGGCTTTATATgtctccacatgcgatatttctt                         >  1:1024163‑M1/1‑27 (MQ=255)
                       aTGCTTATATGCTTACGGCTTTATATgtctccacatgcgatatttctta                        >  1:697542‑M1/1‑26 (MQ=255)
                         gCTTATATGCTTACGGCTTTATATgtctccacatgcgatatttcttaaa                      >  1:173521‑M1/1‑24 (MQ=255)
                         gCTTATATGCTTACGGCTTTATATgtctccacatgcgatatttcttaaa                      >  1:1712532‑M1/1‑24 (MQ=255)
                          cTTATATGCTTACGGCTTTATATgtctccacatgcgatatttcttaaata                    <  1:1399773‑M1/50‑28 (MQ=255)
                          cTTATATGCTTACGGCTTTATATgtctccacatgcgatatttcttaaat                     >  1:738096‑M1/1‑23 (MQ=255)
                           ttATATGCTTACGGCTTTATATgtctccacatgcgatatttcttaaataa                   >  1:246704‑M1/1‑22 (MQ=255)
                                tGCTTACGGCTTTATATgtctccacatgcgatatttcttaaataatgttt              <  1:96138‑M1/50‑34 (MQ=255)
                                tGCTTACGGCTTTATATgtctccacatgcgatatttcttaaataatgttt              >  1:351228‑M1/1‑17 (MQ=255)
                                tGCTTACGGCTTTATATgtctccacatgcgatatttcttaaataatgtt               <  1:658794‑M1/49‑33 (MQ=255)
                                  cTTACGGCTTTATATgtctccacatgcgatatttcttaaataatgtttt             <  1:255595‑M1/49‑35 (MQ=255)
                                    tACGGCTTTATATgtctccacatgcgatatttcttaaataatgttttat           >  1:1922980‑M1/1‑13 (MQ=255)
                                    tACGGCTTTATATgtctccacatgcgatatttcttaaataatgttttat           >  1:986170‑M1/1‑13 (MQ=255)
                                     aCGGCTTTATATgtctccacatgcgatatttcttaaataatgttttatta         >  1:627727‑M1/1‑12 (MQ=255)
                                     aCGGCTTTATATgtctccacatgcgatatttcttaaataatgttttat           <  1:868009‑M1/48‑37 (MQ=255)
                                      cGGCTTTATATgtctccacatgcgatatttcttaaataatgttttattat        >  1:245249‑M1/1‑11 (MQ=255)
                                      cGGCTTTATATgtctccacatgcgatatttcttaaataatgttttattat        <  1:1379675‑M1/50‑40 (MQ=255)
                                      cGGCTTTATATgtctccacatgcgatatttcttaaataatgttttatta         <  1:26075‑M1/49‑39 (MQ=255)
                                      cGGCTTTATATgtctccacatgcgatatttcttaaataatgttttatta         <  1:868637‑M1/49‑39 (MQ=255)
                                      cGGCTTTATATgtctccacatgcgatatttcttaaataatgttttatta         <  1:975299‑M1/49‑39 (MQ=255)
                                       ggCTTTATATgtctccacatgcgatatttcttaaataatgttttattat        >  1:1844864‑M1/1‑10 (MQ=255)
                                         cTTTATATgtctccacatgcgatatttcttaaataatgttttattattac     <  1:48313‑M1/50‑43 (MQ=255)
                                             atatgtctccacatgcgatatttcttaaataatgttttattattaccac  <  1:1359264‑M1/49‑46 (MQ=255)
                                                |                                             
ATTAAAGTATAAATAGCTTATCCATGCTTATATGCTTACGGCTTTATATTACGGGTGAAAAACTGATGAAATTCGATCAAAGTTGCGATTTGAT  >  NC_000913/2101348‑2101441

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