Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,962,566 Δ16,764 bp flhEinsA‑5 17 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1962566 1979329 16764 21 [0] [0] 22 flhE–insA‑5 flhE,flhA,flhB,cheZ,cheY,cheB,cheR,tap,tar,cheW,cheA,motB,motA,flhC,flhD,insB‑5,insA‑5

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 19625650 (0.000)21 (0.780) 16/94 0.3 100% intergenic (+105/+15) yecT/flhE protein YecT/flagellar protein
?NC_000913 1979330 = 0 (0.000)intergenic (‑115/‑423) insA‑5/uspC IS1 protein InsA/universal stress protein C

GACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTAATTTTTTTGTCCTTCAGCGATAATTCACAATCACTT  >  NC_000913/1962518‑1962602
                                               |                                     
gACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTc                                      <  1:1574824‑M1/49‑2 (MQ=255)
  cATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgc                                    >  1:1057646‑M1/1‑46 (MQ=255)
  cATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgc                                    >  1:1333012‑M1/1‑46 (MQ=255)
   aTCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgc                                    >  1:1378281‑M1/1‑45 (MQ=255)
    tCGAGCCTTTTTCACTGAGTTATTAAACATACTCGTGAGCGCGTcgcta                                  <  1:1607127‑M1/49‑6 (MQ=255)
      gAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctac                                 <  1:2059255‑M1/48‑7 (MQ=255)
         ccTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgcta                             <  1:712393‑M1/49‑11 (MQ=255)
         ccTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgcta                             <  1:1932435‑M1/49‑11 (MQ=255)
         ccTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgc                               >  1:981333‑M1/1‑39 (MQ=255)
           tttttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctag                            >  1:1080119‑M1/1‑37 (MQ=255)
             tttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagc                           <  1:713784‑M1/47‑13 (MQ=255)
             tttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagc                           <  1:1035585‑M1/47‑13 (MQ=255)
                cACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagc                       <  1:129959‑M1/48‑17 (MQ=255)
                cACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagc                       <  1:1875743‑M1/48‑17 (MQ=255)
                    gAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaa                  >  1:107603‑M1/1‑28 (MQ=255)
                     aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaat                 >  1:1776521‑M1/1‑27 (MQ=255)
                         attaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaa              <  1:1226292‑M1/48‑26 (MQ=255)
                           taAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaa             >  1:437930‑M1/1‑21 (MQ=255)
                               cATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaa       >  1:285376‑M1/1‑17 (MQ=255)
                               cATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaa       >  1:1512822‑M1/1‑17 (MQ=255)
                                     cgcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatc  >  1:1153398‑M1/1‑11 (MQ=255)
                                               |                                     
GACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTAATTTTTTTGTCCTTCAGCGATAATTCACAATCACTT  >  NC_000913/1962518‑1962602

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