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% 174.5 / NA 54Q382K (CAG→AAG) gabPgamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (24/29);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACA  >  NC_000913_3_pae_tpiA/2795252‑2795535
                                                                                                                                                |                                                                                                                                           
ccGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGtt                                                                                                                                         >  2:808307/1‑149 (MQ=255)
 cGTACGGGGCCGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTg                                                                                                                                        <  2:181670/149‑1 (MQ=255)
  ggaCGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACTGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGc                                                                                                                                       <  2:687763/147‑1 (MQ=255)
  gTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGc                                                                                                                                       <  1:473084/149‑1 (MQ=255)
      gTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGAt                                                                                                                                   >  2:743950/1‑149 (MQ=255)
      gTGGCGCTGTTACTCTCCACCGGCGCGGCATACTTAACGGTGATGGTGACCTATTACGCGCCTGCGTAGCAACTCAAAGCCCTGATCGACAGCTCCGGCCCTATCCACCTGCTGGTTTATTAAGTGATCACCGTTTCAAACTTCCGGAt                                                                                                                                   >  2:670587/1‑149 (MQ=255)
       tGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATg                                                                                                                                  >  2:961233/1‑149 (MQ=255)
              tttACGCTCCCCCGGTCCGTCATTTGTAACGGTTGTGGTGAACTATTACGCCCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCGCTGCTGGTTTATTTATTCATCGCCGTTTCAAAGTTGCGGATGCGCaaaa                                                                                                                           <  2:711916/148‑1 (MQ=255)
                      ccACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGa                                                                                                                   <  1:997337/149‑1 (MQ=255)
                       cACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATAGACAGCTCCATCGCTATCGTCCTGCTGGTTTATTTCGTCATCGCCGTTTCAAAGTTGGGGATGCGCAAAATTCTGCGAg                                                                                                                  >  2:869840/1‑149 (MQ=255)
                           ggCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAga                                                                                                              <  2:27979/149‑1 (MQ=255)
                            gcgcGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaa                                                                                                             <  1:678825/149‑1 (MQ=255)
                                ggCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaaggaa                                                                                                         >  2:126452/1‑149 (MQ=255)
                                             gtggtggtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCtt                                                                                            <  2:824437/149‑1 (MQ=255)
                                                   gtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTGTTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTCCGAAt                                                                                      >  2:144410/1‑149 (MQ=255)
                                                   gtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCGCCGGCGCTATCAACCTGCTGGTTTGTTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAGGGAAGCGAAATTCACttcggagt                                                                                      >  2:298208/1‑143 (MQ=255)
                                                    tgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGAAAGCGAAGTTCGCTTGCGCAtg                                                                                     >  2:160002/1‑149 (MQ=255)
                                                       aCTATTAGGCTCCTGCGAAGGTATTTAAATGTCTGATCGACAGCTCCGGCGCTATCGCCCTGATGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTGTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTgg                                                                                  <  2:426061/149‑1 (MQ=255)
                                                           ttACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATAGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCttt                                                                              >  1:541575/1‑149 (MQ=255)
                                                                cgcCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCg                                                                         >  1:363613/1‑149 (MQ=255)
                                                                 tccTGCGAACGTATTTAAATTTCTGATCTACAGCTCCGGCGCTATCGTCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAATTTGCGGATGCTCAAAATTCTGCTAGCCGAAGTAAGCGAAATTCGCTTGCGAATGTGGCTTTAGCCgc                                                                        <  2:459482/148‑2 (MQ=255)
                                                                 gcCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGt                                                                        <  2:569038/149‑1 (MQ=255)
                                                                  tcggCGAAGGTATTTTGATTTCTGACCGACAGGCCCGGCGCTATCGCCCTGCTGGGTTATTTAGGCCTCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCCGCAGCAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTg                                                                       <  2:231897/146‑1 (MQ=255)
                                                                   cTGCGAACGTATTTAAATTTCTGATCGACAGCGCCCGCCCCATCGCCCTGCTAGTTTATTAACTCCTTGTCGATTGAAAATTTCCTATTTGTAAAATTCTTTGAGCAAAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTgg                                                                      <  1:201336/149‑1 (MQ=255)
                                                                         aGGTATTTAAATTAATGATCGACAGCTAAGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTTCGGATTTTGCTTTATCCGTGGCGCAcc                                                                >  2:967216/1‑149 (MQ=255)
                                                                          ggTATTTAAATGTCTGATCGACAGCTCCGGCGCTATCGCCGTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCt                                                               <  2:691638/149‑1 (MQ=255)
                                                                               ttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTa                                                          <  1:808307/149‑1 (MQ=255)
                                                                                         gATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCt                                                >  1:559537/1‑149 (MQ=255)
                                                                                           gcttcatcTCCGGCGCATTCGCCCTGTTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGATGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCttt                                              <  2:504791/142‑1 (MQ=255)
                                                                                                  cTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTACATTACCt                                       >  2:102938/1‑149 (MQ=255)
                                                                                                   tCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCtt                                      >  1:870924/1‑149 (MQ=255)
                                                                                                         gcTATCGCCCTGCTGGTTTATTTAGTCAGCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTTCGAGCAGAATGAAGTGTAATTTTTTGGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGtt                                <  1:141832/149‑1 (MQ=255)
                                                                                                          ttATGGCCCAGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGTCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTg                               <  2:658703/148‑1 (MQ=255)
                                                                                                          cTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTg                               <  1:806491/149‑1 (MQ=255)
                                                                                                           tATCGTCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGCAGGAGTAAAATTCTGCGATCAGAAGGAAGCGAAATTCTCTTGTGAATGTGGCTTTATCCGTGGCTCACCTGGCTGGTCATCGGCTTTATTACCTTTGTGTTgg                              <  1:390284/149‑1 (MQ=255)
                                                                                                            aTCGCCCCGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTggc                             >  2:502922/1‑148 (MQ=255)
                                                                                                             tCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTa                            <  1:961233/149‑1 (MQ=255)
                                                                                                              cGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAg                           <  2:527455/149‑1 (MQ=255)
                                                                                                                  cttctgGTTTACTTACGCATCGCCGGTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGAt                       <  2:397290/146‑1 (MQ=255)
                                                                                                                    gctgGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGc                     <  2:541575/149‑1 (MQ=255)
                                                                                                                       ggTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTAt                  <  1:730905/149‑1 (MQ=255)
                                                                                                                        gTTTATTTAGTCATTGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGTCTCGTCATCTGCTTTATTACCTGTGTATTGGTAGTGATGCTAtt                 >  2:376807/1‑149 (MQ=255)
                                                                                                                        gTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCt                                                                                                                       <  1:193604/47‑1 (MQ=255)
                                                                                                                        gTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCt                                                                                                                       >  2:193604/1‑47 (MQ=255)
                                                                                                                        gTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCCGCCAGCCGAACGACGCGCAATTCCCTTGCGCCTGTCGCCTCATCCCTGGCCCACTTGGCCACTCTTCGGCTCCACTATCTGTGGGTTGGTAGGGATGCTAtt                 >  2:595593/1‑149 (MQ=255)
                                                                                                                         tttatttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtc                >  1:132968/1‑149 (MQ=255)
                                                                                                                          ttatttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtcc               <  1:517939/149‑1 (MQ=255)
                                                                                                                           tatttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTTTtcct              >  2:496289/1‑148 (MQ=255)
                                                                                                                                 gTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCCTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGgc        >  2:261002/1‑149 (MQ=255)
                                                                                                                                   cATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGgcgc      <  1:743950/149‑1 (MQ=255)
                                                                                                                                    aTCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTATCGTTGTGTTGGTAGTGTTGCTATTCCGTCCGGCGCa     >  2:107126/1‑149 (MQ=255)
                                                                                                                                       gCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACa  >  2:137721/1‑149 (MQ=255)
                                                                                                                                           tttCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGa                        <  1:831327/123‑1 (MQ=255)
                                                                                                                                           tttCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGa                        >  2:831327/1‑123 (MQ=255)
                                                                                                                                                |                                                                                                                                           
CCGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACA  >  NC_000913_3_pae_tpiA/2795252‑2795535

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