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% 140.2 / NA 42Q382K (CAG→AAG) gabPgamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (25/17);  total (25/17)

CCCGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTT  >  NC_000913_3_vch_tpiA/2795251‑2795538
                                                                                                                                                 |                                                                                                                                              
cccGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGt                                                                                                                                             >  1:610975/1‑149 (MQ=255)
                      tcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTggg                                                                                                                       >  1:819994/1‑147 (MQ=255)
                      tcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCg                                                                                                                       >  2:793204/1‑149 (MQ=255)
                                cgGCATTCTTAACGGTGGTGGGGAACTATGTCGCGCCCGCGAAGGTATTTAAAGTTCTGATCGACAGCTCCGGCGCTAGCGCTCGGCTGGTTTATTTAGTCATGGCCGTTTCAAAGGTGCGAATGCGCAAAATTCTGCGAGCAgaagga                                                                                                             <  2:298142/149‑1 (MQ=255)
                                 ggCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaaggaa                                                                                                            >  2:819773/1‑149 (MQ=255)
                                 ggCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaaggaa                                                                                                            <  1:793204/149‑1 (MQ=255)
                                  gCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATATAAATTTCTGATCGACAGCATCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaaggaag                                                                                                           >  1:760362/1‑149 (MQ=255)
                                   cATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGc                                                                                                          >  2:464442/1‑149 (MQ=255)
                                   cATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGc                                                                                                          >  2:262268/1‑149 (MQ=255)
                                    aTTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTGTTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCg                                                                                                         >  2:679881/1‑149 (MQ=255)
                                    aTTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTAATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGTAAGCg                                                                                                         >  2:679803/1‑149 (MQ=255)
                                           aCGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCg                                                                                                  <  2:260211/149‑1 (MQ=255)
                                            cGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGc                                                                                                 >  1:766878/1‑149 (MQ=255)
                                              gtggtggtgATCTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTACTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCtt                                                                                               <  1:494499/149‑1 (MQ=255)
                                                 gtggtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCCAAGTTGCGGATGCGCAAAATTCTGCGGGCAGAAGGAAGCGAACTTCGCTTGCg                                                                                            >  1:427727/1‑149 (MQ=255)
                                                  tggtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCgg                                                                                           >  1:265793/1‑148 (MQ=255)
                                                    gtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAAt                                                                                         <  1:112309/149‑1 (MQ=255)
                                                    gtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTCGCGAAt                                                                                         >  1:125973/1‑149 (MQ=255)
                                                      gAACTATTACGCGCATGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCATAAGGAAGCGAAATTCGCTTGCGAAtgt                                                                                       <  2:829686/149‑1 (MQ=255)
                                                        aCTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAAGGTgg                                                                                     >  1:139698/1‑149 (MQ=255)
                                                         cTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGTCAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGGCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGc                                                                                    <  1:13806/149‑1 (MQ=255)
                                                         cTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGc                                                                                    <  2:567932/149‑1 (MQ=255)
                                                          tATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCt                                                                                   <  1:819773/149‑1 (MQ=255)
                                                              aCGCGCCTGCGAAGGTATTTAAACTACTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTTAAAGTTGCGGATGCGCAAAACTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTAt                                                                               <  1:761127/149‑1 (MQ=255)
                                                                     tGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGCAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGc                                                                        >  1:523311/1‑149 (MQ=255)
                                                                      gCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTGCAAAGTTGAGCTCCCGGATGATTGTTCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCt                                                                       >  1:195122/1‑149 (MQ=255)
                                                                       cGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTc                                                                      >  1:100817/1‑149 (MQ=255)
                                                                       cGAAGGTATTTAAATTGCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTc                                                                      <  2:294231/149‑1 (MQ=255)
                                                                        gAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGAAGAAGTAAGTGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCa                                                                     >  1:578494/1‑149 (MQ=255)
                                                                         aaGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCAc                                                                    >  2:587960/1‑149 (MQ=255)
                                                                           ggTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCt                                                                  >  1:441904/1‑149 (MQ=255)
                                                                                     tttCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTCGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCAt                                                        <  2:125973/149‑1 (MQ=255)
                                                                                                  gCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTAcc                                           <  1:314128/149‑1 (MQ=255)
                                                                                                    tCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTACCATGGCTAGTCTACGGCTTTATTACCtt                                         >  1:570186/1‑149 (MQ=255)
                                                                                                                gCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGt                             >  1:566970/1‑149 (MQ=255)
                                                                                                                  ccTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGa                           <  1:464442/149‑1 (MQ=255)
                                                                                                                     gctgGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATGCGTGGATCACCTGGCTAGTCACCGGCTTTATTGCCTTTGTGCTGGTGGTGTTGc                        >  2:192007/1‑149 (MQ=255)
                                                                                                                       tgGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGGGCAGAAGCAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTTTTAGTGATGCTa                      >  2:693674/1‑149 (MQ=255)
                                                                                                                         gcttatttaCTCATCGCCGTTTCAAAGTTCCGGATGCGCAAAATTCTGCCAGCAAAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTACCCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTAtt                    <  2:495081/147‑1 (MQ=255)
                                                                                                                              tttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtccgtc               <  1:262268/149‑1 (MQ=255)
                                                                                                                                     aTCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTGTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGGGTTGGTAGTGAGGCTATTCCGCCCGGCGCa        >  1:611043/1‑149 (MQ=255)
                                                                                                                                           gTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGtt  <  1:587960/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CCCGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTT  >  NC_000913_3_vch_tpiA/2795251‑2795538

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