Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,978,503 Δ776 bp insB‑5insA‑5 insB‑5, insA‑5

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1978503–1979269 1979278 10–776 85 [0] [0] 80 insB‑5–insA‑5 insB‑5,insA‑5

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 19785020 (0.000)68 (1.450) 28/118 0.1 100% intergenic (‑305/+16) flhD/insB‑5 DNA‑binding transcriptional dual regulator FlhD/IS1 protein InsB
?NC_000913 1979279 = 0 (0.000)intergenic (‑64/‑474) insA‑5/uspC IS1 protein InsA/universal stress protein C

GGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAG  >  NC_000913/1978435‑1978568
                                                                   |                                                                  
ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                    >  1:1184848‑M1/1‑68 (MQ=255)
ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                    >  1:1262701‑M1/1‑68 (MQ=255)
ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                    >  1:3362788‑M1/1‑68 (MQ=255)
ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                    >  1:1599937‑M1/1‑68 (MQ=255)
 gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGt                                                                   <  1:1760820‑M1/68‑2 (MQ=255)
 gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGt                                                                   <  1:2137920‑M1/68‑2 (MQ=255)
 gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGt                                                                   <  1:3641421‑M1/68‑2 (MQ=255)
  ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                  >  1:933196‑M1/1‑66 (MQ=255)
  ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                  >  1:2444338‑M1/1‑66 (MQ=255)
  ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                  >  1:2931359‑M1/1‑66 (MQ=255)
  ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                  >  1:3564491‑M1/1‑66 (MQ=255)
  ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                  >  1:506981‑M1/1‑66 (MQ=255)
  ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                  >  1:2017766‑M1/1‑66 (MQ=255)
  ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGct                                                                  >  1:3225707‑M1/1‑66 (MQ=255)
    ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtata                                                                >  1:1786416‑M1/1‑64 (MQ=255)
     gCTGCGATTTTCCATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                               >  1:2554660‑M1/1‑63 (MQ=255)
     gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                               >  1:500266‑M1/1‑63 (MQ=255)
     gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                               >  1:138734‑M1/1‑63 (MQ=255)
      cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                              <  1:1613140‑M1/68‑7 (MQ=255)
      cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                              <  1:3395938‑M1/68‑7 (MQ=255)
       tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                             >  1:2345969‑M1/1‑61 (MQ=255)
             tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                       <  1:1405816‑M1/68‑14 (MQ=255)
                 aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                   <  1:1682138‑M1/68‑18 (MQ=255)
                    aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttt                                                >  1:1940464‑M1/1‑48 (MQ=255)
                        cGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                            <  1:1654549‑M1/68‑25 (MQ=255)
                         gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                           <  1:956040‑M1/68‑26 (MQ=255)
                         gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                           <  1:1320655‑M1/68‑26 (MQ=255)
                         gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                           <  1:740164‑M1/68‑26 (MQ=255)
                           gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcntactttagttttgggtgt                                         <  1:3093819‑M1/68‑28 (MQ=255)
                             tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                       <  1:2215847‑M1/68‑30 (MQ=255)
                             tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                       <  1:941600‑M1/68‑30 (MQ=255)
                             tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                       <  1:3171131‑M1/68‑30 (MQ=255)
                             tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                       <  1:553504‑M1/68‑30 (MQ=255)
                              gCAGATCACACAAAACACTCAATTACTTATCATAAATGtataagtcatacttttgctttgggtgtatt                                      >  1:3591117‑M1/1‑38 (MQ=255)
                              gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                      >  1:2680484‑M1/1‑38 (MQ=255)
                              gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                      >  1:3762938‑M1/1‑38 (MQ=255)
                              gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                      >  1:1958089‑M1/1‑38 (MQ=255)
                                aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                    >  1:3337903‑M1/1‑36 (MQ=255)
                                   tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                 >  1:1644574‑M1/1‑33 (MQ=255)
                                   tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                 >  1:1646669‑M1/1‑33 (MQ=255)
                                   tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                 >  1:1653909‑M1/1‑33 (MQ=255)
                                   tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                 >  1:1729542‑M1/1‑33 (MQ=255)
                                   tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                 >  1:3758831‑M1/1‑33 (MQ=255)
                                   tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                 >  1:3251915‑M1/1‑33 (MQ=255)
                                   tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                 >  1:2523531‑M1/1‑33 (MQ=255)
                                   tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                 >  1:1804633‑M1/1‑33 (MQ=255)
                                   tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                 >  1:3693614‑M1/1‑33 (MQ=255)
                                   tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtacttcaat                                 >  1:2814318‑M1/1‑33 (MQ=255)
                                   tCACACAAAACACTCAATTACTTAACATAAATGtataagacatacttctgtcctggctgtatttcact                                 >  1:3578388‑M1/1‑33 (MQ=255)
                                     acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                               >  1:2932412‑M1/1‑31 (MQ=255)
                                     acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                               >  1:3721411‑M1/1‑31 (MQ=255)
                                       acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                             >  1:3491295‑M1/1‑29 (MQ=255)
                                       acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                             >  1:2739543‑M1/1‑29 (MQ=255)
                                       acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                             >  1:322320‑M1/1‑29 (MQ=255)
                                       acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                             <  1:1459957‑M1/68‑40 (MQ=255)
                                       acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                             >  1:1751830‑M1/1‑29 (MQ=255)
                                        caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                            >  1:65954‑M1/1‑28 (MQ=255)
                                            acacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                        >  1:1101676‑M1/1‑24 (MQ=255)
                                            acacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                        >  1:1348961‑M1/1‑24 (MQ=255)
                                            acacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                        >  1:3771944‑M1/1‑24 (MQ=255)
                                            acacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                        >  1:3042069‑M1/1‑24 (MQ=255)
                                             cacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcgatctgttaaaaa                       >  1:3593599‑M1/1‑23 (MQ=255)
                                               cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagt                     >  1:1962418‑M1/1‑21 (MQ=255)
                                               cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagt                     >  1:610374‑M1/1‑21 (MQ=255)
                                                 cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagcct                   >  1:1778232‑M1/1‑19 (MQ=255)
                                                   aTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                 >  1:3854686‑M1/1‑17 (MQ=255)
                                                      aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc              >  1:1751954‑M1/1‑14 (MQ=255)
                                                      aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc              >  1:912632‑M1/1‑14 (MQ=255)
                                                      aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc              >  1:70598‑M1/1‑14 (MQ=255)
                                                       cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct             <  1:1437735‑M1/68‑56 (MQ=255)
                                                        ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta            >  1:1878599‑M1/1‑12 (MQ=255)
                                                        ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta            >  1:58138‑M1/1‑12 (MQ=255)
                                                         tAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac           <  1:721829‑M1/68‑58 (MQ=255)
                                                             aTAAATGtataagtcatanttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta       >  1:3098100‑M1/1‑7 (MQ=255)
                                                             aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta       >  1:2990215‑M1/1‑7 (MQ=255)
                                                             aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta       >  1:205100‑M1/1‑7 (MQ=255)
                                                                  tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaacgtttttcgctacgctagcacg  >  1:3115628‑M1/1‑2 (MQ=255)
                                                                  tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag  >  1:971013‑M1/1‑2 (MQ=255)
                                                                  tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag  >  1:746685‑M1/1‑2 (MQ=255)
                                                                  tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag  >  1:473136‑M1/1‑2 (MQ=255)
                                                                  tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag  >  1:3631358‑M1/1‑2 (MQ=255)
                                                                  tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag  >  1:1915260‑M1/1‑2 (MQ=255)
                                                                  tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag  >  1:1538010‑M1/1‑2 (MQ=255)
                                                                  tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag  >  1:1158194‑M1/1‑2 (MQ=255)
                                                                  tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaac  >  1:81951‑M1/1‑2 (MQ=255)
                                                                   |                                                                  
GGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAG  >  NC_000913/1978435‑1978568

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