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–1979269 1979278 10–776 41 [0] [0] 41 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)39 (0.850) 35/274 0.3 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

AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913_3_pae_tpiA/1978357‑1978644
                                                                                                                                                 |                                                                                                                                              
aaaTGGTGCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                             <  2:68266‑M1/149‑4 (MQ=255)
aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATAGACGCAACTGTACTCGTCACTAACAGCACACACAACGGAGAGGGGCTGCGATTTTCAATAATGCGTGATGCACATCAAACAAAACACCCAATTACTTAACATCAATGtct                                                                                                                                             >  1:375102‑M1/1‑146 (MQ=255)
 aaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtata                                                                                                                                            >  2:285526‑M1/1‑145 (MQ=255)
     ggTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtc                                                                                                                                        <  1:122584‑M1/149‑9 (MQ=255)
      gTCTTTACACTTATCCAAGAGTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGAACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                       <  2:158951‑M1/149‑10 (MQ=255)
              acTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACCAGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTGCTTAACATAAATGtataagtcatacttttg                                                                                                                               <  1:374322‑M1/149‑18 (MQ=255)
                  aTCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                                                           <  1:268485‑M1/149‑22 (MQ=255)
                            tttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                                                                 <  1:749453‑M1/149‑32 (MQ=255)
                                   aaccgacGCAACTGTACTCGGCACCACACGGAAATACAACGGAGGGCGGCGGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctg                                                                                                          <  1:175664‑M1/146‑39 (MQ=255)
                                        aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                                                                     >  1:394158‑M1/1‑106 (MQ=255)
                                               tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                                                              >  1:313591‑M1/1‑99 (MQ=255)
                                                    tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                         >  1:96441‑M1/1‑94 (MQ=255)
                                                              acGCGCATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataattcatacttttgttttgggtgtatttctatctgttaaacagttattagctagtctagcaa                                                                               >  1:272881‑M1/1‑84 (MQ=255)
                                                              acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                               >  1:387892‑M1/1‑84 (MQ=255)
                                                                      aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                       >  1:744070‑M1/1‑76 (MQ=255)
                                                                        aaCGGAGGGGCGCTGCTATTTTAAATAAGGAGTGAGGCAGATCACAAAAAACACTCAAGTACTTAACAAAAATGtataagtcattcttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatga                                                                     <  2:444082‑M1/149‑76 (MQ=255)
                                                                          cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                         <  1:102184‑M1/111‑40 (MQ=255)
                                                                          cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                         >  2:102184‑M1/1‑72 (MQ=255)
                                                                              ggggggCTGCGATTTTCAATAATCCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaacaatgaaacagg                                                               <  1:521499‑M1/149‑82 (MQ=255)
                                                                               gggggCGGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagga                                                              <  1:479760‑M1/149‑83 (MQ=255)
                                                                                   gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa                                                          <  2:96441‑M1/149‑87 (MQ=255)
                                                                                   gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa                                                          <  2:387892‑M1/149‑87 (MQ=255)
                                                                                   gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgataaaaagtttttcgctacgctagcaaccaaaaatgaaacaggaataa                                                          >  2:728451‑M1/1‑63 (MQ=255)
                                                                                    cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                         <  1:568926‑M1/149‑88 (MQ=255)
                                                                                        gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaacatgaaaaagaaataatagaa                                                     >  2:344863‑M1/1‑58 (MQ=255)
                                                                                         aTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa                                                    >  1:220642‑M1/1‑57 (MQ=255)
                                                                                                  aatGCGTGATTCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgtt                                           <  1:640004‑M1/149‑102 (MQ=255)
                                                                                                      cGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaagctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                       >  1:583726‑M1/1‑44 (MQ=255)
                                                                                                           tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                  >  1:585281‑M1/1‑39 (MQ=255)
                                                                                                              aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaa                               <  1:715579‑M1/149‑114 (MQ=255)
                                                                                                                   acaaaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                          <  2:356543‑M1/149‑119 (MQ=255)
                                                                                                                     gcaaAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcccacagtcaaaataact                        >  2:462709‑M1/2‑29 (MQ=255)
                                                                                                                       aaaaCAATCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactca                      <  2:641394‑M1/149‑123 (MQ=255)
                                                                                                                            acTCAATTACTTAACATAAATGtataagtgatacttttgttttggctgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgtgatgttgcgcacagtacaaatacctccccgta                 >  2:773472‑M1/1‑22 (MQ=255)
                                                                                                                            acTCAATTACTTAACATAAATGtataagtcatacttttgttttggatgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgtaatgttgcgcacagtcaaaataattcaccgta                 >  2:309306‑M1/1‑22 (MQ=255)
                                                                                                                                      ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagcatttagcttctgtagcgtgctaaaatgaaacaagagtaattgaaatgggttgttgagcacagtcaaaataactcaccgtaaataatcatc       >  2:775800‑M1/1‑12 (MQ=255)
                                                                                                                                       tAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatct      <  2:744070‑M1/149‑139 (MQ=255)
                                                                                                                                        tatataaATGtataagtcatacttttgttttgggtttatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaacgaaacaagtatgttctagatgggattttgcgcacagtcaaaataactcaccgtaaataatcatctg     <  2:393441‑M1/146‑140 (MQ=255)
                                                                                                                                        aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgatacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaactaatcatctg     >  1:22857‑M1/1‑10 (MQ=255)
                                                                                                                                        aaCATAAATGtataagtcatactatcgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg     <  1:644581‑M1/149‑140 (MQ=255)
                                                                                                                                           aTAAATGtataagtcatactgttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta  <  2:411408‑M1/149‑143 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913_3_pae_tpiA/1978357‑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: