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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913_3_pae_tpiA 1978503–1979266 1979278 13–776 28 [0] [0] 28 insB1–insA insB1,insA

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

GCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTA  >  NC_000913_3_pae_tpiA/1979131‑1979420
                                                                                                                                                    |                                                                                                                                             
agaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgt                                                                                                                                               >  1:490358‑M2/149‑149 (MQ=255)
  aaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacacaacactcaattacttaacataaatgtat                                                                                                                                             >  1:303769‑M2/147‑149 (MQ=255)
    atgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                                                                           >  1:156586‑M2/145‑149 (MQ=255)
               cacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtttt                                                                                                                                <  1:303569‑M2/16‑1 (MQ=255)
                acttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATActtttt                                                                                                                               >  1:224263‑M2/133‑148 (MQ=255)
                                 cctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCaa                                                                                                              >  1:210856‑M2/116‑149 (MQ=255)
                                  ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTTGGTGTATTTCAAt                                                                                                             >  2:453453‑M2/115‑149 (MQ=255)
                                  ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                             >  2:245339‑M2/115‑149 (MQ=255)
                                       gcgccggaacggtactcgtcactacacgcacatacaacggaggggggctgcgattttcaagaatgcgtgatgcaggtcacacaaaacactcaagtacttaacataaatgTATAAGGCATACTTTTGTTTTGGGTGTATTTCAATCTGtt                                                                                                        <  1:94386‑M2/40‑1 (MQ=255)
                                          acgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaa                                                                                                     >  1:416880‑M2/107‑149 (MQ=255)
                                                 tgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                                                                           >  2:308558‑M2/100‑104 (MQ=255)
                                                 tgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                                                                           <  1:308558‑M2/5‑1 (MQ=255)
                                                      tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTCTATTTCTATCTGTTAAAAATTTTTTtgct                                                                                         >  2:10820‑M2/95‑149 (MQ=255)
                                                            ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                                                                   <  2:490358‑M2/61‑1 (MQ=255)
                                                                      atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaa                                                                         >  2:290191‑M2/79‑149 (MQ=255)
                                                                      atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaa                                                                         <  1:263719‑M2/71‑1 (MQ=255)
                                                                                        gcgagtttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCAGACTTTTGTTTTGGGTCTATTTCAATTTGTTAAAACGTTTTTAGCTACACTAGCAACGAAAAATAAAAGAGGAATAATCg                                                       >  2:194322‑M2/61‑149 (MQ=255)
                                                                                                         gtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAACAAGTTTTTCGCTACTCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGca                                      <  1:436768‑M2/106‑1 (MQ=255)
                                                                                                              gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGGGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATTAAACAGGAATAAGCGGAATGGGATGTTGCGCACGGTc                                 >  1:353772‑M2/39‑149 (MQ=255)
                                                                                                              gaagataacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATATTTCGCACAGTc                                 >  2:335308‑M2/39‑149 (MQ=255)
                                                                                                                 gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAAAAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaa                              <  2:306468‑M2/114‑1 (MQ=255)
                                                                                                                   tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAAt                            >  1:335205‑M2/34‑149 (MQ=255)
                                                                                                                     acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                          <  2:224263‑M2/118‑1 (MQ=255)
                                                                                                                         aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCa                      >  1:122851‑M2/28‑149 (MQ=255)
                                                                                                                                tcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaaa               >  1:21110‑M2/21‑149 (MQ=255)
                                                                                                                                   attacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaataa            >  1:399734‑M2/18‑149 (MQ=255)
                                                                                                                                        ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcatc       >  1:384251‑M2/13‑149 (MQ=255)
                                                                                                                                             ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACTCTAGCAAGCAAAAATGAAACAGGAATAATGAAAATAGGATTTTCCGCACAGTCAAAATAACTTACCATAAATAACTATCTGCta  >  2:34288‑M2/8‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTA  >  NC_000913_3_pae_tpiA/1979131‑1979420

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: