Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	1,748,487	A→T	37.3%	I453F (ATC→TTC)	ydhS →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	1,748,487	0	A	T	37.3%	48.3 / 66.2	51	I453F (ATC→TTC)	ydhS	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (18/14);  new base T (9/10);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 5.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-01

GCGGGAATCATTCATATTCTCGCCCTCGGTGAAGACTACAAAATGGAAATTAATGAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAACCCGAAGATATTCGTGGGCGCGTAGCGTTCGGCGCGTTGCCCTGGTTGATGCACGACCAGCCTTTC  >  NC_000913_3_hsa_tpiA/1748340‑1748636                                                                                                                                                    |                                                                                                                                                      gCGGGAATCATTCATATTCTCGCCCTCGGTGAAGACTACAAAATGGAAATTAATGAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATTTc                                                                                                                                                     <  3:93190/150‑1 (MQ=255)            tcatATTCTCGCCCTCGGTGAAGACTACAAAATGGAAATTAATGAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTgg                                                                                                                                          >  1:500844/1‑150 (MQ=255)             catATTCTCGCCCTCGGTGAAGACTACAAAATTGAAATTAATGAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTggg                                                                                                                                         <  2:382284/150‑1 (MQ=255)             catATTCTCGCCCTCGGTGAAGACTACAAAATGGAAATTAATGAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGTCGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATTTCCCTTTCCCTggg                                                                                                                                         <  1:132736/150‑1 (MQ=255)                atTCTCGCCCTCGGTGAAGACTACAAAATGGAAATTAATGAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGATTAAAGTGAAAGATATCCCTTTCCCTGGGCTa                                                                                                                                      >  2:98494/1‑150 (MQ=255)                      gCCCTCGGTGAAGACTACAAAATGGAAATTAATGAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAg                                                                                                                                <  1:349649/150‑1 (MQ=255)                      gCCCTCGGTGAAGACTACAAAATGGAAATTAATGAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAg                                                                                                                                <  1:42228/150‑1 (MQ=255)                         cTCGGTGAAGACTACAAAATGGAAATTAATGAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATTTCCCTTTCCCTGGGCTACGCGAGCaa                                                                                                                             >  4:325233/1‑150 (MQ=255)                          tCGGTGAAGACTACAAAATGGAAATTAATGAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATTTCCCTTTCCCTGGGCTACGCGAGCAAt                                                                                                                            <  4:92257/150‑1 (MQ=255)                            ggTGAAGACTACAAAATGGAAATTAATGAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTa                                                                                                                          <  2:465191/150‑1 (MQ=255)                              tGAAGACTACAAAATGGAAATTAATGAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATTTCCCTTTCCCTGGGCTACGCGAGCAATTACa                                                                                                                        >  2:83656/1‑150 (MQ=255)                              tGAAGACTACAAAATGGAAATTAATGAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACa                                                                                                                        <  4:13342/150‑1 (MQ=255)                              tGAAGACTACAAAATGGAAATTAATGAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCCACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATTTCCCTTTCCCTGGGCTACGCGAGCAATTACa                                                                                                                        >  2:83569/1‑150 (MQ=255)                                aaGACTACAAAATGGAAATTAATGAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCGTTTCCCTGGGCTACGCGAGCAATTACAGa                                                                                                                      >  2:367689/1‑150 (MQ=255)                                aaGACTACAAAATGGAAATTAATGAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGa                                                                                                                      >  1:277847/1‑150 (MQ=255)                                 aGACTAAAAAATGGAAATTAATGAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGaa                                                                                                                     <  4:223295/150‑1 (MQ=255)                                  gACTACAAAATGGAAATTAATGAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCCTATAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGaaa                                                                                                                    >  3:196534/1‑150 (MQ=255)                                   aCTACAAAATGGAAATTAATGAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGaaaa                                                                                                                   >  3:266735/1‑150 (MQ=255)                                      aCAAAATGGAAATTAATGAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATTTCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAg                                                                                                                <  1:457108/150‑1 (MQ=255)                                                atttAATGAGTCGCGCACCGTCCTGAAAACGGAATACAACAGTTACTTGTGTGACGTTTTTATTGATGCCTGCGGACAGCGTTCGCTTAAAGTGAATGATATCTCTTTCCTTGGGCTACGCGAGCAATTATAGAAAACAGTGTCTGaaaa                                                                                                      >  4:287922/3‑149 (MQ=255)                                                       gAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATTTCCCTTTCCCTGTGCTACGCGAGGAATTACAGAAAACAGGGGATGAAATCCCTGac                                                                                               >  4:359725/1‑149 (MQ=255)                                                       gAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATTTCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGAt                                                                                               >  3:106202/1‑150 (MQ=255)                                                               cACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATTTCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGa                                                                                       <  2:304913/150‑1 (MQ=255)                                                                aCCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATTTCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGaa                                                                                      <  3:216146/150‑1 (MQ=255)                                                                  cGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATTTCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa                                                                                    <  2:207245/150‑1 (MQ=255)                                                                  cGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATTTCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGa                                                                                    >  1:433373/1‑150 (MQ=255)                                                                       tGAAAAGGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGGGGACAACGTCCGCTTAAAGTGAAAGATTTCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATtata                                                                               <  2:99169/150‑1 (MQ=255)                                                                       tGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATTTCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATtata                                                                               <  2:132770/150‑1 (MQ=255)                                                                       tGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATTTCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATtata                                                                               <  4:373810/150‑1 (MQ=255)                                                                           aaCGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGtt                                                                           >  2:576263/1‑150 (MQ=255)                                                                            aCGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTa                                                                          <  4:273056/150‑1 (MQ=255)                                                                                  gACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCaa                                                                    <  3:297209/150‑1 (MQ=255)                                                                                    caacaGCTACTCGTTTGACGTTTTTATTCATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAAcc                                                                  <  2:79621/150‑1 (MQ=255)                                                                                     aacaGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAAccc                                                                 >  4:202355/1‑150 (MQ=255)                                                                                      acaGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAACCCg                                                                <  1:299243/150‑1 (MQ=255)                                                                                         gctgctCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCGTAAAGTGAAAGATATCCCGTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAACCCGAAg                                                             <  2:420397/150‑1 (MQ=255)                                                                                         gCTACTCGGTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAACCCGAAg                                                             >  4:82284/1‑150 (MQ=255)                                                                                            aCTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATTTCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAACCCGAAGata                                                          >  2:144855/1‑150 (MQ=255)                                                                                              tCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAACCCGAAGATAtt                                                        >  3:69705/1‑150 (MQ=255)                                                                                                    gACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAACCCGAAGATATTCGTggg                                                  <  1:164573/150‑1 (MQ=255)                                                                                                     aCGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAACCCGAAGATATTCGTGGgc                                                 >  1:552410/1‑150 (MQ=255)                                                                                                        tttttATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATTTCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAACCCGAAGATATTCGTGGGcgcg                                              >  2:32417/1‑150 (MQ=255)                                                                                                               tGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATTTCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAACCCGAAGATATTCGTGGGCGCGTAGCGtt                                       >  2:17226/1‑150 (MQ=255)                                                                                                               tGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAACCCGAAGATATTCGTGGGCGCGTAGCGtt                                       >  2:44165/1‑150 (MQ=255)                                                                                                                          gACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAACCCGAAGATATTCGTGGGCGCGTAGCGTTCGGCGCGTTGc                            >  2:299205/1‑150 (MQ=255)                                                                                                                            cAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAACCCGAAGATATTCGTGGGCGCGTAGCGTTCGGCGCGTTGccc                          >  3:123128/1‑150 (MQ=255)                                                                                                                                 tCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAACCCGAAGATATTCGTGGGCGCGTAGCGTTCGGCGCGTTGCCCTGGtt                     >  3:317239/1‑150 (MQ=255)                                                                                                                                      ttAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAACCCGAAGATATTCGTGGGCGCGTAGCGTTCGGCGCGTTGCCCTGGTTGATGc                >  4:26714/1‑150 (MQ=255)                                                                                                                                            tGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAACCCGAAGATATTCGTGGGCGCGTAGCGTTCGGCGCGTTGGCCTGGTTGATGCACGAcc          >  4:160006/1‑150 (MQ=255)                                                                                                                                               aaGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAACCCGAAGATATTCGTGGGAGACTAGCGTTCGGCGCGTTGCCCTGGTTGATGCACGACCAGc       <  4:387922/150‑1 (MQ=255)                                                                                                                                                 gATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAACCCGAAGATATTCGTGGGCGCGTAGCGTTCGGCGCGTTGCCCTGGTTGATGCACGACCAGCCt     <  3:59173/150‑1 (MQ=255)                                                                                                                                                    atCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAACCCGAAGATATTCGTGGGCGCGTAGCGTTCGGCGCGTTGCCCTGGTTGATGCACGACCAGCCTTTc  <  2:525613/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GCGGGAATCATTCATATTCTCGCCCTCGGTGAAGACTACAAAATGGAAATTAATGAGTCGCGCACCGTCCTGAAAACGGAAGACAACAGCTACTCGTTTGACGTTTTTATTGATGCCCGCGGACAACGTCCGCTTAAAGTGAAAGATATCCCTTTCCCTGGGCTACGCGAGCAATTACAGAAAACAGGGGATGAAATCCCTGATGTTGGCGAAGATTATACGTTACAGCAACCCGAAGATATTCGTGGGCGCGTAGCGTTCGGCGCGTTGCCCTGGTTGATGCACGACCAGCCTTTC  >  NC_000913_3_hsa_tpiA/1748340‑1748636

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