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.2% 161.3 / NA 54Q382K (CAG→AAG) gabPgamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (26/27);  total (27/27)
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

CCGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTA  >  NC_000913_3_pae_tpiA/2795252‑2795539
                                                                                                                                                |                                                                                                                                               
ccGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGtt                                                                                                                                             >  2:746375/1‑149 (MQ=255)
 cGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTg                                                                                                                                            >  1:772760/1‑149 (MQ=255)
   tACGTGGCGGTGTTACTCTCCACCGGCGCGGCACTCTGATCTGTGGTTGTGAACTAGTGTGCGCCTGCGCAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCTCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCg                                                                                                                                          <  2:521239/149‑1 (MQ=255)
   tACCTGGCCGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGAACGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCg                                                                                                                                          <  1:368458/149‑1 (MQ=255)
              gtTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCCCCGTTTCAAAGTTGCGGATGCGCaaaa                                                                                                                               >  1:76979/1‑149 (MQ=255)
                   tctcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTg                                                                                                                          >  1:270586/1‑149 (MQ=255)
                    ctcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGc                                                                                                                         <  1:935645/149‑1 (MQ=255)
                            gcgcGGCATTCTGAACGGTGGTGGTGAACTATTACGCGCCTGCGTAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaa                                                                                                                 <  1:789844/149‑1 (MQ=255)
                                ggCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAATCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgatgaaa                                                                                                             >  1:449860/1‑144 (MQ=255)
                                ggCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGCCGCTATCGCCCTGCTGGTTTATTTAGTAATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGCGCAgaaggaa                                                                                                             >  1:549520/1‑149 (MQ=255)
                                 gCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTATGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTATCAAAGTTCCGCCGGCTCACAATTCAGCAGGAAgaaggaaa                                                                                                            >  2:195170/1‑148 (MQ=255)
                                         aaCGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACACCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTTCGGATGCCCAAAATTCTGCGAGCAGAAGGAAGCGAAATTc                                                                                                    >  2:294554/1‑149 (MQ=255)
                                          aCGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCATTTCAAAGTTGCGGATGCGCCAAATTCTGCGAGCAGAAGGAAGCGAAATTCg                                                                                                   >  2:130544/1‑149 (MQ=255)
                                           cGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCAGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCACGCAGAAGGAAGCGAAATTCGc                                                                                                  >  1:498487/1‑149 (MQ=255)
                                            ggtggtggtGAACTATTACGCGCCTGCGAAGGTATTTAAAGTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCt                                                                                                 >  2:1008655/1‑149 (MQ=255)
                                              tggtggtgAACTGTTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGAt                                                                                                                                       >  2:646660/1‑109 (MQ=255)
                                              tggtggtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGAt                                                                                                                                       <  1:646660/109‑1 (MQ=255)
                                                     gAACTATTACGCTCCTGCGAAGGTATTTAAATTTCTGATCTACAGCTCCGGCGCTATCGCCCTTCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAAtgt                                                                                        <  1:790205/149‑1 (MQ=255)
                                                      aaCGATTACTCGCCTGCTAAGGTATTTAAATTTCTGATCGCCAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAAtgtg                                                                                       <  2:141372/149‑1 (MQ=255)
                                                           ttACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCttt                                                                                  >  1:526671/1‑149 (MQ=255)
                                                                   cTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGaa                                                                                           <  2:437697/132‑1 (MQ=255)
                                                                   cTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGaa                                                                                           >  1:437697/1‑132 (MQ=255)
                                                                   cTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTgg                                                                          >  1:869723/1‑149 (MQ=255)
                                                                        aaGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGACGCACAATTAGCTTGAGAATGTCGCTTTATCCGTGGCTccc                                                                     >  2:821315/1‑147 (MQ=255)
                                                                             aTTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATTCGCAAAATTCTGTGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGc                                                                <  1:555677/149‑1 (MQ=255)
                                                                              tttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCt                                                               >  1:1031447/1‑149 (MQ=255)
                                                                               ttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAATTTCGCGTGCGACTGTGGCTTTATCCGTGGCTCACCTGGCTa                                                              <  2:270586/149‑1 (MQ=255)
                                                                               ttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTa                                                              <  2:1024540/149‑1 (MQ=255)
                                                                               tcaaaTTTCTGATATACACCTCAGGCGCTATGGTCCTGCGCTTTGATTTAGTCAGCGCCTTTTCAAAGTGGCGGATGCGCAAAACTCTGCGAGCAGAAGGAAGCGAAAGTCGCTTGCGAATGTGGCTTTACCCGTGGCTCACCTGGCTa                                                              <  2:22766/147‑1 (MQ=255)
                                                                                      tctcaacGACAGCTCCGGCGCCATCGCCCTGCTGGTTTATTTAGTCATCGCCGTATCAAAGTTGCGGATGCGCAAAATTCTGCCAGCAGAAGGAATCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCg                                                       <  2:790069/143‑1 (MQ=255)
                                                                                         gATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCAGTTTCAAAGTTACGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCCCACATGGCTACTCATCGCCt                                                    >  2:197648/1‑149 (MQ=255)
                                                                                            cGACAGTCCCGGCGCTATCGCCCTGCGGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAAGTCGCGTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTtta                                                 <  2:48729/149‑1 (MQ=255)
                                                                                              aCAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGc                                                                                                  <  2:378894/98‑1 (MQ=255)
                                                                                              aCAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGc                                                                                                  >  1:378894/1‑98 (MQ=255)
                                                                                              aCAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAATCATCGGCTttatt                                               >  1:936807/1‑149 (MQ=255)
                                                                                               cAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTttatta                                              <  1:1008655/149‑1 (MQ=255)
                                                                                                    ccGGCGCTATCGGCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCttt                                         <  2:542856/149‑1 (MQ=255)
                                                                                                       gcgcTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGATTTATCCGTGGCTCACCTGGCTAGTCCTCGGCTTTATTAACTTtgtg                                      >  1:65382/1‑149 (MQ=255)
                                                                                                          ctactgtccTGCTTGTTTATTTAGTCATCGTCGTTTCAAAGTTGCGCATGCGCAAAATTCTGCGAGCAGAAGGAAGACCACTTCACTATCGCCTGGTTGGGCATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTg                                   <  2:119799/142‑1 (MQ=255)
                                                                                                               gCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACATTTGAGCTGCTAGt                              >  1:710178/1‑149 (MQ=255)
                                                                                                                 ccTGCTTGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGa                            <  1:134027/149‑1 (MQ=255)
                                                                                                                  ctgctgGTTTATTTAGTCATCGCCGGTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGAt                           <  2:96445/149‑1 (MQ=255)
                                                                                                                   tgctgGTTTATTTAGTCATCGACGTTTAACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGAGAAATTCGCTTGCGCATATCGCTCACTCCGTTGCTCACCTGCCTAGTCAACGGCTTTATACCCTTGGTGTTGGCAGtagag                          >  2:425214/1‑145 (MQ=255)
                                                                                                                             attatgcATCGCTGTTTCAAAGTGGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtccgtc                <  2:76979/143‑1 (MQ=255)
                                                                                                                              ttaGTCATCGCCTTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtccgtcc               <  1:856580/149‑1 (MQ=255)
                                                                                                                                aGTCATCGCCGTTTCAAAGTTGCCGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCgg             <  1:625411/149‑1 (MQ=255)
                                                                                                                                 gTCATCGCCTTTTCAAAGTAGCTGATGCGCAAAATTCTGCGAGCCGGAGTTATCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCGGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGgc            <  1:140425/149‑1 (MQ=255)
                                                                                                                                 gTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGgc            <  1:128162/149‑1 (MQ=255)
                                                                                                                                   cATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGgcgc          <  1:800327/149‑1 (MQ=255)
                                                                                                                                     tCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCaa        >  2:466859/1‑149 (MQ=255)
                                                                                                                                        ccGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAg     >  1:347555/1‑149 (MQ=255)
                                                                                                                                        ccGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGATGGTAGTCATGCTATTCCGTCCGGCGCAACAg     >  2:75428/1‑149 (MQ=255)
                                                                                                                                         cGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACATTTGTGTTGGTAGTGATGCTATTCCGGCCGGCGCAACAGt    <  2:324614/149‑1 (MQ=255)
                                                                                                                                           tttCAAAGTTGCGGATGCGCAAAATTCTGCGAGCACACGTCAGCGAAATTCGCTTGCGAATGTGCCTTTATCCGTCGCTCACCTGGCTAGTCATCGGCCTCATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGGCACCGTTa  >  2:270324/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                               
CCGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTA  >  NC_000913_3_pae_tpiA/2795252‑2795539

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑