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 36 [0] [0] 32 [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)32 (1.070) 27/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

TCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTT  >  NC_000913/2137669‑2137765
                                                 |                                               
attaaagtataaatagcttatccatgcttatatgcttacggctttatatg                                                 >  1:1708622‑M2/50‑50 (MQ=255)
attaaagtataaatagcttatccatgcttatatgcttacggctttatatg                                                 >  1:1731614‑M2/50‑50 (MQ=255)
  taaagtataaatagcttatccatgcttatatgcttacggctttatatGtc                                               >  1:1603986‑M2/48‑50 (MQ=255)
     agtataaatagcttatccatgcttatatgcttacggctttatatGtctc                                             >  1:1350215‑M2/45‑49 (MQ=255)
       tataaatagcttatccatgcttatatgcttacggctttatatGTCTTca                                           >  1:2298993‑M2/43‑49 (MQ=255)
       tataaatagcttatccatgcttatatgcttacggctttatatGTCTCcac                                          <  1:2115989‑M2/8‑1 (MQ=255)
        ataaatagcttatccatgcttatatgcttacggctttatatGTCTCcaca                                         >  1:2427921‑M2/42‑50 (MQ=255)
         taaatagcttatccatgcttatatgcttacggctttatatGTCTCCACAt                                        <  1:876300‑M2/10‑1 (MQ=255)
           aatagcttatccatgcttatatgcttacggctttatatGTCTCCACATg                                       <  1:605038‑M2/11‑1 (MQ=255)
            atagcttatccatgcttatatgcttacggctttatatGTCTCCACATGc                                      >  1:1230130‑M2/38‑49 (MQ=255)
                  tatccatgcttatatgcttacggctttatatGTCTCCACATGCGATAttt                               >  1:348678‑M2/32‑50 (MQ=255)
                   atccatgcttatatgcttacggctttatatGTCTCCACATGCGATAtt                                <  1:1771537‑M2/18‑1 (MQ=255)
                    tccatgcttatatgcttacggctttatatGTCTCCACATGCGATATTTCt                             <  1:2154333‑M2/21‑1 (MQ=255)
                       atgcttatatgcttacggctttatatGTCTCCACATGCGATATTTCTTaa                          <  1:961359‑M2/24‑1 (MQ=255)
                         gcttatatgcttacggctttatatGTCTCCACATGCGATATTTCTTaaa                         >  1:852799‑M2/25‑49 (MQ=255)
                         gcttatatgcttacggctttatatGTCTCCACATGCGATATTTCTTaaa                         >  1:816894‑M2/25‑49 (MQ=255)
                          cttatatgcttacggctttatatGTCTCCACATGCGATATTTCTTAaata                       <  1:2067857‑M2/27‑1 (MQ=255)
                          cttatatgcttacggctttatatGTCTCCACATGCGATATTTCTTAaat                        <  1:749222‑M2/26‑1 (MQ=255)
                          cttatatgcttacggctttatatGTCTCCACATGCGATATTTCTTAaat                        <  1:1455089‑M2/26‑1 (MQ=255)
                            tatatgcttacggctttatatGTCTCCACATGCGATATTTCTTAaataat                     >  1:2402617‑M2/22‑50 (MQ=255)
                              tatgcttacggctttatatGTCTCCACATGCGATATTTCTTAAATAATg                    <  1:2153172‑M2/30‑1 (MQ=255)
                                tgcttacggctttatatGTCTCCACATGCGATATTTCTTAAATAATGtt                  <  1:1399269‑M2/32‑1 (MQ=255)
                                 gcttacggctttatatGTCTCCACATGCGATATTTCTTAAATAATGttt                 >  1:2250774‑M2/17‑49 (MQ=255)
                                  cttacggctttatatGTCTCCACATGCGATATTTCTTAAATAATGtttt                <  1:340685‑M2/34‑1 (MQ=255)
                                  cttacggctttatatGTCTCCACATGCGATATTTCTTAAATAATGTTtta               <  1:2632818‑M2/35‑1 (MQ=255)
                                     acggctttatatGTCTCCACATGCGATATTTCTTAAATAATGTTttatt             <  1:432935‑M2/37‑1 (MQ=255)
                                     acggctttatatGTCTCCACATGCGATATTTCTTAAATAATGTTttatt             <  1:52991‑M2/37‑1 (MQ=255)
                                     acggctttatatGTCTCCACATGCGATATTTCTTAAATAATGTTttatt             <  1:108597‑M2/37‑1 (MQ=255)
                                      cggctttatatGTCTCCACATGCGATATTTCTTAAATAATGTTttattat           >  1:1926375‑M2/12‑50 (MQ=255)
                                       ggctttatatGTCTCCACATGCGATATTTCTTAAATAATGTTttattatt          <  1:2009963‑M2/40‑1 (MQ=255)
                                       ggctttatatGTCTCCACATGCGATATTTCTTAAATAATGTTttattat           >  1:1671341‑M2/11‑49 (MQ=255)
                                               atGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACttt  >  1:1125123‑M2/3‑50 (MQ=255)
                                                 |                                               
TCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTT  >  NC_000913/2137669‑2137765

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: