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,5760GA97.7% 132.5 / NA 44G263E (GGG→GAG) ydiBquinate/shikimate 5‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (16/27);  total (16/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 = 1.00e+00

TCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGATCTCGTCATTGG  >  NC_000913_3_hsa_tpiA/1774431‑1774720
                                                                                                                                                 |                                                                                                                                                
tCATTGGTTAATGATATGAGTTTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGCAGTGGCATGTTGTTGTGGCAAGAGGc                                                                                                                                               <  2:56729/149‑1 (MQ=255)
  tttGGTTAATGATATAAGGCTGTTACATCCGGGAATTCTGGTCACTGAATGCGTGTCTAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGAAAAACGATTAATGGTTACGGCATGTTGTTGTGGCAAGAGGCTg                                                                                                                                             <  2:97112/148‑1 (MQ=255)
       ttAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAg                                                                                                                                        >  2:739433/1‑149 (MQ=255)
            gATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTcac                                                                                                                                   <  1:176273/149‑1 (MQ=255)
                       ggtACACCCGGGACTTCTGGTCACCGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTg                                                                                                                        <  1:381991/147‑1 (MQ=255)
                             tCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAg                                                                                                                  <  2:712174/149‑1 (MQ=255)
                                gggACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGAtt                                                                                                               >  1:772728/1‑149 (MQ=255)
                                gggACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGCCAAAGAtt                                                                                                               >  1:173852/1‑149 (MQ=255)
                                gggACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGAAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGAtt                                                                                                               >  2:503224/1‑149 (MQ=255)
                                                    gCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTaaa                                                                                           <  1:583032/149‑1 (MQ=255)
                                                      gtgtATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCATCTGGAATATGTTAAACa                                                                                         >  2:607107/1‑149 (MQ=255)
                                                      gtgtATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGAAAAACGCTGATTGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACa                                                                                         <  1:404893/149‑1 (MQ=255)
                                                         tataACCGGCATATGACGAAGTTATGGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTGGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGt                                                                                      <  2:591130/149‑1 (MQ=255)
                                                          ataACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTc                                                                                     <  1:66149/149‑1 (MQ=255)
                                                            aaCCCGCATATGACGACGTTACTGCAGCAGTCGCAACCAGCTCGCCCCCAAACCATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCAt                                                                                   <  1:401202/149‑1 (MQ=255)
                                                              cccGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATgg                                                                                 >  1:33740/1‑149 (MQ=255)
                                                               ccGCATATGACGAAGTTATTGCAGCAGGAGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGGCCTggg                                                                                >  1:393742/1‑149 (MQ=255)
                                                                    tatGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCg                                                                           >  2:732859/1‑149 (MQ=255)
                                                                      tGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGt                                                                         >  2:627371/1‑149 (MQ=255)
                                                                        aCGAAGTTGTTGGAGCAGGCGCAACAAGCTGGGTGCAAAACGATTGATGGGTACGGCATGTTGTTGTGGCAAGGGGCTGAACAGTTCACATTTTGGAGTGGGAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGc                                                                       <  2:555511/149‑1 (MQ=255)
                                                                         cGAAGTTATTGCAGCAGGCGTAACGAGTTGGTGGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGcc                                                                      <  2:173852/149‑1 (MQ=255)
                                                                         cGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGcc                                                                      >  2:464854/1‑149 (MQ=255)
                                                                          gAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCt                                                                     <  1:781624/149‑1 (MQ=255)
                                                                          gAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCt                                                                     <  1:349693/149‑1 (MQ=255)
                                                                            aGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGa                                                                   <  1:143663/149‑1 (MQ=255)
                                                                                          gcgcAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGt                                                     >  1:164557/1‑149 (MQ=255)
                                                                                                  aGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAg                                             <  1:732859/149‑1 (MQ=255)
                                                                                                   gCTTGTTGGAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAgg                                            <  2:772728/149‑1 (MQ=255)
                                                                                                   gCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGa                                                           >  2:663539/1‑134 (MQ=255)
                                                                                                   gCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGa                                                           <  1:663539/134‑1 (MQ=255)
                                                                                                   gCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAgg                                            <  1:44153/149‑1 (MQ=255)
                                                                                                            aaaaCGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACGTTATGTACAGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTGATGGGGCTCAGTGCCTGACAGGCTGACCGCGTGCAGAAAGAGTaaaaaag                                   >  2:324659/1‑148 (MQ=255)
                                                                                                            aaaaCGATTGATGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAAt                                   >  2:237772/1‑149 (MQ=255)
                                                                                                                       tGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTaa                        <  2:109582/149‑1 (MQ=255)
                                                                                                                       tGGATACGGCATGTTGTTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTaa                        <  1:820151/149‑1 (MQ=255)
                                                                                                                       tGGATACGGCATGTTAGTGTGGCAAGAGGCTGAACAGTTCACATTATGGACTTGTAAAGATTTCCCGCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTaa                        <  2:63808/149‑1 (MQ=255)
                                                                                                                              ggCATGTTGTTGTGGCGAGAGGAGGAAGAGTTCAAATTATGGACTGGCAAAGAGTTGCATCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGGAAAAAATGAAAACCGTAACTGTaaa                 <  2:393742/149‑1 (MQ=255)
                                                                                                                                  tgttgttgtGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGAt             >  2:488259/1‑149 (MQ=255)
                                                                                                                                     tgttgtGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGAtctc          <  1:754350/149‑1 (MQ=255)
                                                                                                                                     tgttgtGGCAAGAGGCTGAACAGTTAACATTGTGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGAtctc          <  1:347231/149‑1 (MQ=255)
                                                                                                                                     tgttgtGGCAAGAGGCTGAACAGGTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGAtctc          <  1:428586/149‑1 (MQ=255)
                                                                                                                                      gttgtGGCAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGATCTCg         <  1:267588/149‑1 (MQ=255)
                                                                                                                                            gCAAGAAGCTGAACAGTTCACATTATGGACTGGCAAAGAGTTGCCTCAGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGAAAAAAAAAAAAAAAATGAAAACCGTAACTGTAAAAGATCTCGTCATTg   <  2:164557/149‑1 (MQ=255)
                                                                                                                                             cAAGAGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGGTCGGTGCCTGACAGGCTGAACGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGATCTCGTCATTgg  >  2:760636/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGGGGCTGAACAGTTCACATTATGGACTGGCAAAGATTTCCCTCTGGAATATGTTAAACAGGTCATGGGGTTCGGTGCCTGACAGGCTGACCGCGTGCAGAAAGGGTAAAAAATGAAAACCGTAACTGTAAAAGATCTCGTCATTGG  >  NC_000913_3_hsa_tpiA/1774431‑1774720

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: