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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1978503–1979257 1979278 22–776 11 [0] [0] 11 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)11 (1.140) 10/248 NT 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

ACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATT  >  NC_000913/1978371‑1978631
                                                                                                                                   |                                                                                                                                 
acTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcata                                                                                                                        <  2:43680‑M1/143‑12 (MQ=255)
         aGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTCCCCGCACATACAACGGAGTGGGGCTGAGATTTTCAATAATACGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtct                                                                                                               >  1:75736‑M1/1‑123 (MQ=255)
                      aTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttca                                                                                                  >  2:166924‑M1/1‑110 (MQ=255)
                          aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactattgtttatggggtatttcaatct                                                                                              >  1:74199‑M1/1‑106 (MQ=255)
                                 tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAAGAATGCCTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                                                                                       <  1:19776‑M1/143‑45 (MQ=255)
                                                                     gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                   <  2:118289‑M1/143‑81 (MQ=255)
                                                                                      tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAACTGtataagtcatactttggttttggatgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga                                  >  2:33594‑M1/1‑46 (MQ=255)
                                                                                                 gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                       <  1:142386‑M1/143‑109 (MQ=255)
                                                                                                     atacaagACACTCAATTACTAAACATAAATGtataagtcatacttttgttttggttgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtca                   <  1:122863‑M1/141‑113 (MQ=255)
                                                                                                     acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtca                   <  1:23479‑M1/143‑113 (MQ=255)
                                                                                                                      aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttacgctacggaagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaa  <  2:110078‑M1/143‑130 (MQ=255)
                                                                                                                                   |                                                                                                                                 
ACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATT  >  NC_000913/1978371‑1978631

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