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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1978503 1979278 776 17 [0] [0] 19 insB1–insA insB1,insA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 19785020 (0.000)16 (0.400) 14/260 NT 100% intergenic (‑305/+16) flhD/insB1 flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
?NC_000913 1979279 = 0 (0.000)intergenic (‑64/‑474) insA/uspC IS1 repressor TnpA/universal stress protein

ATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCA  >  NC_000913/1978359‑1978633
                                                                                                                                               |                                                                                                                                   
aTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                <  1:230419‑M1/149‑6 (MQ=255)
 tGGGTCTTTACTCTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATCCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                                                                                               <  2:254161‑M1/149‑7 (MQ=255)
                  cTAAGATTTTGCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTTAAGTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                              <  1:294631‑M1/149‑24 (MQ=255)
                                 aaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctg                                                                                               <  2:141761‑M1/149‑39 (MQ=255)
                                  aTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctat                                                                                              >  1:162410‑M1/1‑110 (MQ=255)
                                             tGTACTCGTCACTTCACGTACATACAACTGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcattcttttgttttgggtgtatttcaatctgttaaaaagtttt                                                                                   <  2:373489‑M1/149‑51 (MQ=255)
                                              gTACTCGTCACTCCACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                  <  1:656765‑M1/149‑52 (MQ=255)
                                                       aCTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgct                                                                         >  1:497178‑M1/1‑89 (MQ=255)
                                                         tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag                                                                       >  1:433619‑M1/1‑87 (MQ=255)
                                                                       aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                                                         >  1:620534‑M1/1‑73 (MQ=255)
                                                                          gatggtggctgcgaTTTTCAATTATGCGTGATGCAGATCACACAAAACACTCAATTACGTAACATAAATGtataagtcatacttctgttttgggtgtatttcaacctgtcaaaaagtttttcgctacgctagcaagcaaaaatgaaaca                                                      <  2:475809‑M1/143‑80 (MQ=255)
                                                                          gAGGGGGTCTGCGATTTTCCAGAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaca                                                      <  2:127306‑M1/149‑80 (MQ=255)
                                                                          gAGGGGGGCGGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaca                                                      <  2:285709‑M1/149‑80 (MQ=255)
                                                                                  cagcgATTTTAAATAATGCGTGATACAGATCAAACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttcggtgtatttcaatctgttaaaaagtttttcactaaaataaaaaacaaaaatgaaacaggaataat                                              <  1:6446‑M1/147‑88 (MQ=255)
                                                                                                 atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagacatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatggaacaggaataatcaaaatgggatgttg                               >  2:48002‑M1/1‑47 (MQ=255)
                                                                                                            aGATCACAAAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaa                    <  1:310481‑M1/149‑114 (MQ=255)
                                                                                                                              aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata  >  1:141162‑M1/1‑18 (MQ=255)
                                                                                                                                               |                                                                                                                                   
ATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCA  >  NC_000913/1978359‑1978633

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: