Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 2,817,170 C→A E112* (GAA→TAA)  yqaB ← fructose‑1‑P and 6‑phosphogluconate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA2,817,1700CA97.8% 140.0 / NA 45E112* (GAA→TAA) yqaBfructose‑1‑P and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (20/24);  minor base G (1/0);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGA‑AGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCC  >  NC_000913_3_hsa_tpiA/2817026‑2817315
                                                                                                                                                |                                                                                                                                                  
tCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTc                                                                                                                                                >  1:293626/1‑149 (MQ=255)
tCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTc                                                                                                                                                >  1:446094/1‑149 (MQ=255)
tCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGAAGAGATGGGCTGGGACGACGGCGTCAAAATAAGAGCGCAATCCCAGGTGCCCCGGCAATGCCTCCGCGCTGTCGCTTTGACTc                                                                                                                                                >  2:41622/1‑149 (MQ=255)
   gctGCACGCCCATACGCTGCGCGCAAAACAAAAATGTGTCTGGCGCGGGGTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGAGACAGGTAGGCCGCAGCGATGGCGCTTTAACTcccc                                                                                                                                             >  2:697747/3‑149 (MQ=255)
        aCGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCTTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTcc                                                                                                                                        <  2:247982/149‑1 (MQ=255)
        aCGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTcc                                                                                                                                        <  1:520745/149‑1 (MQ=255)
          gCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTa                                                                                                                                      <  2:647526/149‑1 (MQ=255)
            ccATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGAAATGCCCCAGCGATGGCGCTTTAACTCCCTGTTCATACa                                                                                                                                    >  1:12205/1‑149 (MQ=255)
                aCGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCa                                                                                                                                <  1:411622/149‑1 (MQ=255)
                   cTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTg                                                                                                                             <  1:601897/149‑1 (MQ=255)
                         gcACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGcgac                                                                                                                       >  2:723138/1‑149 (MQ=255)
                             aacaaAAATGTGTCTGGCGCGGGTTTAAAGAGTCAGAGGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGAcc                                                                                                                   >  2:525736/1‑149 (MQ=255)
                                   aaTGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCa                                                                                                             >  2:66838/1‑149 (MQ=255)
                                    aTGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCaa                                                                                                            <  2:276918/149‑1 (MQ=255)
                                         tCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACtttt                                                                                                       >  2:626860/1‑149 (MQ=255)
                                           tGGCGCGGGTTTACGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTca                                                                                                     <  1:708706/149‑1 (MQ=255)
                                            ggCGCGGGTTTATGTTTTTTGACGTGATCAGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTcac                                                                                                    <  2:212160/149‑1 (MQ=255)
                                                 gggTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACAt                                                                                               <  1:114057/149‑1 (MQ=255)
                                                   gTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTcacc                                                                                                   <  2:599517/143‑1 (MQ=255)
                                                   gTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTcacc                                                                                                   >  1:599517/1‑143 (MQ=255)
                                                     ttATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATcaac                                                                                           <  1:626860/149‑1 (MQ=255)
                                                           tgtTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAgg                                                                                     >  1:578442/1‑149 (MQ=255)
                                                                gACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGA‑AGCg                                                                               <  2:723950/149‑1 (MQ=255)
                                                                gACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGA‑AGCg                                                                               <  1:401368/149‑1 (MQ=255)
                                                                          gCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACCACAGCA‑AGCGGTTCGAcaca                                                                     >  2:98550/1‑146 (MQ=255)
                                                                           cAGCGACGACGGCGTCAAAATACTGGCGTAATATCAGGAGTGATAGCAATGCCTCAGCGATGGCGCTTTACCTCCCCGTTCCTACAGCCATTGGGCGACGCCCATGCCAACTTTTCACCACATCAACAAGCGGA‑AGCGGTTCCCCGCTa                                                                    >  2:492897/1‑149 (MQ=255)
                                                                            aGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGTCGACAATGCCCACTTTTCACCACATCAACAA‑AGGATACCGGTTCAACGCTAt                                                                   >  2:588042/1‑149 (MQ=255)
                                                                                  gacgGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCTCCACATCAACAAGAGGA‑AGCGGTTCGACGCTATCcagca                                                             <  2:524830/149‑1 (MQ=255)
                                                                                   acgGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGA‑AGCGGTTCGACGCTATCcagcag                                                            <  1:723138/149‑1 (MQ=255)
                                                                                           aaaaTAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGA‑AGCGGTTCGACGCTATCCAGCAGCATActtc                                                    <  2:693950/149‑1 (MQ=255)
                                                                                              ataatGGCGTAATCCCAGGTGCGCCAGCAATGCCTAAGGGAGCGCTCTTTAACGCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGA‑AGCGGTTCGACGCTATCCAGCAGCATACTTCTTa                                                 <  2:446094/149‑1 (MQ=255)
                                                                                              agaatGGCGTACTCCCAGCTGGGCGATCGATGCGCCAGCGAAGGCGCGTCAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGA‑AGCGGTTCGACGCTATCCAGCAGCATACTTCTTa                                                 <  1:697747/147‑1 (MQ=255)
                                                                                               taatGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGA‑AGCGGTTCGACGCTATCCAGCAGCATACTTCTTAc                                                <  2:612809/149‑1 (MQ=255)
                                                                                                          tCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGA‑AGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGttt                                     >  1:626353/1‑149 (MQ=255)
                                                                                                          tCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGA‑AGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGttt                                     <  2:578442/149‑1 (MQ=255)
                                                                                                            ccAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGA‑AGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGttttt                                   <  2:293626/149‑1 (MQ=255)
                                                                                                             cAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGTA‑AGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGtttttt                                  <  1:336075/149‑1 (MQ=255)
                                                                                                             cAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGA‑AGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGcttttt                                  >  2:188922/1‑149 (MQ=255)
                                                                                                                  gcgcCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGA‑AGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCAcgc                             <  1:658862/149‑1 (MQ=255)
                                                                                                                       aGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTCCAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGA‑AGCGGTTCGACGCTATCCAGCACCATACTTCTTCCTGCTTCTGTTTTTTCCCGCGCCaa                        >  2:357750/1‑149 (MQ=255)
                                                                                                                               cTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTCCACCACATCAACAAGAGGA‑AGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTACTGTTTTTTAACACGCTAAAGCATgag                >  1:64363/1‑147 (MQ=255)
                                                                                                                                tCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGA‑AGCGGTTCGACGCTATCCAGCAGCATActtc                                                    <  1:351745/112‑1 (MQ=255)
                                                                                                                                tCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGA‑AGCGGTTCGACGCTATCCAGCAGCATActtc                                                    >  2:351745/1‑112 (MQ=255)
                                                                                                                                      aTGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGA‑AGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTGgag         >  2:92667/1‑149 (MQ=255)
                                                                                                                                             tttAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGA‑AGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGcc  >  2:187258/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGA‑AGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCC  >  NC_000913_3_hsa_tpiA/2817026‑2817315

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