Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_tpiA 2,817,170 C→A 100% 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,1700CA100.0% 96.8 / NA 33E112* (GAA→TAA) yqaBfructose‑1‑P and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (12/21);  total (12/21)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCC  >  NC_000913_3_hsa_tpiA/2817026‑2817315
                                                                                                                                                 |                                                                                                                                                 
tCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGACAGAAATGCCTCCGGGATGCCGCTTTAACTc                                                                                                                                               >  2:349735/1‑149 (MQ=255)
   gTTGCACGCCCATCCGCCGCGCGCAAAACAAAAATGTGTCTGGCGCGGGTTTATGGTGGTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTcccc                                                                                                                                            <  1:231650/149‑1 (MQ=255)
   gTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTcccc                                                                                                                                            <  1:629451/149‑1 (MQ=255)
          gCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTa                                                                                                                                     <  1:226341/149‑1 (MQ=255)
            ccATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACa                                                                                                                                   <  1:378128/149‑1 (MQ=255)
            ccATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTGTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAACATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCACCGATGGCGCTTTAACTCGCCGTTCCTACa                                                                                                                                   >  2:445057/1‑149 (MQ=255)
                         gcACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCAAGGTGTGCCTGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGcgac                                                                                                                      >  1:647786/1‑149 (MQ=255)
                         gcACAACAAAAATGTGTCTGGCGCGGGGTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGcgac                                                                                                                      <  1:630445/149‑1 (MQ=255)
                           acaacaAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGcgacga                                                                                                                    <  1:559929/149‑1 (MQ=255)
                             aacaaAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGAcc                                                                                                                  >  1:443594/1‑149 (MQ=255)
                                  aaaTGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGGTCGGCAGCGACGACGGCGGCAAAATAATGGCGTAATCCCAGGTGCGCCAGTAATGCCTCAGCGAGGGCGCTTTAACTCCCCGTTCCTACAGCCATTGTGGGACTCGCATGcc                                                                                                             >  2:221044/1‑149 (MQ=255)
                                                    tttATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATcaa                                                                                           <  1:535108/149‑1 (MQ=255)
                                                                gcggTGATCGGCAGGGAGCGAGGCGTCAAAAGAAGGGCTTAATCCCAGGTGAGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCg                                                                               <  2:87229/146‑1 (MQ=255)
                                                                 aCGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCgg                                                                              <  2:149938/149‑1 (MQ=255)
                                                                 aCGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCCTCATCCCGGTTTCACCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACATCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCgg                                                                              <  1:122497/149‑1 (MQ=255)
                                                                   gTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAAGCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACCCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGtt                                                                            <  2:461863/149‑1 (MQ=255)
                                                                    ttaTCAGCAGTGGCGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCTCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACTTCAACAAGAGGAAGCGGTTc                                                                           <  2:247958/147‑1 (MQ=255)
                                                                      aTCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGa                                                                         <  1:265117/149‑1 (MQ=255)
                                                                           cAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTa                                                                    <  1:494695/149‑1 (MQ=255)
                                                                                   tacgGCGTCAAAATAATTGTGTTATTCTAGGTGTGTCAGTATTGCTTCAGCGATGGCTTTTTATCTTTCCGTTCTTTCAGTCATTGGGCGACGACCATTTCAACTTTTCACCTCATTAATAATAGTAAGCGGTTCGACGCTATCcagca                                                             <  2:265461/148‑1 (MQ=255)
                                                                                   gacgGCGTCAAAATAATGGCGTACTCCCAGGTGCGCCGGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCCTGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCcagca                                                             <  2:143490/149‑1 (MQ=255)
                                                                                      ggCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCCCCACACCAACCACAGCAAGCGATTCGACCCTATCCAgcacca                                                          >  2:306838/1‑149 (MQ=255)
                                                                                       gCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCAt                                                         >  1:599998/1‑149 (MQ=255)
                                                                                              aataatGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATActtctt                                                  <  1:188225/149‑1 (MQ=255)
                                                                                               ataatGGCGCAATCCCAGGTGCGCCAGCAATGCCTCAGGGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTa                                                 <  2:576262/149‑1 (MQ=255)
                                                                                                taatGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTAc                                                >  1:207850/1‑149 (MQ=255)
                                                                                                 aatGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTCTCCAGCAGCATACTTCTTACt                                               >  1:396275/1‑149 (MQ=255)
                                                                                                     gCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGGCATTGCGCGGCGGACATGCCAACTTTTCACCACATCAGCAAGACGAAGCGTTTCGAGACTATCCACCACCAGGCTTCTTACTGCtt                                           >  2:452030/1‑149 (MQ=255)
                                                                                                             ccAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGtttat                                   <  2:611047/149‑3 (MQ=255)
                                                                                                                    cgcCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCAcgcg                            >  2:1310/1‑149 (MQ=255)
                                                                                                                                 tCAGCGGTGGCGCTTGAATCCCCCGTTCTTAGAGCCATTGGGCTACGACCATGCCAACTTTTCACCACATCGACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTGACTGCTTCTGTTTTTTCACGCGCTAACGCATGCgg               <  2:457228/149‑1 (MQ=255)
                                                                                                                                 tCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTACCGCATGCgg               >  2:281729/1‑149 (MQ=255)
                                                                                                                                              tttAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGcc  <  1:566085/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCC  >  NC_000913_3_hsa_tpiA/2817026‑2817315

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