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,3960CA92.9% 179.8 / 6.8 56Q382K (CAG→AAG) gabPgamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (26/26);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

CCCGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTA  >  NC_000913_3_hsa_tpiA/2795251‑2795539
                                                                                                                                                 |                                                                                                                                               
cccGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGtt                                                                                                                                             >  1:548831/1‑150 (MQ=255)
cccGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGATTTATTTAGTCACCGCCGTTTCAAAGtc                                                                                                                                             >  2:350832/1‑149 (MQ=255)
                ttACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAAtt                                                                                                                             >  4:2496/1‑150 (MQ=255)
                ttACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGAAGCGCAAAAtt                                                                                                                             <  4:70912/150‑1 (MQ=255)
                    tctcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGc                                                                                                                         >  1:510817/1‑150 (MQ=255)
                    tcccccccGGCGCGGCATTCTGAACGGTGGTGGTGAACTATTCCGCGCCTGCGGAGGTATTTAAATTTGTGCTCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGc                                                                                                                         <  1:350832/144‑1 (MQ=255)
                      tcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGa                                                                                                                       >  3:459614/1‑150 (MQ=255)
                       ccACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTCTTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAg                                                                                                                      <  1:401448/150‑1 (MQ=255)
                         aCCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCa                                                                                                                    >  3:340491/1‑150 (MQ=255)
                              cgcgGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaagg                                                                                                               <  3:2496/150‑1 (MQ=255)
                                 ggCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAgaaggaag                                                                                                            >  4:139645/1‑150 (MQ=255)
                                        ttAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAAtt                                                                                                     <  4:120566/150‑1 (MQ=255)
                                            cGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCt                                                                                                 <  2:699447/150‑1 (MQ=255)
                                               tggtggtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGc                                                                                              >  4:231943/1‑150 (MQ=255)
                                                 cctatgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGa                                                                                            <  3:165886/146‑1 (MQ=255)
                                                       aaCTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTgg                                                                                      <  4:44073/150‑1 (MQ=255)
                                                                 cgcCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGt                                                                            <  3:231943/150‑1 (MQ=255)
                                                                  gcCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTg                                                                           <  4:190578/150‑1 (MQ=255)
                                                                   ccTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTgg                                                                          <  2:548831/150‑1 (MQ=255)
                                                                     tGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTgg                                                                          >  1:206382/1‑148 (MQ=255)
                                                                     tGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTgg                                                                          <  2:206382/148‑1 (MQ=255)
                                                                        gAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTg                                                                                                                          <  2:328595/97‑1 (MQ=255)
                                                                        gAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTg                                                                                                                          >  1:328595/1‑97 (MQ=255)
                                                                             tATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCaaa                                                                                                                                >  4:225850/1‑86 (MQ=255)
                                                                             tATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCaaa                                                                                                                                <  3:225850/86‑1 (MQ=255)
                                                                             tATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGc                                                                <  1:193427/150‑1 (MQ=255)
                                                                               tttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTa                                                              <  1:302934/150‑1 (MQ=255)
                                                                                   aaTTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGc                                                                >  2:500209/1‑144 (MQ=255)
                                                                                   aaTTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGc                                                                <  1:500209/144‑1 (MQ=255)
                                                                                    aTTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCACAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCAt                                                         >  4:23965/1‑150 (MQ=255)
                                                                                    aTTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCAt                                                         >  3:29289/1‑150 (MQ=255)
                                                                                        cTGATCGACAGCTCCGGCGCTATGGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGc                                                     <  4:266439/150‑1 (MQ=255)
                                                                                         tGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCt                                                    >  3:335490/1‑150 (MQ=255)
                                                                                              gACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTttatt                                               >  3:238491/1‑150 (MQ=255)
                                                                                                   cTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCtt                                          <  2:720804/150‑1 (MQ=255)
                                                                                                    tCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCttt                                         <  3:139645/150‑1 (MQ=255)
                                                                                                    tCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCttt                                         <  2:510817/150‑1 (MQ=255)
                                                                                                       ggCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTtgtg                                      >  2:347342/1‑150 (MQ=255)
                                                                                                         cgcTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGtt                                    >  2:210811/1‑150 (MQ=255)
                                                                                                          gcTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTg                                   <  3:353592/150‑1 (MQ=255)
                                                                                                           cTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTgg                                  >  1:154569/1‑150 (MQ=255)
                                                                                                             aTCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTa                                <  2:437437/150‑1 (MQ=255)
                                                                                                                 cccTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGa                            >  4:27703/1‑150 (MQ=255)
                                                                                                                        ggTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTAtt                     <  4:340491/150‑1 (MQ=255)
                                                                                                                        ggTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTAtt                     >  4:164803/1‑150 (MQ=255)
                                                                                                                        ggTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGCGATGCTAtt                     <  2:154569/150‑1 (MQ=255)
                                                                                                                           ttatttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTa                                                              <  1:148985/106‑1 (MQ=255)
                                                                                                                           ttatttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTa                                                              >  2:148985/1‑106 (MQ=255)
                                                                                                                            tatttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtccgt                 <  1:535884/150‑1 (MQ=255)
                                                                                                                               ttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtccgtccg              >  2:309798/1‑150 (MQ=255)
                                                                                                                               ttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtccgtccg              <  3:164803/150‑1 (MQ=255)
                                                                                                                                 aGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTTGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGgc            >  3:76938/1‑150 (MQ=255)
                                                                                                                                   tCATCGCCGTTTCAAAGTCGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGgcgc          >  3:55244/1‑150 (MQ=255)
                                                                                                                                      tCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAAc       >  2:636207/1‑150 (MQ=255)
                                                                                                                                         ccGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGt    >  2:356931/1‑150 (MQ=255)
                                                                                                                                           gcttCAAAATTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCTAAATTCCCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCCTTGTGTTGGTAGTGAAGCTATTCCGTCCGGCGCAACAGTTa  >  2:639926/3‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CCCGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTA  >  NC_000913_3_hsa_tpiA/2795251‑2795539

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