Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	264,098	G→T	L50L (CTC→CTA)	ykfH ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	264,098	0	G	T	95.8%	63.3 / NA	24	L50L (CTC→CTA)	ykfH	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (12/11);  minor base A (1/0);  total (13/11)
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

GCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAT  >  NC_000913_3_hsa_tpiA/263963‑264241                                                                                                                                        |                                                                                                                                                gcgcAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCAAGGACTCCGGCGa                                                                                                                                    >  2:207873/1‑149 (MQ=255)            tctcCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACACTTCCCCACTATTATCAGTTCAGGGTGATGCTCATGAGGCAGGTATAGGGACCATTGCGGTCCTGGCGGCGTTCGCCGTATACCGCTAGGACTCCGGCGATCTTCGGACCg                                                                                                                         >  2:271894/1‑149 (MQ=255)               cccGCTGGCGGCCCCTGGTAGGGGTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGGCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTcc                                                                                                                      <  1:396460/149‑1 (MQ=255)                ccGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTccc                                                                                                                     >  1:286427/1‑149 (MQ=255)                 cgggggaGGCCCTTGGTCGTGTTTCTCATTACAGTTCTCAACTATTGTCAGTTGAGGGTGTTGCTCAGCAGGCAGGCAAAGGGCTCATTGCGGTCCTGGCGGCGTTCGGCGTATATCGCTAGGACTCCGGCGATATCCGGAACGTCCCt                                                                                                                    <  2:38039/143‑1 (MQ=255)                        ggCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgtg                                                                                                             >  1:72200/1‑149 (MQ=255)                        ggCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgtg                                                                                                             >  1:168611/1‑149 (MQ=255)                               ggTAGGGTTGCTCATTACAGTTCGACACTACTGTCAGTTTCGGGTGTTCCGTACGAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGAc                                                                                                      <  1:96001/149‑1 (MQ=25)                                   gggTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGAc                                                                                                  <  1:418094/149‑1 (MQ=255)                                        gCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGCATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGACGGTGTAATGACAGACGCTAc                                                                                             >  1:275816/1‑149 (MQ=255)                                                  ttaccacACTAAAGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCTCATTGCGGTCCTGGCGGCGTTCGGCGAATACCGCTAGAACCCCGGCGATACCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTg                                                                                   <  2:398075/143‑1 (MQ=21)                                                    tctcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGt                                                                                 <  1:281323/149‑1 (MQ=255)                                                      tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCTGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTAt                                                                               <  1:24832/149‑1 (MQ=255)                                                      tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCGGGAGCGTCTCTACGGGTCTAAGGACACACGCTCCCCATCCACTGGAAt                                                                               >  2:188048/1‑149 (MQ=255)                                                                     tCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGTCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTACCACTGGTATTTGCCGGTGCAGTAg                                                                >  1:379112/1‑149 (MQ=255)                                                                           tGATGCTCATCAGGCAGGTATAGGGCCCATTGCGATCCTGACGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCGGGAACGTCCCTGGCGGTGTAATGACAGACGCTACCCCGCCACTGGTATTTGCCGGTGGAGTAGCGATAg                                                          >  2:221660/1‑149 (MQ=21)                                                                               gCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCAGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTc                                                      >  1:333844/1‑149 (MQ=255)                                                                                                       aTTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTACCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGtata                              >  2:387482/1‑149 (MQ=255)                                                                                                                    cggcgTTCGGCGTATACCCCTAGGACTCCGGCGCTATCCGGAACTTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCt                 <  2:275816/149‑1 (MQ=255)                                                                                                                     ggcgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTg                <  1:252045/149‑1 (MQ=255)                                                                                                                        gTTCGGCGTATACCGCTAGGTATCCGGGGATATTCGGAACGTACCTGCCGGTGGAATGACAGACGCTACCGTGCCACTGGTAGTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTGATTGCCCTGCGt             <  2:59782/149‑1 (MQ=14)                                                                                                                        gTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGt             >  2:327549/1‑149 (MQ=255)                                                                                                                            ggCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCGACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTAc         <  1:387482/149‑1 (MQ=255)                                                                                                                                   aCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATact  >  1:206752/1‑147 (MQ=255)                                                                                                                                        |                                                                                                                                                GCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAT  >  NC_000913_3_hsa_tpiA/263963‑264241

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