Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	1,544,022	C→A	L13F (TTG→TTT)	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,022	0	C	A	88.5%	73.8 / ‑6.1	27	L13F (TTG→TTT)	narU	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base A (11/12);  total (13/14)
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

ACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAA  >  NC_000913_3_hsa_tpiA/1543894‑1544157                                                                                                                                  |                                                                                                                                         aCTGCGCTAAATAGCATCTAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGCTTTCTTCGAGCGATATGTTTTCCTTTATTTGCCCATAAGGCCGGATTTTCTAGTTTCCAGTCGCG‑CAAAAAGATAACGATTAtttt                                                                                                                      >  2:430410/1‑149 (MQ=255)                        acaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCTAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTg                                                                                               >  1:247105/1‑149 (MQ=255)                           gagggCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                            <  2:527339/146‑1 (MQ=255)                           gAAGGAAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                            <  1:494314/149‑1 (MQ=255)                              ggCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTc                                                                                         <  2:445812/149‑1 (MQ=255)                                  agaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTGTTTTCCCAGAAGGCCGGATTTGCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAt                                                                                     >  1:540493/1‑149 (MQ=255)                                         gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacaca                                                                              >  1:534789/1‑149 (MQ=255)                                           cAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAtt                                                                            >  2:112613/1‑149 (MQ=255)                                                          cAGAGATGTCTGCTAGCAATATGTTGTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTGTCATTTTGCAGTGGCCAAGTTTGCCTCACATGCACCCATTGGTAATGAAAAAAAg                                                             >  2:559957/1‑149 (MQ=255)                                                           agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAACAGGCCGGATTTTCTGGACTCCTGCCGCG‑CAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATTAAACAAAGa                                                            >  1:522508/1‑149 (MQ=255)                                                            gagaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAc                                                           >  2:226356/1‑149 (MQ=255)                                                                tttCTTAGAGCAAAATGTTTTCCTTTATTTTACCAGAAGAAAGAAAAATCTGGTTAACAGTCGCG‑AAAAAAGATAACGACTATTTTCCTCATTTTGCAGGGCCATATTGTTCCTCACATGCACACAAAAAAAAAAAAAAAAAGACaaaa                                                       <  2:491052/149‑1 (MQ=255)                                                                 ttcttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCA‑TAAAAAGATAACGACTACGTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAac                                                      <  1:527819/149‑1 (MQ=255)                                                                        aGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGCAAAAAAGACAAAACACGAGGt                                               >  1:58605/1‑149 (MQ=255)                                                                         gCAATTTGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTATGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTa                                              <  2:352398/149‑1 (MQ=255)                                                                         gCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTa                                              <  1:218305/149‑1 (MQ=255)                                                                                gTTTTCCTTTATTTCCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCa                                       <  1:76507/149‑1 (MQ=255)                                                                                       tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAAACCAAACACGAGGTAAGGCGCAATAGCCa                                >  2:24556/1‑149 (MQ=255)                                                                                          aTTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGtt                             >  1:283818/1‑149 (MQ=255)                                                                                                  aGACGGCCGGATTCTCTGGTTTCCAGTCGCGAC‑ACAAGATCCCGTCTTTTTTTCTCATTTTGCAGTGCCACATTGTCCCTCACATGCACAAATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAAt                     <  1:559957/149‑1 (MQ=255)                                                                                                      ggCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg                 <  1:332710/149‑1 (MQ=255)                                                                                                      ggCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg                 <  1:306418/149‑1 (MQ=255)                                                                                                       gccgttttttCTTTTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGTTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAgg                >  2:36115/7‑149 (MQ=255)                                                                                                                     tttCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCTTCACATGCACACATTGGAAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGac  <  2:540493/149‑2 (MQ=255)                                                                                                                     tttCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGc                                  >  2:410741/1‑117 (MQ=255)                                                                                                                     tttCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGc                                  <  1:410741/117‑1 (MQ=255)                                                                                                                     tttCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGaa  <  1:112613/149‑1 (MQ=255)                                                                                                                                  |                                                                                                                                         ACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAA  >  NC_000913_3_hsa_tpiA/1543894‑1544157

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