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,3960CA98.0% 160.8 / NA 49Q382K (CAG→AAG) gabPgamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (28/20);  minor base T (1/0);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTA  >  NC_000913_3_pae_tpiA/2795264‑2795539
                                                                                                                                    |                                                                                                                                               
gtgtTACTCTCCACCGGCGCGGCATTATTAACGGTGGTGGTGAACTATTACGCGCCTGTGAAGGTATGTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTTTTTGGTCATCGAGGTTTCATAGTTTCGGGTGGGCaa                                                                                                                                 >  2:42683/1‑149 (MQ=255)
  gtTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGGATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAATGTTACGGATGCGCaaac                                                                                                                               >  1:384635/1‑148 (MQ=255)
       tctcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAACGTTGCGGATGCGCAAAATTCTg                                                                                                                          >  1:3956/1‑149 (MQ=255)
          ccACCGGCGCGGCATTCTTAACGGTGGTGGTGAACCCTTACGCGGCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGa                                                                                                                       <  1:459366/149‑1 (MQ=255)
                       aTTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGTAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCg                                                                                                          >  1:621703/1‑149 (MQ=255)
                           ttAACGGTGGTGGTGAACTATTACGCGCGTGCGAAGGTATTTAAATTTCTGATAGAAAGATCAGGCGGGATGGGACTGATGGTTTATTTAGTCATCGCCGTTTCAAAGTTGAGGATGCGCAAAATCTTACGAGCAGAAGGACGCGAAAt                                                                                                      >  2:45452/1‑149 (MQ=255)
                                ggtggtggtGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAACTTCTGCGAGCAGAAGGAAGCGAAATTCGCt                                                                                                 <  2:425782/149‑1 (MQ=255)
                                     tggcgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGAACGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCTTCACAGTTTAAAAGTTCTGTATTCGCAAAATTCTGCGTACATCAGGAAGCACAATTCGCTTGcca                                                                                            >  2:653955/1‑147 (MQ=255)
                                       gtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGacg                                                                                          >  2:87385/1‑147 (MQ=255)
                                                      ccTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTg                                                                           <  1:854433/149‑1 (MQ=255)
                                                       cTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTgg                                                                          <  1:873560/149‑1 (MQ=255)
                                                            aaGGAATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCAc                                                                     >  1:999816/1‑149 (MQ=255)
                                                               gTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCCAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTg                                                                  >  2:104195/1‑149 (MQ=255)
                                                                       aTTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGGAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTTGCTCACCTGGCTAGTCa                                                          >  1:99169/1‑149 (MQ=255)
                                                                         ttCTTAGCGACAGCTCCGGCGCTTTCGCGATGCTGGTTTATTTAGTCATCGCCGTTTCAAATTTGCGGATGTGCATATTTCTGCGAGCAGAAGGAAGTGCAATTCGCTTGCGAATGTGGCTTTATCCTTGGTTCACCTTGCTAGTCATc                                                        <  2:858939/149‑1 (MQ=255)
                                                                            tGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGc                                                     >  2:440760/1‑149 (MQ=255)
                                                                             ggccGACAGCTCCGTCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCt                                                    <  2:1050583/146‑1 (MQ=255)
                                                                             gATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCt                                                    >  1:645759/1‑149 (MQ=255)
                                                                             gATCGACAGCCCCGGCGCTATCGCCCTGCTGGTTTATTTCGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTTGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCt                                                    <  2:472121/149‑1 (MQ=255)
                                                                              aTCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGAGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCtt                                                   >  2:474070/1‑149 (MQ=255)
                                                                                    aGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTAc                                             <  2:757796/149‑1 (MQ=255)
                                                                                      cTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGc                                                                                     <  1:975945/107‑1 (MQ=255)
                                                                                      cTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGc                                                                                     >  2:975945/1‑107 (MQ=255)
                                                                                      cTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCt                                           >  2:507755/1‑149 (MQ=255)
                                                                                            cgcTATCCCACTGCTAGTTTATTAAGTCATCTCCGTTTCAAAATTCCGGATTCCCAAAATTCTCTAAGCAGAAGGAACCAAAAATCGCTTGCGAATGTGGCTTTATCCTTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTtgttt                                     <  2:976886/149‑3 (MQ=255)
                                                                                            ccatATCGCCCGGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTCGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTgtgt                                     <  1:638657/146‑1 (MQ=255)
                                                                                               tATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCg                                                                             >  2:53195/1‑106 (MQ=255)
                                                                                               tATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCg                                                                             <  1:53195/106‑1 (MQ=255)
                                                                                                 tCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTACCCTGGCTAGTCATCGGCTTTATTACATTTGTGTTGGTa                                >  2:961004/1‑149 (MQ=255)
                                                                                                  cGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCGCCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAg                               >  1:982778/1‑149 (MQ=255)
                                                                                                   gCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCAc                                                                     <  2:981519/110‑1 (MQ=255)
                                                                                                   gCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCAc                                                                     >  1:981519/1‑110 (MQ=255)
                                                                                                      ctgctgGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGAt                           <  1:7646/149‑1 (MQ=255)
                                                                                                          tgGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGGGATGCTa                       >  2:744144/1‑149 (MQ=255)
                                                                                                            gTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTAtt                     >  1:42204/1‑149 (MQ=255)
                                                                                                             tttatttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtc                    <  1:636610/149‑1 (MQ=255)
                                                                                                               tatttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTttat                                                >  2:147278/1‑119 (MQ=255)
                                                                                                               tatttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTttat                                                <  1:147278/119‑1 (MQ=255)
                                                                                                                 tttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtccgtc                <  1:582086/149‑1 (MQ=255)
                                                                                                                 ttctgTCACCGCCGTTTCAAAGGTGCGGATGCGCACAATTCTCCGATCAGAAGGAAGCGAAATTCGCTTGCGAAGGTGGCGTTATCCGCGGCTCCCCTTGCTAGTCATCGGCTTTATTCCCTTTGTGTTGGTAGTGATGCTATtccgtc                <  2:943445/145‑1 (MQ=255)
                                                                                                                  ttaGTCATCGCCGTTTCAAAGGTGCGGATGCGCAAAATTCTGCGAGTAGGAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtccgtcc               <  2:275763/149‑1 (MQ=255)
                                                                                                                   taGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtccgtccg              <  1:440760/149‑1 (MQ=255)
                                                                                                                    aGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTTATGCTATTCCGTCCgg             >  2:278712/1‑149 (MQ=255)
                                                                                                                        aTCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTTCGTTCGTAGTGCTGCTATTCCGTCCAGCGCa         >  1:138070/1‑149 (MQ=255)
                                                                                                                        aTCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCa         >  2:850214/1‑149 (MQ=255)
                                                                                                                        aTCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCa         >  2:849573/1‑149 (MQ=255)
                                                                                                                        aTCGCCGTTTCAAAGTTGCAGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCa         >  1:851569/1‑149 (MQ=255)
                                                                                                                               tttCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTa  >  1:763231/1‑149 (MQ=255)
                                                                                                                               tttCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTa  >  1:790309/1‑149 (MQ=255)
                                                                                                                                    |                                                                                                                                               
GTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTA  >  NC_000913_3_pae_tpiA/2795264‑2795539

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