Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	1,544,021:1	+A	76.5%	coding (40/1389 nt)	narU ←	nitrate/nitrite transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	1,544,021	1	.	A	76.5%	145.9 / 43.2	51	coding (40/1389 nt)	narU	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (8/4);  new base A (19/20);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTG  >  NC_000913_3_hsa_tpiA/1543872‑1544163                                                                                                                                                       |                                                                                                                                                 gATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTcgcg                                                                                                                                                   >  1:291194/1‑150 (MQ=255)   ttttATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTcgcgaa                                                                                                                                                 <  1:175151/150‑3 (MQ=255)        ttGAGATTAACGTTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑a‑aaaag                                                                                                                                         <  2:578585/150‑3 (MQ=255)            gaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACTACTGCCTGATATACAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATaa                                                                                                                                      <  1:201721/150‑1 (MQ=255)              ttAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACg                                                                                                                                    >  2:217752/1‑150 (MQ=255)                       aCTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTAttttt                                                                                                                           <  3:276119/150‑1 (MQ=255)                                aaTAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTg                                                                                                                  >  1:518276/1‑150 (MQ=255)                                aaTAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTg                                                                                                                  <  4:285667/150‑1 (MQ=255)                                  tAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                <  1:398692/150‑1 (MQ=255)                                  tAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAA‑AAAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                >  2:109775/1‑150 (MQ=255)                                     gaTCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTg                                                                                                             <  2:387453/149‑1 (MQ=255)                                     cATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTg                                                                                                             <  1:468907/150‑1 (MQ=255)                                         cAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCat                                                                                                         >  1:221170/1‑150 (MQ=255)                                          aGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGCTTCCAGTCGCG‑A‑ACAAAGATAACGACTATTTTTCTCATTTTGCAGTTCCaga                                                                                                        >  2:318388/1‑148 (MQ=255)                                             cacaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTg                                                                                                     <  4:7182/150‑1 (MQ=255)                                             aacaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTTTGGTTTCAAAAAAAAAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTg                                                                                                     <  2:393698/149‑1 (MQ=255)                                                     gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGTTTTTCTGGTTTCCAGTTGCG‑AAA‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcac                                                                                             >  1:285022/1‑150 (MQ=255)                                                     gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcac                                                                                             >  4:116223/1‑150 (MQ=255)                                                        agaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATg                                                                                          <  1:104312/150‑1 (MQ=255)                                                            gTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGAG‑AAA‑AAAGATAACGACTATTTTTTTCATTTTGCAGTGACATATTATTCCTCACATGcaca                                                                                      >  2:474967/1‑150 (MQ=255)                                                              aGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacaca                                                                                    >  1:91341/1‑150 (MQ=255)                                                                 cAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                                 <  2:248677/150‑1 (MQ=255)                                                                  aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                                >  1:352675/1‑150 (MQ=255)                                                                   actgactgATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAA‑AAAGATAACGACTATTTTTATAATTTTGCAGTGCCATATTGTTTCTCACATGCACACATTGGt                                                                               >  3:13788/1‑150 (MQ=255)                                                                     tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                             >  1:416891/1‑150 (MQ=255)                                                                       actgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATg                                                                           >  1:403/1‑150 (MQ=255)                                                                              tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                    >  4:64633/1‑150 (MQ=255)                                                                              tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                    <  1:363955/150‑1 (MQ=255)                                                                                 agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTTTTTTTTTTCTGTTCCATATTGTTCGTCACATGCATACATTGGTAATGAAAAAAAGAc                                                                 >  2:65473/1‑150 (MQ=255)                                                                                 agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAc                                                                 >  3:329465/1‑150 (MQ=255)                                                                                       ttcttcGAGCAATATGTTTTCCTTTATTTTCCAAGAAGGACGGATATTCTGGTTTCCAGTCGCG‑AAA‑AAAGATAACGACTAATTTTCTCACTTTGCAGTCCCATATTGTTCCTCACTTCCACACAAAAAAAAAAAAAAAAAGACAAAaca                                                           <  4:178050/150‑1 (MQ=255)                                                                                          ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑CAA‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                        >  2:75194/1‑150 (MQ=255)                                                                                              aGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTa                                                    >  3:134594/1‑150 (MQ=255)                                                                                              aGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTa                                                    <  1:53133/150‑1 (MQ=255)                                                                                                cAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAg                                                  >  4:204912/1‑150 (MQ=255)                                                                                                cAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAg                                                  <  1:513780/150‑1 (MQ=255)                                                                                                  atatGTTTTCATTTATTTTACCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAATAGAAAAAAAAAAAAGACAAAACACGAGGTAAGgc                                                <  4:89277/150‑1 (MQ=255)                                                                                                  atatGTGTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTGTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGgc                                                >  4:158729/1‑150 (MQ=255)                                                                                                            cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                      <  2:291194/150‑1 (MQ=255)                                                                                                            cTTTATCTTCCCAGAAGGCCGGATTCTCTGGTTTCCAGTCGCG‑AAA‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCAATATCGTTCCACCCATGAAAACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                      <  2:499178/150‑1 (MQ=255)                                                                                                            cTTTAGTTGCCCAAGAGGCCGAATTTTCTGGTTTCCAGTCGCG‑AAA‑AAAAATAACGACTATTTTTCTCATTTTGCAGTGCCAAATTGTTCCTCACATGCACACATTGGTAATGAGAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                      <  2:391031/150‑1 (MQ=255)                                                                                                             tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTTTCTTTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                     >  2:36814/1‑150 (MQ=255)                                                                                                             tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCTCACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                     <  1:427324/150‑1 (MQ=255)                                                                                                             tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                     >  2:10151/1‑150 (MQ=255)                                                                                                                     cccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGa                             >  1:518765/1‑150 (MQ=255)                                                                                                                      ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGaa                            <  2:352675/150‑1 (MQ=255)                                                                                                                          aaGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGTCTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaa                        <  1:174321/150‑1 (MQ=255)                                                                                                                              ccGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGAt                    >  2:105650/1‑150 (MQ=255)                                                                                                                                  aTTTTCTGGTTTCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAAt                <  1:226915/150‑1 (MQ=255)                                                                                                                                     ttCTGGTTTCCAGTCGCG‑AAA‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAACAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg             >  3:379545/1‑150 (MQ=255)                                                                                                                                           tttCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAg       <  2:328743/150‑1 (MQ=255)                                                                                                                                             tCCAGTCGCGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACCTTGGTAATGAAGAAAAGACAAAACACGCGGTAAGGCGCAATAGCCAGTTATTAGAATTACGGATGGATTGGGTGAAGTg     >  1:312087/1‑150 (MQ=255)                                                                                                                                                aGTCGCG‑AAA‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTg  >  1:49498/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                 GATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTG  >  NC_000913_3_hsa_tpiA/1543872‑1544163

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