Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	2,418,678	G→A	7.8%	R167C (CGC→TGC)	yfcD ←	putative NUDIX hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	2,418,678	0	G	A	7.8%	126.8 / 2.5	52	R167C (CGC→TGC)	yfcD	putative NUDIX hydrolase
Reads supporting (aligned to +/- strand):  ref base G (29/18);  new base A (2/2);  total (32/20)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.39e-01

CCTACTTGGACTGCTGTTTATTATGTCCAGCGTCGTAGTGATCGGCGCTCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCG  >  NC_000913_3_hsa_tpiA/2418532‑2418826                                                                                                                                                   |                                                                                                                                                     ccTACTTGGACTGCTGTTTATTATGTCCAGCGTCGTAGTGATCGGCGCTCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCACt                                                                                                                                                    >  1:558264/1‑149 (MQ=255) ccTACTTGGACTGCTGTTTATTATGTCCAGCGTCGTAGTGATCGGCGCTCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCACt                                                                                                                                                    >  1:558261/1‑149 (MQ=255)   tagttGGACTGCCGTTTATTATGTCCAGCGTCGTAGTGATCGGCGCTCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTACGACTTTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCCTTCTTGCAATTGCGCTTc                                                                                                                                                  >  2:116882/4‑149 (MQ=255)      ttGGACTGCTGTTTATTATGTCCAGCGTCGTAGTGATCGGCGCTCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATc                                                                                                                                               <  1:669969/149‑1 (MQ=255)      ttGGACTGCTGTTTATTATGTCCAGCGTCGTAGTGATCGGCGCTCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCCCGACTGCTTCATTTTTGGCATTGCGCTTCATc                                                                                                                                               <  1:833970/149‑1 (MQ=255)              ctgTTTATTATGTCCAGCGTCGTAGTGATCGGCGCTCAGTTGATGGGCTACCGCTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAActc                                                                                                                                       >  2:50275/1‑147 (MQ=255)              ctgTTTATTATGTCCAGCGTCGTAGTGATCGGCGCTCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTCCGGCTTCATTTTTGGAGTTGCGCTTCATCCACAAcac                                                                                                                                       >  1:384889/1‑147 (MQ=255)                gTTGATTACGGCCGGCGTGCTGGTGATGGGCGCCCGGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCCACTGCGGCCCTTTTTTTATTCTGCCGTCTCAGTCGCTGCGGCTTCATTTTTGGCATTGCGCTTCATCCACCACGCTa                                                                                                                                     <  2:48147/149‑1 (MQ=255)                 tttATTATGTCCAGCGTCGTAGTGATCGGCGCTCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAg                                                                                                                                    >  1:820811/1‑149 (MQ=255)                      tatGTCCAGCGTCGTAGTGATCGGCGCTCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACCCTAACGCtt                                                                                                                               >  2:352400/1‑149 (MQ=255)                           ccAGCGGCGTAGTGATCGGCGCTCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGc                                                                                                                          <  1:105401/149‑1 (MQ=255)                                   gtAGTGATCGGCGCTCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCt                                                                                                                                >  2:733628/1‑135 (MQ=255)                                   gtAGTGATCGGCGCTCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCt                                                                                                                                <  1:733628/135‑1 (MQ=255)                                    tAGTGATCGGCGCTCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGa                                                                                                                 <  1:752470/149‑1 (MQ=255)                                         aTCGGCTCGTAGTTGATGGGCTACCACTAAAAATGTTAAGGGCCGCAATTTCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGaa                                                                                                            <  2:653896/149‑1 (MQ=255)                                         aTCGGCGCTCAGTTGCTGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCCTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGaa                                                                                                            <  1:237291/149‑1 (MQ=255)                                         aTCGGCGCTCAGTGGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATACTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGATTCATCCGCAACGCTAGCGCTTTCAGCGAGTCCGGAGTgca                                                                                                            >  2:67125/1‑147 (MQ=255)                                              cgcTCAGTTGATGGGCTACCACTAACAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCCGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCAt                                                                                                       <  1:324081/149‑1 (MQ=255)                                                cTCAGTTGATGGCCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTTGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTTGGGTGAACTCATCg                                                                                                     >  2:345057/1‑149 (MQ=255)                                                   aGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCAAATAGCACACGTCTGTTCATCCACAACGCTAGCGCTTTCACCGAGTCTGGAGTGACCTCATCGCAg                                                                                                  >  1:400930/1‑149 (MQ=255)                                                      tGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGt                                                                                               <  1:41579/149‑1 (MQ=255)                                                                  aCTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGGTttct                                                                                   >  1:532204/1‑149 (MQ=255)                                                                      aaaaTGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCg                                                                               >  1:444155/1‑149 (MQ=255)                                                                      aaaaTGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGACTTCTTCCg                                                                               >  2:269722/1‑149 (MQ=255)                                                                        aagtttAAGAGCCGCAATTGCGGCTCCTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGc                                                                             <  1:94285/145‑1 (MQ=255)                                                                        aaTGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTTCGCTTCATCCACAACGCCAGCGCTTTCAGCGCGTCTGGAGTGAACTCATCGCAGCGTGCGCTTATTTCTTCCGGc                                                                             >  2:394662/1‑149 (MQ=255)                                                                              aaGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTcggc                                                                       >  1:736441/1‑149 (MQ=255)                                                                                gagCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTATGGAGTGAACTCATCGCAGCGTGCGGTTATTTCTTCCGGCGTcaccca                                                                     >  1:493519/1‑149 (MQ=255)                                                                                   ccGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGCGTCTGGAGTGAACTCATCGCAGCGTGCGGTGCTTTCTTACAGCGTcacgctgca                                                                  >  2:660624/1‑142 (MQ=255)                                                                                   ccGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGCGTCTACAGTGAACTCATCGCAGCGTGCGGTGCTTTCTTCCGGCGTCAGCCAGCa                                                                  >  1:112417/1‑149 (MQ=255)                                                                                    cGCAATTGCGGCTCTTTTTTTTTTTTTTTTTTTCAGTCTAAACGGCTTCATTTTTGGCATTGGGCTTCATCCACAACGCTAGCGCTTTCACCGACTCTGGACTGACCTCATCGCACTTTGTGGTTATTTCTTCCGCCGTCACCCCGCaa                                                                 >  2:451992/1‑149 (MQ=255)                                                                                     gctaTTGCGGCTTTTTTTTCATTCTGCCGTTTTAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCaaa                                                                <  2:530128/146‑1 (MQ=255)                                                                                     gCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAATCTGGAATGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCaac                                                                >  2:426253/1‑148 (MQ=255)                                                                                             ggCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAACCTTCACTg                                                        >  1:601194/1‑149 (MQ=255)                                                                                              gCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACACCGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTACGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGa                                                       >  1:123323/1‑149 (MQ=255)                                                                                                        cATTCTGCCGTTTCAGTCTTTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCTTGCGGTCATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTc                                             <  2:635083/149‑1 (MQ=255)                                                                                                         aTTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTcc                                            >  1:390354/1‑149 (MQ=255)                                                                                                                       gTCTCTACGGCTTCATTTTTGGCATTGCACTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAggg                              <  2:558264/149‑1 (MQ=255)                                                                                                                       gTCTCTACGGCTTCATTTTTGGCATTGCACTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAggg                              <  2:558261/149‑1 (MQ=255)                                                                                                                          tctACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTTCGCAGGGAcc                           >  1:8909/1‑149 (MQ=255)                                                                                                                          tctACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGCCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGAcc                           >  2:147612/1‑149 (MQ=255)                                                                                                                                 cTTGATTTTTGGCAGTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGGGTGAACTCATCGCAGCGTGCGGTGATTGCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTgagag                    <  2:255845/149‑1 (MQ=255)                                                                                                                                 cTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGGGTCTGGAGTGAACTCATCGCAGCGTGCGGTGGTTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGATGGGACCGTgagag                    >  2:678057/1‑149 (MQ=255)                                                                                                                                 cTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGACCTCATCGCAGCGTGCGGTGCTTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCGTCTTCCTGTAGTGCGAACGGCCCGagagag                    >  2:253338/1‑149 (MQ=255)                                                                                                                                  ttCATTGTTTGCCTCGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTTTCGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTgagaga                   <  2:272513/149‑1 (MQ=255)                                                                                                                                      tttttGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCGTCTTCCGGTAGTGCGAAGGGGCCGTGAGAGACGCa               >  1:541701/1‑149 (MQ=255)                                                                                                                                      tttttGGCATTGCGCTGCATCCACAACGCTAGCGCTTTCAGCGAGTCTTGAGTTAACTCCTCGCGGCGCGCGGTGGTTTCTTCCGCTGCCTTGCCCTATACTGCTCTGACGTCACCTTCCTGTAGTTCCAAGGGACCGTAGGAGAAGCa               >  2:174663/1‑149 (MQ=255)                                                                                                                                         ttGGCATTGCGTTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCCTCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCt            <  2:684170/149‑1 (MQ=255)                                                                                                                                          tGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTg           <  2:788537/149‑1 (MQ=255)                                                                                                                                            gCATGGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTCGTCCGAAGGTCCCCTGAGCGACGCAGCTGaa         >  2:639752/1‑149 (MQ=255)                                                                                                                                              aTTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAGGGGACCGTGAGAGACGCAGCTGAACa       >  1:352005/1‑149 (MQ=255)                                                                                                                                                ttcgcTTCATCAACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAAt     <  2:411984/147‑1 (MQ=255)                                                                                                                                                tGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATGGCAGCGTGCGGTGATTTCTTCTGGCGTCAGCCAGCAAACTTCACTGACTTCATCATCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAAt     <  2:798443/149‑1 (MQ=255)                                                                                                                                                 gcgcTTCATCCACAACGCTAGCGCTTTCAGCTAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATg    <  1:288182/149‑1 (MQ=255)                                                                                                                                                   gcTTCATCCACAACGCTAGGGCTTTCGGCGAGTCTGGAGTGAACGCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCTAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGGTGAACAATgcg  <  2:313514/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CCTACTTGGACTGCTGTTTATTATGTCCAGCGTCGTAGTGATCGGCGCTCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCG  >  NC_000913_3_hsa_tpiA/2418532‑2418826

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: ‑