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

GTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTA  >  NC_000913_3_hsa_tpiA/2795266‑2795539
                                                                                                                                  |                                                                                                                                                 
gtTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCaaaa                                                                                                                                 >  1:693142/1‑149 (MQ=255)
   aCTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGCACTATTACGCGCCTGCGACGGTATTTAAATGTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTCGATTATTTAGTCATCGCCGTTTCAAAGTTGCGGCGGCCCACAATTc                                                                                                                              >  2:85815/1‑149 (MQ=255)
   aCTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTc                                                                                                                              <  1:164970/149‑1 (MQ=255)
    ctctCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGAGTAATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTAAAAGTTTCGACTGCGCAAAATACt                                                                                                                             >  1:205916/1‑149 (MQ=255)
       tcCACCGGCGCCGCATTCTTAACGGTGGTTGTGAACGATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCg                                                                                                                         <  1:14810/149‑1 (MQ=255)
          aCCGGCGCGGCATTCTTAACGGTGGTGCTGCATTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGAGGCGCAAAATTCTgcgcgc                                                                                                                      >  2:600218/1‑149 (MQ=255)
                gcgGCATTCTTAACGGTGGTGGGGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaagg                                                                                                                <  2:627089/149‑1 (MQ=255)
                 cgGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaagga                                                                                                               <  1:130407/149‑1 (MQ=255)
                  ggCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaaggaa                                                                                                              <  1:74099/149‑1 (MQ=255)
                  ggCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaaggaa                                                                                                              <  1:595252/149‑1 (MQ=255)
                   gCATTCTTATCGGTGGTGGTGAACCATTACGCGCATGTGAAGGTATTTAAATTTCTGATCGACAGTTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCTTCGCCTTTTCAAAGTTGCGGATGCGCAAAATTCTGCTCGCAgaattttg                                                                                                             <  2:402350/149‑6 (MQ=255)
                   gCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaaggaag                                                                                                             >  1:693009/1‑149 (MQ=255)
                   gCAGTCTTAACGGTGGTCGTGAACTGTTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGTTGGTTTATTTAGTCATCGCCGTGTCAAAGTTGCTGATGCGCAAAATTCTGCGAGCAgaaggaag                                                                                                             <  2:52215/149‑1 (MQ=255)
                      ttCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGa                                                                                                                                          <  1:740035/118‑1 (MQ=255)
                      ttCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGa                                                                                                                                          >  2:740035/1‑118 (MQ=255)
                         ttAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAAt                                                                                                      >  1:260408/1‑149 (MQ=255)
                          tAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTCGTCATCTCCCCTTCAAAGTTGCGGATGGCCAAATCGTAGCTAGCAGAAGGACGCGAAAtt                                                                                                     >  1:189587/1‑149 (MQ=255)
                                 ggtggtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGAACGACATCTCCTGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGCTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAACTCCCTTGc                                                                                              >  2:4150/1‑149 (MQ=255)
                                      tgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGa                                                                                                                        <  1:458492/119‑1 (MQ=255)
                                      tgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGa                                                                                                                        >  2:458492/1‑119 (MQ=255)
                                          cTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGc                                                                                     >  1:780222/1‑149 (MQ=255)
                                          cTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGATGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTgga                                                                                     >  2:491881/1‑148 (MQ=255)
                                             ttACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCttt                                                                                  >  2:647708/1‑149 (MQ=255)
                                               aCGGGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTAt                                                                                <  1:795946/149‑1 (MQ=255)
                                                   gcCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGGTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAAAAGGAGCGCAAATCCCCTTGCGAATGTGGCTTTATCCGt                                                                            >  2:343597/1‑149 (MQ=255)
                                                     cTGCGAAGGTATTTAAATTTCTGATCGACCGCTCCGGCGCTATCGCCCTGCTGGCTTATTTAGTCATCGCTGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCCTGCGAATGTGGCTTTTTCCGTgg                                                                          >  1:704223/1‑149 (MQ=255)
                                                     cTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGAATAAGGTAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTgg                                                                          <  1:410948/149‑1 (MQ=255)
                                                                 ttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTAGTATTTCTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTa                                                              <  1:65391/149‑1 (MQ=255)
                                                                  tAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAg                                                             <  2:780222/149‑1 (MQ=255)
                                                                       ttCGGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATc                                                        <  1:631281/149‑1 (MQ=255)
                                                                          tGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCTGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGc                                                     >  1:743120/1‑149 (MQ=255)
                                                                           gATCGACCGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCt                                                    <  2:498028/149‑1 (MQ=255)
                                                                           gATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCt                                                    >  1:631262/1‑149 (MQ=255)
                                                                             tCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCttt                                                  >  2:602459/1‑149 (MQ=255)
                                                                               gACAGCTCCGGCGCTATCGCCCTGCTGGTTTCTTTAGTCATCGCCGTTCCAAACTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCGGTGGCTCACCTGGCTAGTCATCGGCTttat                                                >  2:452309/1‑149 (MQ=255)
                                                                               gACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTttat                                                >  2:415655/1‑149 (MQ=255)
                                                                                         gcgcTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTtgtg                                      <  2:693142/149‑1 (MQ=255)
                                                                                             tATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTgg                                  <  1:415655/149‑1 (MQ=255)
                                                                                             gaTCGCCCTGCTGGTTTATTTGATCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTgg                                  <  1:600218/148‑1 (MQ=255)
                                                                                                  cccTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCt                                                                        <  1:238325/106‑1 (MQ=255)
                                                                                                  cccTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCt                                                                        >  2:238325/1‑106 (MQ=255)
                                                                                                   ccTGCTGCTTTATTTAGTCATCGCCGTTTCAAACTTGCGCATGCGCAAAATTCTGCTAGCAGAAGGAATCGAAATTCGCTTGCGAATTTTGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGa                            <  1:491881/149‑1 (MQ=255)
                                                                                                    cctcttatttatttaGTCACCGCCGTTTCAAAGTTGCTAAATCCCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCTCTTGCCAATTTGTCTTTATCCTTAGCTCACCTAGCTAGTCATCTGCTTTATTACCTTTGTCTTAGTAGTCAt                           <  2:777927/142‑1 (MQ=255)
                                                                                                        tgGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTGGTCATCGCCTCTATTACCTTTGTGTTGGTCGTGATGCtc                       >  2:679879/1‑148 (MQ=255)
                                                                                                          gTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTCATTACCTATCGGGTGGTATTGATGATAtt                     >  2:440100/1‑149 (MQ=255)
                                                                                                             tagttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCGGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGATAGTGAGTTGCTTTATTCCCTTTGTGTTGGTAGTGATGCTATtccg                  <  1:452309/146‑1 (MQ=255)
                                                                                                                   gTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGgc            >  1:673165/1‑149 (MQ=255)
                                                                                                                    tCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGGGATGCTATTCCGTCCGgcg           >  1:132928/1‑149 (MQ=255)
                                                                                                                          ccGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAg     >  1:668407/1‑149 (MQ=255)
                                                                                                                           cGTTTCAAAGTTGCGGCTGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGt    <  1:647708/149‑1 (MQ=255)
                                                                                                                             tttCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTa  <  1:89310/149‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                                 
GTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTA  >  NC_000913_3_hsa_tpiA/2795266‑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: