Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,962,903 Δ15,600 bp IS1‑mediated [flhE]flhD 15 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1962903 1979165–1978505 15603–16263 126 [0] [12] 14 [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 = 19629020 (0.000)126 (2.030) 93/388 0.0 100% coding (71/393 nt) flhE flagellar protein
?NC_000913 = 3584195 NA (NA)noncoding (768/768 nt) IS1 repeat region

ACGCCATTCATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGC‑‑TTTTGCCGTTACGCACCA‑‑CCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTC‑‑‑AAAAACACCATCATACACTAAATCAGT‑‑AAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACG‑‑CTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAAT  >  NC_000913/1979025‑1979408
                                                                                                                                                    |                                                                                                                                                                                                                                                                  
aCGCCATTCATGGCCAGTATCAAGATTTCTGGTGCGTACCGGGGTTGAGAAGCGGTGAAGTGAACTGCAGTGCCATGTTTACGGCAGTGAGAGCAGAGATAGCGCGTGATGTCCGGCGGTGC‑‑TTTTGCCGTTACGCACCA‑‑CCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTC‑‑‑AAAAACACCATCATACACTAAATCAGT‑‑AAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCa                                                                                                     <  1:378411/285‑1 (MQ=14)
  gCCATTCATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGC‑‑TTTTGCCGTTACGCACCA‑‑CCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTC‑‑‑AAAAACACCATCATACACTAAATCAGT‑‑AAGTTGGCAGCATCACCCATAAATGta                                                                                                                                     <  1:1214096/254‑1 (MQ=33)
                    cAATGATTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGTGC‑‑TTTTGCCGTTACGCACCA‑‑CCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTC‑‑‑AAAAACACCATCATACACTAAATCAGT‑‑AAGTTGGCAGCATCACCCATAAATGTATAAGt                                                                                                                                >  1:1214110/1‑240 (MQ=255)
                                     gTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGC‑‑TTTTGCCGTTACGCACA‑‑CCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTC‑‑‑AAAAACACCATCATACACTAAATCAGT‑‑AAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaaaagttt                                                                                   <  1:1214106/268‑7 (MQ=255)
                                                                                                  cAGAGATAGCGCTGATGTCCGGCGGTGC‑‑TTTTGCCGTTACGCACCA‑‑CCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAAGAAACAGAAGCACCTCAGAGCACCTCA‑‑AAAAACACCATCATACACTAAATCAGT‑‑AAGTTGGCAGCATCACCCATAAATg                                                                                                                                       <  1:378439/161‑1 (MQ=11)
                                                                                                             cTGATGTCCGGCGGTGC‑‑TTTTGCCGTTACGCACCA‑‑CCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTC‑‑‑AAAAACACCATCATACACTAAATCAGT‑‑AAGTTGGCAGCATCACCCATAAATGTATAAGTCATACttttgtt                                                                                                                   <  1:378487/166‑1 (MQ=255)
                                                                                                               ggATGTCCGGCGGGTGCTTTTTTGCCCGTTACGCACCACCCCCCGTCAGTAGCTGAAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCA‑‑AAAAACACCATCATACACTAAAATCAGT‑‑AAGTTGGCAGCATCAcccc                                                                                                                                             <  1:378791/148‑2 (MQ=11)
                                                                                                                gATGTCCGGCGTGCTTTTTGCCGTTACGCACCA‑‑CCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTC‑‑‑AAAAACACCATCATACACTAAATCAGT‑‑AAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTGTTTTGGGTGTAtt                                                                                                        >  1:1214200/1‑174 (MQ=255)
                                                                                                                            gTGC‑‑TTTTGCCGTTACGCACCA‑‑CCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTC‑‑‑AAAAACACCATCATACACTAAATCAGT‑‑AAGTTGGCAGCATCACCCATAAATGTATAAGTCAATACTTTGTTTTGGgtgt                                                                                                           >  1:1214449/1‑161 (MQ=255)
                                                                                                                               c‑‑TTTTGCCGTTACGCACCA‑‑CCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTC‑‑‑AAAAACACCATCATACACTAAATCAGT‑‑AAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTCGCTACGTCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaat  >  1:1214453/1‑260 (MQ=255)
                                                                                                                                          gTTACGCACCACCCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAGCCACTGGAGCACCTC‑‑‑AAAAACACCATCATACACTAAATCAGT‑‑AAGTTGGCAGCATCACCCAt                                                                                                                                            <  1:378580/118‑1 (MQ=255)
                                                                                                                                                  cca‑‑ccCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCCTCAAAAAAAACACCATCATACACTAAATCAGTTAAAGTTGGCAGGCATCAcccga                                                                                                                                            >  1:378808/1‑115 (MQ=11)
                                                                                                                                                    a‑‑ccCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTC‑‑‑AAAAACACCATCATACACTAAATCAGT‑‑AAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACG‑‑CTAGCAAGCAAAAATGa                                                         <  1:1214447/188‑1 (MQ=255)
                                                                                                                                                    a‑‑ccCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTC‑‑‑AAAAACACCATCATACACTAAATCAGT‑‑AAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTGTTT‑‑GGTGTATTTCAATCTGTTAAAAGTTTTTCGCTACGTACTAGCAAGCaaaatg                                                           >  1:1214472/1‑182 (MQ=18)
                                                                                                                                                    |                                                                                                                                                                                                                                                                  
ACGCCATTCATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGC‑‑TTTTGCCGTTACGCACCA‑‑CCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTC‑‑‑AAAAACACCATCATACACTAAATCAGT‑‑AAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACG‑‑CTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAAT  >  NC_000913/1979025‑1979408

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