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	97.3%	104.3 / NA	37	L50L (CTC→CTA)	ykfH	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (18/18);  minor base C (1/0);  total (19/18)
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

CGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTAC  >  NC_000913_3_hsa_tpiA/263952‑264234                                                                                                                                                   |                                                                                                                                         cGGGCTCGCAGGGGCAGTGTAATCCCCCGCTGCAGGCCACTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAgg                                                                                                                                       <  1:532792/150‑1 (MQ=255)      tCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTcc                                                                                                                                  >  1:6021/1‑150 (MQ=255)                    tAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGaa                                                                                                                    <  2:38426/150‑1 (MQ=255)                       tctcCCGATGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGGACgc                                                                                                                 >  1:636019/1‑149 (MQ=255)                         tcCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGGTCAGGGTGATGCTCATCAGACAGGTATCGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTACACCGCTAGGACTCCGGCGCTATCCGGGACGTcc                                                                                                               >  1:70750/1‑150 (MQ=255)                                   ggCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgtgt                                                                                                     >  1:640291/1‑150 (MQ=255)                                       ccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGaa                                                                                                                    <  2:566058/131‑1 (MQ=255)                                       ccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGaa                                                                                                                    >  1:566058/1‑131 (MQ=255)                                             ggggTTGCTTATTACAGTTCTTCATTTTTGTCAGTTCAGGGTGATGCTCATCGGGCAGGTATAGGGCCCTTTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTTCGGCGATATCCGGAACGTCGCTGCCGGTGTAATGACAGAc                                                                                           <  2:70249/149‑1 (MQ=25)                                                 ttGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATCCCGCTAGGACTCCGGCGCTACCCGGAACGTCCCTGCCGGTGTACGGACAGACGCTa                                                                                       >  2:25597/1‑150 (MQ=255)                                                   gCTCATTAAAGTTCTCCACTATTGTCAGGTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTAcc                                                                                     <  1:30436/150‑1 (MQ=255)                                                            aGTTCTCCACTATTGTCAGGTGAGGGTGATGCTCATCAGGCAGGGATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTg                                                                            <  2:455221/150‑1 (MQ=255)                                                                 tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTAtt                                                                       >  1:263730/1‑150 (MQ=255)                                                                 tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGCGCCGTGACAGCAGCTCGCGTTCCACTGGTAtt                                                                       >  2:121033/1‑150 (MQ=255)                                                                  ccACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGTAACGTCCCTGCCGATGTAATGACAGACGCTACCGTGCCACTGGTAttt                                                                      >  2:211529/1‑150 (MQ=255)                                                                   cACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCCTGCCACTGGTAtttc                                                                     >  2:425495/1‑149 (MQ=255)                                                                      tataGTCAGTTGAGGGTGCAGATCATACGGTAGGTATAGGGCCCATTGCGGGCCTGGCGTCGTTCGTTGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTGGCCg                                                                  <  2:479212/150‑1 (MQ=21)                                                                          gTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTAGGTGTATAGCGCCAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTCTAATGACACACGCTACCCTGCCACTGGTATTTGCCGGTGc                                                              >  2:284636/1‑150 (MQ=255)                                                                               ttCAGGGTGATGCTCAACAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAg                                                         >  1:600682/1‑150 (MQ=255)                                                                                    ggTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCCGTAGCGcaa                                                    >  2:255562/1‑150 (MQ=255)                                                                                      tGATGCTCATCAGGCAGGTATAGGGTCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGCAGCGGCCCTGCCGGTGACGTAGCATTCGCGACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGa                                                  <  1:677641/150‑1 (MQ=21)                                                                                          gCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCg                                              >  2:618285/1‑150 (MQ=255)                                                                                          gCTAATTAGGCAGGGATAGGGCCCCTTGCGGTCCTGGGGGCGTTCGGCGTAGACCGCTAGGACTCGGGCGATATGAGGACCGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCg                                              <  2:70750/150‑1 (MQ=18)                                                                                           cTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGGGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCgg                                             <  2:6021/150‑1 (MQ=255)                                                                                           cTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCgg                                             <  2:61973/150‑1 (MQ=255)                                                                                                   gcaggTATAGTGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAgg                                     >  1:114884/1‑150 (MQ=255)                                                                                                   gcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAgg                                     >  1:285817/1‑150 (MQ=255)                                                                                                    cagtTATAGGGCCCATTGCGGTCCGGGCGGCGTTCGGCGTATACGGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGCCAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGa                                    <  1:211529/150‑1 (MQ=18)                                                                                                       gTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGc                                 >  1:232969/1‑150 (MQ=255)                                                                                                                ccATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGtat                        <  1:618285/150‑1 (MQ=255)                                                                                                                ccATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGtat                        <  2:367889/150‑1 (MQ=255)                                                                                                                ccATGGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGtat                        <  1:137688/150‑1 (MQ=255)                                                                                                                             ggcggcGCTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGtcc           <  1:185993/150‑1 (MQ=255)                                                                                                                              gcggcgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTGCTTGCCCt          >  2:422388/1‑150 (MQ=255)                                                                                                                                      cGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTAc  <  2:114884/150‑1 (MQ=255)                                                                                                                                              aCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGc                                                                    >  2:202935/1‑76 (MQ=255)                                                                                                                                              aCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGc                                                                    <  1:202935/76‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                         CGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTAC  >  NC_000913_3_hsa_tpiA/263952‑264234

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