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

GACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1962518‑1962565
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATC  >  NC_000913/1979330‑1979366
                                                                                     
GACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTC                                      <  1:1574824/49‑1
  CATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTCGC                                    >  1:1057646/1‑49
  CATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTCGC                                    >  1:1333012/1‑49
   ATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTCGC                                    >  1:1378281/1‑48
    TCGAGCCTTTTTCACTGAGTTATTAAACATACTCGTGAGCGCGTCGCTA                                  <  1:1607127/49‑1
      GAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTCGCTAC                                 <  1:2059255/48‑1
         CCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTCGCTACGCTA                             <  1:712393/49‑1
         CCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTCGCTACGCTA                             <  1:1932435/49‑1
         CCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTCGCTACGC                               >  1:981333/1‑47
           TTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTCGCTACGCTAG                            >  1:1080119/1‑48
             TTTCACTGAGTTATTAAACATACTCGCGAGCGCGTCGCTACGCTAGC                           <  1:713784/47‑1
             TTTCACTGAGTTATTAAACATACTCGCGAGCGCGTCGCTACGCTAGC                           <  1:1035585/47‑1
                CACTGAGTTATTAAACATACTCGCGAGCGCGTCGCTACGCTAGCAAGC                       <  1:129959/48‑1
                CACTGAGTTATTAAACATACTCGCGAGCGCGTCGCTACGCTAGCAAGC                       <  1:1875743/48‑1
                    GAGTTATTAAACATACTCGCGAGCGCGTCGCTACGCTAGCAAGCAAAAA                  >  1:107603/1‑49
                     AGTTATTAAACATACTCGCGAGCGCGTCGCTACGCTAGCAAGCAAAAAT                 >  1:1776521/1‑49
                         ATTAAACATACTCGCGAGCGCGTCGCTACGCTAGCAAGCAAAAATGAA              <  1:1226292/48‑1
                           TAAACATACTCGCGAGCGCGTCGCTACGCTAGCAAGCAAAAATGAAA             >  1:437930/1‑47
                               CATACTCGCGAGCGCGTCGCTACGCTAGCAAGCAAAAATGAAACAGGAA       >  1:285376/1‑49
                               CATACTCGCGAGCGCGTCGCTACGCTAGCAAGCAAAAATGAAACAGGAA       >  1:1512822/1‑49
                                     CGCGAGCGCGTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATC  >  1:1153398/1‑48
                                                                                     
GACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1962518‑1962565
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATC  >  NC_000913/1979330‑1979366

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: 
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.