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

AGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913_3_pae_tpiA/1978355‑1978644
                                                                                                                                                   |                                                                                                                                              
aGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGt                                                                                                                                               >  1:490358‑M1/1‑148 (MQ=255)
  aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACACAACACTCAATTACTTAACATAAATGtat                                                                                                                                             >  1:303769‑M1/1‑146 (MQ=255)
    aTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                           >  1:156586‑M1/1‑144 (MQ=255)
               cacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                                <  1:303569‑M1/149‑17 (MQ=255)
                acTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttt                                                                                                                               >  1:224263‑M1/1‑132 (MQ=255)
                                 ccTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaa                                                                                                              >  1:210856‑M1/1‑115 (MQ=255)
                                  cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttttggtgtatttcaat                                                                                                             >  2:453453‑M1/1‑114 (MQ=255)
                                  cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                             >  2:245339‑M1/1‑114 (MQ=255)
                                       gcgccggaacggtactCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAAGAATGCGTGATGCAGGTCACACAAAACACTCAAGTACTTAACATAAATGtataaggcatacttttgttttgggtgtatttcaatctgtt                                                                                                        <  1:94386‑M1/142‑41 (MQ=255)
                                          aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                                                                     >  1:416880‑M1/1‑106 (MQ=255)
                                                 tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                           >  2:308558‑M1/1‑99 (MQ=255)
                                                 tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                           <  1:308558‑M1/104‑6 (MQ=255)
                                                      tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtctatttctatctgttaaaaatttttttgct                                                                                         >  2:10820‑M1/1‑94 (MQ=255)
                                                            cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                                                                   <  2:490358‑M1/149‑62 (MQ=255)
                                                                      aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaa                                                                         >  2:290191‑M1/1‑78 (MQ=255)
                                                                      aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaa                                                                         <  1:263719‑M1/149‑72 (MQ=255)
                                                                                        gCGAGTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcagacttttgttttgggtctatttcaatttgttaaaacgtttttagctacactagcaacgaaaaataaaagaggaataatcg                                                       >  2:194322‑M1/1‑60 (MQ=255)
                                                                                                         gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaacaagtttttcgctactctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgca                                      <  1:436768‑M1/149‑107 (MQ=255)
                                                                                                              gaagaTAACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatatttcgcacagtc                                 >  2:335308‑M1/3‑38 (MQ=255)
                                                                                                              gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggggtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaattaaacaggaataagcggaatgggatgttgcgcacggtc                                 >  1:353772‑M1/1‑38 (MQ=255)
                                                                                                                 gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctaaaaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaa                              <  2:306468‑M1/149‑115 (MQ=255)
                                                                                                                   tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                            >  1:335205‑M1/1‑33 (MQ=255)
                                                                                                                     acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                          <  2:224263‑M1/149‑119 (MQ=255)
                                                                                                                         aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactca                      >  1:122851‑M1/1‑27 (MQ=255)
                                                                                                                                tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaa               >  1:21110‑M1/1‑20 (MQ=255)
                                                                                                                                   aTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataa            >  1:399734‑M1/1‑17 (MQ=255)
                                                                                                                                        ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatc       >  1:384251‑M1/1‑12 (MQ=255)
                                                                                                                                             aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctactctagcaagcaaaaatgaaacaggaataatgaaaataggattttccgcacagtcaaaataacttaccataaataactatctgcta  >  2:34288‑M1/1‑7 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
AGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913_3_pae_tpiA/1978355‑1978644

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: