Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 2,795,396 C→A Q382K (CAG→AAG)  gabP → gamma‑aminobutyrate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA2,795,3960CA100.0% 136.8 / NA 40Q382K (CAG→AAG) gabPgamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (19/21);  total (19/21)

GTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCA  >  NC_000913_3_pae_tpiA/2795266‑2795532
                                                                                                                                  |                                                                                                                                        
gtTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCaaaa                                                                                                                        <  2:87350/149‑1 (MQ=255)
gtTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCCGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCaaaa                                                                                                                        <  2:72777/149‑1 (MQ=255)
 ttACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAAt                                                                                                                       >  2:96398/1‑149 (MQ=255)
 ttACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAAt                                                                                                                       >  2:232640/1‑149 (MQ=255)
 ttACTCTCCACCGGCGCGGCATTCTTAACGGTGCCGGTCCACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAAt                                                                                                                       >  2:431598/1‑149 (MQ=255)
 ttACTCCCCACCGGCGCGGCATTCTAAACGATGATGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTCGTTTATTTAGTCATCGCCGTTTCAAAGTAGCGGATGCGCAAAAt                                                                                                                       >  2:431616/1‑149 (MQ=255)
      ctcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGc                                                                                                                  <  2:245247/149‑1 (MQ=255)
        ccACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGa                                                                                                                <  2:93759/149‑1 (MQ=255)
              gcgcGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGGTGCGCAAAATTCTGCGAGCAgaa                                                                                                          >  1:198347/1‑149 (MQ=255)
                     aTTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCg                                                                                                   <  1:232640/149‑1 (MQ=255)
                          tAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAAtt                                                                                              <  1:85407/149‑1 (MQ=255)
                             aagtgatggtgAACTAGTACGCGCCGGCGAAGGTATCTAAATTTCGGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTAAAAAAAAAAAAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGc                                                                                           <  2:44532/147‑1 (MQ=255)
                              ggtggtggtGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCt                                                                                          <  2:151209/149‑1 (MQ=255)
                                      tgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGACGGAAGCGAAATTCGCTTTCGAAtc                                                                                  >  2:102093/1‑148 (MQ=255)
                                          cTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGc                                                                              >  2:505757/1‑149 (MQ=255)
                                           tATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCt                                                                             <  1:96398/149‑1 (MQ=255)
                                             ttACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGCGCAGACGGAAGCAAAATTCGCTTGCGAATGTGGCttt                                                                           >  2:41069/1‑149 (MQ=255)
                                                 gcgcCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATcc                                                                       >  1:375938/1‑149 (MQ=255)
                                                   gcCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGt                                                                     >  1:394085/1‑149 (MQ=255)
                                                    ccTGCGAAGGTATTGAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTg                                                                    <  1:55991/149‑1 (MQ=255)
                                                          aaGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCAc                                                              <  1:505757/149‑1 (MQ=255)
                                                            ggTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCt                                                            <  1:344625/149‑1 (MQ=255)
                                                             gTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGGCGCGGCGTGCGCTAGCGAATGTGGCTTTATACGTGGCTCACCTg                                                           >  2:298890/1‑149 (MQ=255)
                                                                      tttCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCAt                                                  >  1:107112/1‑149 (MQ=255)
                                                                           gATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCt                                             >  2:476984/1‑149 (MQ=255)
                                                                            aTCGTCAGCTCCGGCGCTTTCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCtt                                            <  2:60362/149‑1 (MQ=255)
                                                                             tCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTCTTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGTGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTCTCCGTGGCCCACCTGGCTAGTCATCGGCatt                                           >  2:139708/1‑149 (MQ=255)
                                                                              cGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGCGGCTAAAGCCCTTTCTCACCTGGCTAGTCATCGGCTtta                                          <  1:476984/149‑1 (MQ=255)
                                                                                  aGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTACGGAGGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTAc                                      >  1:503986/1‑149 (MQ=255)
                                                                                      ccGGCGGTATGCCCCTGGTGGTTTATTTAGTCATCGCCGGTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCttt                                  <  2:256889/149‑1 (MQ=255)
                                                                                       cGGCGCTATCGCCCTGATGGTTTATTTAGTCATCGCCGTTTCAAAAATGCGGATGCGCAAAATTCTGAGATCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGCCTTACATGGCTAGCAATCGGCTTTATTGCCtact                                 >  2:189123/1‑146 (MQ=255)
                                                                                                cGCCCTGCTGGGTTATTTAGTCAGCGCCGTTTCAAAGTTGCGGATGCGCAAAATTGTGGGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAg                        <  2:583051/149‑1 (MQ=255)
                                                                                                    ctgctgGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGAt                    >  1:205861/1‑149 (MQ=255)
                                                                                                          gTTTATTTAGGCAACGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTAtt              <  2:210821/149‑1 (MQ=255)
                                                                                                              agttaGTCATCGCCGTTTAAAAGGTGCGGAGGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCCCACCGGGCGAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtccgt          <  1:603238/147‑1 (MQ=255)
                                                                                                                ttaTTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtccgtcc        <  1:249619/149‑1 (MQ=255)
                                                                                                                 taGTCATCGCCGTTTCAAAGTTGCGGATGCGAAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtccgtccg       <  2:448977/149‑1 (MQ=255)
                                                                                                                      aTCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCTGTCCGgcgct  >  1:425218/1‑148 (MQ=255)
                                                                                                                       tCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTa                <  1:563741/134‑1 (MQ=255)
                                                                                                                       tCGCCGTTTCAAAGTTGCGGAGGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTa                >  2:563741/1‑134 (MQ=255)
                                                                                                                                  |                                                                                                                                        
GTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCA  >  NC_000913_3_pae_tpiA/2795266‑2795532

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: ‑