Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_vch_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_vch_tpiA2,795,3960CA100.0% 168.9 / NA 49Q382K (CAG→AAG) gabPgamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (24/25);  total (24/25)

CCCGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTA  >  NC_000913_3_vch_tpiA/2795251‑2795539
                                                                                                                                                 |                                                                                                                                               
cccGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGt                                                                                                                                              >  2:604377/1‑149 (MQ=255)
       gTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTATGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGAt                                                                                                                                       >  2:721608/1‑149 (MQ=255)
        ttgtggTGTTCCACTTCACCGGCGCGGTATTCTTTACGGTGTTGGTGAACTATTACGCGCCTGCGAAGTTAGTTAAATTTCTGATCGACTGCTCCGGCTATCTTGCTCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATg                                                                                                                                      <  2:191949/145‑1 (MQ=255)
        tGGCGGTGTTACTCTCCACCGGCGCGGCATTCCTAACGGTGTTGGTGAACTATTACGCGCCTGCGAAGGTATTTAACTTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCTGATg                                                                                                                                      <  2:192505/149‑1 (MQ=255)
         gggggTGTTACTCTCCCCCGGCGCGGCATTCTGAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATgc                                                                                                                                     <  1:345032/146‑1 (MQ=255)
           cGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTAGTTCAATTTCTGATCGACAGATGCGCCGCTATAGCCCTGCTGGTTTATTTAGTCATCGCCATTTCAAACTAGCGGATgcgc                                                                                                                                   >  2:505920/1‑149 (MQ=255)
                         aCCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGc                                                                                                                     >  1:78117/1‑149 (MQ=255)
                             gcgcGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaa                                                                                                                 <  2:757257/149‑1 (MQ=255)
                             gcgcGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaa                                                                                                                 <  1:602049/149‑1 (MQ=255)
                             gcgcGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaa                                                                                                                 <  1:461040/149‑1 (MQ=255)
                                        ttAACGGTGGTGGTGAACTATTACGCGCCGGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGCGCAGAAGGAAGCGAAAt                                                                                                      >  1:533292/1‑149 (MQ=255)
                                          aaCGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTc                                                                                                    >  2:44446/1‑149 (MQ=255)
                                          aaCGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTc                                                                                                    <  1:656384/149‑1 (MQ=255)
                                            cGGTGGTGGTGAACTATTACGGCCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGCAGGAAGCGAAATTCGc                                                                                                  <  1:726851/149‑1 (MQ=255)
                                             ggtggtggtGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCt                                                                                                 <  1:49373/149‑1 (MQ=255)
                                                  tggtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCgg                                                                                            >  1:582327/1‑148 (MQ=255)
                                                      gAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAAtgt                                                                                        >  1:779109/1‑149 (MQ=255)
                                                       aaCTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAAtgtg                                                                                       <  2:582327/149‑1 (MQ=255)
                                                         cTATTACGCGCCTGCGAAGTTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGc                                                                                     <  1:338042/149‑1 (MQ=255)
                                                          tATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTg                                                                                                                                            <  1:290203/93‑1 (MQ=255)
                                                          tATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTg                                                                                                                                            >  2:290203/1‑93 (MQ=255)
                                                          tATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCt                                                                                    >  1:531297/1‑149 (MQ=255)
                                                           aTTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGTGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGTGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCtt                                                                                   >  2:522040/1‑149 (MQ=255)
                                                              aCGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTAt                                                                                >  1:390311/1‑149 (MQ=255)
                                                                       cGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGAGAATGTGGCTTTATCCGTGGCtg                                                                       >  1:129274/1‑148 (MQ=255)
                                                                           ggTTTTTAGAGTTGTGAGCGACAGCGCCGGCGCTATCGCTCGGCTGGTTTATTTAGGCGTCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGGAGAAGGAAGCGAAAGTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCt                                                                   <  2:733913/149‑1 (MQ=255)
                                                                              agttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGc                                                                <  1:604377/147‑1 (MQ=255)
                                                                                  aaaTTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGt                                                            >  1:331407/1‑149 (MQ=255)
                                                                                      ttCTGGTCAACATCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATc                                                        <  2:338682/149‑1 (MQ=255)
                                                                                      ttCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATc                                                        <  2:331407/149‑1 (MQ=255)
                                                                                      ttCTAATCGACAGCTGCTTCGGTAGCGCCCGGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATc                                                        <  2:339423/149‑1 (MQ=255)
                                                                                         tGATCGACAGCTCCGGCGCTATCGCCATGCTGGTTTATTTAGTCATCGCCGTTTCAACGTTGCGGATGCGCAAACTTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCCTCGGc                                                     >  2:736039/1‑149 (MQ=255)
                                                                                            cctaCAGCTCCGGCCCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCttt                                                  <  1:605424/146‑1 (MQ=255)
                                                                                                     ccGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCttt                                         >  1:796588/1‑149 (MQ=255)
                                                                                                       ggCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGACTGTGGCTTGTGCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTtgt                                       >  1:482008/1‑149 (MQ=255)
                                                                                                       ggCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTtgt                                       >  1:428068/1‑149 (MQ=255)
                                                                                                            tATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATATGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTgg                                  >  1:692456/1‑149 (MQ=255)
                                                                                                             aTCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGATAGTCATCGGCATTATTACCCTtgattttgt                                 >  2:504648/1‑142 (MQ=255)
                                                                                                              tCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTa                                <  1:44446/149‑1 (MQ=255)
                                                                                                                gCCCTGCTGGTTTATTTAGTCATCGGCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGt                              <  1:745823/149‑1 (MQ=255)
                                                                                                                gCCCTGCTGGGTTATTTAGTCACCACCGTTTCAAAGTTGCGGATGTGCAAAATTGTTCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTTTTACCTTTGTGTTGGTAGt                              <  2:704067/149‑1 (MQ=255)
                                                                                                                        ggTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTAt                      <  1:170618/149‑1 (MQ=255)
                                                                                                                                  gTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGgc            <  2:259132/149‑1 (MQ=255)
                                                                                                                                     aTCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGACTGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCa         >  2:500680/1‑149 (MQ=255)
                                                                                                                                     aTCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTTGCTCACCTGGCTAGTCATCGGCTTTATTTCCTTTTCGTTGGTAATGATGCTATTCCGTCCGGCGCa         >  2:512507/1‑149 (MQ=255)
                                                                                                                                      tCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCaa        >  2:592503/1‑149 (MQ=255)
                                                                                                                                        gCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGTTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACa      <  1:8182/149‑1 (MQ=255)
                                                                                                                                        gCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACa      <  1:592503/149‑1 (MQ=255)
                                                                                                                                            tttCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGCTGCTATTCCGTCCGGCGCAACAGTTa  >  1:806987/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CCCGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTA  >  NC_000913_3_vch_tpiA/2795251‑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: ‑