Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 1,774,576 G→A G263E (GGG→GAG)  ydiB → quinate/shikimate 5‑dehydrogenase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA1,774,5760GA100.0% 148.1 / NA 44G263E (GGG→GAG) ydiBquinate/shikimate 5‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/26);  total (18/26)

TTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGATCTCGTCATTGG  >  NC_000913_3_hsa_tpiA/1774434‑1774720
                                                                                                                                               |                                                                                                                                                
ttGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGa                                                                                                                                            >  2:355203/1‑149 (MQ=255)
   gTTAATGATAACAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACa                                                                                                                                         <  1:414419/149‑1 (MQ=255)
    ttAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGACGTTATTGCAGCAGGCG‑AACAACGAAAATTTCATATCGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAg                                                                                                                                        >  1:142509/1‑149 (MQ=255)
       aTGATATCAGTCCGTTACATCCGGGACGTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTc                                                                                                                                     <  1:178016/149‑1 (MQ=255)
        tGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTca                                                                                                                                    <  1:688562/149‑1 (MQ=255)
                     ttACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTgg                                                                                                                       <  1:194413/149‑1 (MQ=255)
                      tACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACTAAGTTATTGCAACAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGGCTGGc                                                                                                                      >  1:244936/1‑149 (MQ=255)
                      tACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGc                                                                                                                      <  1:580734/149‑1 (MQ=255)
                           cgtggACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGa                                                                                                                 <  2:733178/146‑1 (MQ=255)
                              ggACTTCTTGTCACTGAATGCGTGTATAATCCGCATAGGACGAAGTTATTGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGAttt                                                                                                              <  2:559230/149‑1 (MQ=255)
                              ggACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGCTCACCTTATGGACTAGGAAAGGttt                                                                                                              >  2:378709/1‑149 (MQ=255)
                                aCTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAT‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTcc                                                                                                            <  2:678952/149‑1 (MQ=255)
                                          aCTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTg                                                                                                       >  2:426329/1‑144 (MQ=255)
                                          aCTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTg                                                                                                       <  1:426329/144‑1 (MQ=255)
                                          aCTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAata                                                                                                  <  1:424126/149‑1 (MQ=255)
                                            tGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATg                                                                                                <  2:357987/149‑1 (MQ=255)
                                                 gCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAA‑GCCGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTaaa                                                                                           <  2:412441/149‑1 (MQ=255)
                                                     gtATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAgg                                                                                       >  2:417578/1‑149 (MQ=255)
                                                         aaCCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAACGATTTCCCTCTCGAATATTTTAACCAGATcgt                                                                                   >  2:300640/1‑147 (MQ=255)
                                                           ccggCATAGGACGAAGTTATTGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGCCTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATgg                                                                                 <  2:588381/146‑1 (MQ=255)
                                                                     aCGAAGTTATTGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGc                                                                       >  1:168518/1‑149 (MQ=255)
                                                                       gAAGTTAGTGCAGCAGGCGCCACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCt                                                                     <  2:187219/149‑1 (MQ=255)
                                                                         aGTTATTGGAGTAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGa                                                                   <  2:34951/149‑1 (MQ=255)
                                                                         aGTTATTGCAGCAGGCGCAGCAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGa                                                                   <  2:35903/149‑1 (MQ=255)
                                                                              ttGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTg                                                                                                       <  2:77647/108‑1 (MQ=255)
                                                                              ttGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTg                                                                                                       >  1:77647/1‑108 (MQ=255)
                                                                               tGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGCTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCt                                                             <  1:261993/149‑1 (MQ=255)
                                                                                gcagcaGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTg                                                            >  1:147457/1‑149 (MQ=255)
                                                                                     agGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACcgc                                                       <  1:417578/149‑1 (MQ=255)
                                                                                         gcAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGAGGTTCGGTGCCTGACAGGCTGACCGCGTgg                                                   >  2:180463/1‑148 (MQ=255)
                                                                                           aacaa‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAg                                                 >  2:191592/1‑149 (MQ=255)
                                                                                            acaa‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGa                                                >  2:256656/1‑149 (MQ=255)
                                                                                                   tGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGt                                          <  1:256656/149‑1 (MQ=255)
                                                                                                          aaaaCGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAgg                                                               <  2:375070/121‑1 (MQ=255)
                                                                                                          aaaaCGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAgg                                                               >  1:375070/1‑121 (MQ=255)
                                                                                                          aaaaCGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAAt                                   >  1:270918/1‑149 (MQ=255)
                                                                                                          aaaaCGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAAt                                   <  1:191592/149‑1 (MQ=255)
                                                                                                                 ttGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAAcc                            <  2:147457/149‑1 (MQ=255)
                                                                                                                       gATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAAcg                      >  1:528431/1‑148 (MQ=255)
                                                                                                                             gCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTaaaa                <  1:355203/149‑1 (MQ=255)
                                                                                                                                     ttgtGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATg                                                                                                >  1:605061/1‑61 (MQ=255)
                                                                                                                                     ttgtGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATg                                                                                                <  2:605061/61‑1 (MQ=255)
                                                                                                                                          gCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGATCTCGTCATTg   >  1:70606/1‑149 (MQ=255)
                                                                                                                                           cAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGATCTCGTCATTgg  <  1:242607/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
TTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAA‑GCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGATCTCGTCATTGG  >  NC_000913_3_hsa_tpiA/1774434‑1774720

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