Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 3,560,455 +G intergenic (‑2/+1) glpR ← / ← glpR pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA3,560,4551.G97.6% 132.4 / ‑3.9 41intergenic (‑2/+1)glpR/glpRpseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (15/25);  total (16/25)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTT  >  NC_000913_3_hsa_tpiA/3560311‑3560597
                                                                                                                                                 |                                                                                                                                               
atatCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGcgcag                                                                                                                                              >  2:800365/1‑145 (MQ=255)
     gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGCACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCa                                                                                                                                        >  2:245560/1‑149 (MQ=255)
       tAAACAGCGTCGAGCCATTGGGGATTTGCTACGACACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGTCCTTGCGAAAGTGCCACGGCGTGTTAACCAGACTGGCAGGCAGCGCCGGCCCGCACAGCTGGC‑CACGGATCAgg                                                                                                                                      >  2:736847/1‑149 (MQ=255)
          aCAGCGTCGGGCCGTTGGGGATTTGCCCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCCGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGttt                                                                                                                                   <  2:96923/149‑1 (MQ=255)
          aCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGttt                                                                                                                                   <  1:179286/149‑1 (MQ=255)
          aCAGCGTCGAGCCATTGGCGATTTGCTCCGCCGCGTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGttt                                                                                                                                   <  2:527440/149‑1 (MQ=255)
          aCAGCGCCGGGCCATTGTGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGttt                                                                                                                                   <  2:665323/149‑1 (MQ=255)
            aGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGGTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTg                                                                                                                                 <  2:108625/149‑1 (MQ=255)
                tccaGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTcc                                                                                                                             <  1:54735/146‑1 (MQ=255)
                tCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTcc                                                                                                                             >  1:54597/1‑149 (MQ=255)
                      cATTGGGGATTTGCTCCGCCACTTTTCGGGCGATGCGCTCTTTTTCTTCGTTCTGGGTGGCCCTGCGATCGTTCCACGGCGTGTTAACCGCACTGGAAGGCAGCGCCGCACCGCCATGCTCGCGAA‑GGATCACGTCTGGCTCCACCAGc                                                                                                                       >  2:62695/1‑149 (MQ=255)
                                   cTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgc                                                                                                          >  1:624843/1‑149 (MQ=255)
                                   cTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgc                                                                                                          >  1:214866/1‑149 (MQ=255)
                                            tttGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAg                                                                                                 >  1:526029/1‑149 (MQ=255)
                                                 gggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAg                                                                                                                                       <  1:520194/106‑1 (MQ=255)
                                                 gggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAg                                                                                                                                       >  2:520194/1‑106 (MQ=255)
                                                            tctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGga                                                                                 <  2:809391/149‑1 (MQ=255)
                                                               ttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaa                                                                              <  1:600297/149‑1 (MQ=255)
                                                                tttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGcaa                                                                             <  2:640573/149‑1 (MQ=255)
                                                                    ttcGGTTTGGGTGGCCGTGCGATCGGGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                                         <  2:670763/149‑1 (MQ=255)
                                                                       ggTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTa                                                                      <  1:71844/149‑1 (MQ=255)
                                                                          cTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCa                                                                   >  1:698417/1‑149 (MQ=255)
                                                                                 gCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctc                                                                        >  1:640356/1‑137 (MQ=255)
                                                                                 gCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctc                                                                        <  2:640356/137‑1 (MQ=255)
                                                                                 gCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTc                                                            <  1:479816/149‑1 (MQ=255)
                                                                                   cTTGCGATGGGCCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCgg                                                          <  2:130085/149‑1 (MQ=255)
                                                                                   cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCgg                                                          >  2:265836/1‑149 (MQ=255)
                                                                                    ttGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGGTGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTGCGGt                                                         >  1:173695/1‑149 (MQ=255)
                                                                                      gggATCGGGCCAGGGCGTGTTGACCGAACGGGAAGGCCGCGCGGCAACGGGATGATGGCGCA‑GGATCAGGTTTTGCTGCGCTAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAc                                                       <  2:176288/147‑1 (MQ=255)
                                                                                            gTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAt                                                 <  1:265836/149‑1 (MQ=255)
                                                                                               ccACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAc                                              <  1:93401/149‑1 (MQ=255)
                                                                                                     cGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgct                                        >  1:443726/1‑149 (MQ=255)
                                                                                                      gtgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGACCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgctg                                       >  1:33340/1‑149 (MQ=255)
                                                                                                           gtaccaaCTGGATGGCAGCGCCGCACCGCCATCATGGCGCA‑GCAGCGGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGTGGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGTTCTTCGTTACTGACATAACCCTGCTGTTTaa                                  <  2:331696/144‑1 (MQ=255)
                                                                                                                 aCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcccatt                            >  2:424973/1‑145 (MQ=255)
                                                                                                                  cgggAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                           <  1:585294/147‑1 (MQ=255)
                                                                                                                    ggAAGGCAGCGCCGCACCGCCATTATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGa                         <  1:61823/149‑1 (MQ=255)
                                                                                                                     gAAGGCAGCGCCGCACCGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGat                        <  2:526029/149‑1 (MQ=255)
                                                                                                                                   aCCGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGAc          <  1:800365/149‑1 (MQ=255)
                                                                                                                                    ccGCCATGATGGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACg         <  2:443726/149‑1 (MQ=255)
                                                                                                                                        cacgatgGCGCA‑GGATCAGGTTTTGCTCCTCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCATTTCGATAATACCGTTGTGACGTtgt     <  2:173695/146‑1 (MQ=255)
                                                                                                                                           gatgGCGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtt  >  1:627621/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
ATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTT  >  NC_000913_3_hsa_tpiA/3560311‑3560597

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