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% 178.5 / NA 55G263E (GGG→GAG) ydiBquinate/shikimate 5‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/26);  total (29/26)

GTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGATCTCGTCATTGGT  >  NC_000913_3_hsa_tpiA/1774437‑1774721
                                                                                                                                            |                                                                                                                                                 
gTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACa                                                                                                                                          <  1:208618/149‑1 (MQ=255)
   aaTGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGtt                                                                                                                                       >  1:164059/1‑149 (MQ=255)
     tGATATCAGTCTGTTACATCCGGGACTTCTGGTCCCTGACTGCGTGTTTTCATCCTCACAGACTAAAGTATTGCAGCAGGCGCAACAAGCTGGGTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTca                                                                                                                                     <  2:688431/149‑1 (MQ=255)
           cAGTCCGTTCCAGTCGGGACTGCTGTTCACTGAATGCGTGTATACCCTTCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGCTTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTAt                                                                                                                               <  2:418137/149‑1 (MQ=255)
           cAGCCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTAt                                                                                                                               <  1:509972/149‑1 (MQ=255)
            aGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATg                                                                                                                              <  1:652250/149‑1 (MQ=255)
            aGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATg                                                                                                                              >  1:178659/1‑149 (MQ=255)
                        ccGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGa                                                                                                                  >  1:482060/1‑149 (MQ=255)
                         cGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGAt                                                                                                                 >  1:237903/1‑149 (MQ=255)
                               ttCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTca                                                                                                                                     <  1:14990/123‑1 (MQ=255)
                               ttCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTca                                                                                                                                     >  2:14990/1‑123 (MQ=255)
                                tCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCctc                                                                                                          <  2:178659/149‑1 (MQ=255)
                                 cgtgtaaCTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCctct                                                                                                         <  2:845969/143‑1 (MQ=255)
                                 cTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCctct                                                                                                         <  2:164059/149‑1 (MQ=255)
                                   ggTGACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTgg                                                                                                       <  2:705588/149‑1 (MQ=255)
                                     tCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGaa                                                                                                     >  2:303491/1‑149 (MQ=255)
                                        cTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCCGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTTTGGCAAGAGGTTGTACAGTTCACATTATGGACTGGCAAAGATTTCCCTATTGAatat                                                                                                  >  2:524316/1‑149 (MQ=255)
                                         tGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATg                                                                                                                              >  2:160413/1‑120 (MQ=255)
                                         tGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCCGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATg                                                                                                                              <  1:160413/120‑1 (MQ=255)
                                               cGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAAc                                                                                           <  1:292685/149‑1 (MQ=255)
                                                gtgtATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGCACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACa                                                                                          >  1:752728/1‑149 (MQ=255)
                                                 tgtATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAg                                                                                         >  2:135499/1‑149 (MQ=255)
                                                  gtATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGCCTGGCAAAGATTTCCCTCTGGAATATGTTAAACAgg                                                                                        >  1:335638/1‑149 (MQ=255)
                                                  gtATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAgg                                                                                        >  1:431029/1‑149 (MQ=255)
                                                         cccGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATgg                                                                                  >  1:532466/1‑149 (MQ=255)
                                                          ccGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATggg                                                                                 >  2:183454/1‑149 (MQ=255)
                                                                 tGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGt                                                                          <  2:237903/149‑1 (MQ=255)
                                                                  gACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGCTAACATTCTGGACTGGCAAAGCTTTCCCTCTGGAACATGTTAACCAGGACATGGGGTTCGGTg                                                                         >  1:682007/1‑149 (MQ=255)
                                                                   aCGGAGTCATTGCAGCAGTCGCATCAACCGTGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTTTGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCTGTGc                                                                        <  1:336321/149‑1 (MQ=255)
                                                                     gAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGTGTTCGGTGCCt                                                                      <  1:608789/149‑1 (MQ=255)
                                                                     gAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGAGTTCCCTCTGGAATATGTTAAACAGGTAGTGGGGGTCCGTGCCt                                                                      >  2:849114/1‑149 (MQ=255)
                                                                         ttATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTcaca                                                                                                                                   >  2:419108/1‑84 (MQ=255)
                                                                         ttATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTcaca                                                                                                                                   <  1:419108/84‑1 (MQ=255)
                                                                               cagcagGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGa                                                            >  1:581530/1‑149 (MQ=255)
                                                                                     gcgcAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGt                                                      >  1:268727/1‑149 (MQ=255)
                                                                                     gcgcAACAAGCTGGTTGCAAAACGATTGATGGATAAGGCATGTCGTTTTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGt                                                      <  1:511368/149‑1 (MQ=255)
                                                                                       gcAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGCCCGCGTGc                                                    >  1:645148/1‑149 (MQ=255)
                                                                                            aaGGTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGaaa                                               <  2:955125/149‑1 (MQ=255)
                                                                                            aaGCTGGGTGGAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGaaa                                               <  2:547068/149‑1 (MQ=255)
                                                                                             aGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGCCAGGCTGACCGCGTGAAGAAAg                                              >  2:612734/1‑149 (MQ=255)
                                                                                                       aaaaCGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGAttt                                                                                                               >  1:563543/1‑74 (MQ=255)
                                                                                                       aaaaCGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGAttt                                                                                                               <  2:563543/74‑1 (MQ=255)
                                                                                                       aaaaCGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAAt                                    <  2:236581/149‑1 (MQ=255)
                                                                                                       aaaaCGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGAAGAAAGGGTAAAAAAt                                    >  1:752894/1‑149 (MQ=255)
                                                                                                            gATTGATGGATACGGCATGTTGTTGTGGCGAGAGGCTGAACAGTTCACGTTATGGACTGGCCAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGCGCAGAAAGGGTAAAAAATGaaaa                               >  1:812702/1‑149 (MQ=255)
                                                                                                            gATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGaaaa                               >  1:137241/1‑149 (MQ=255)
                                                                                                                 aTGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTa                          >  1:361142/1‑149 (MQ=255)
                                                                                                                      tACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGAt                                                                                                                 >  1:734380/1‑57 (MQ=255)
                                                                                                                      tACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGAt                                                                                                                 <  2:734380/57‑1 (MQ=255)
                                                                                                                       aCGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTa                    <  1:389202/149‑1 (MQ=255)
                                                                                                                          gCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTaaaa                 <  1:756920/149‑1 (MQ=255)
                                                                                                                          gCAGGTTGTTGGGGCAAGAGGCTGAACAGTTCACATTATTGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCCGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTaaaa                 <  1:40871/149‑1 (MQ=255)
                                                                                                                              gttgttgtGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAAATGTAAAAGAtc             >  1:607794/1‑149 (MQ=255)
                                                                                                                                 gctgtgGCAAGAGGGTGAACAGTTCATATTATGGACTGGCAAAGATTCCCCTCTGGAATATGTTAAACAGGTCAGGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGATCTCg          <  2:752894/147‑1 (MQ=255)
                                                                                                                                         aaGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGATCTCGTAATTgct  >  2:284017/1‑147 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
GTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGATCTCGTCATTGGT  >  NC_000913_3_hsa_tpiA/1774437‑1774721

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: