Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	264,098	G→T	100%	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	100.0%	67.5 / NA	28	L50L (CTC→CTA)	ykfH	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (16/12);  total (16/12)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTAT  >  NC_000913_3_hsa_tpiA/263956‑264212                                                                                                                                               |                                                                                                                   cTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACt                                                                                                              >  2:589018/1‑149 (MQ=255)  tCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTc                                                                                                             >  1:98025/1‑149 (MQ=255)     cAGGCGGAGTGGAATCTCCCGCTGCAGGTCCGTGGTAGAGTTTCTCGTGAGAGGTCGCCACTATTGTCAGTTCAGGGAGGTACTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCgg                                                                                                          <  2:495713/149‑1 (MQ=255)            aGTGTAATCTCCCGCTGCACGCCCCTGATAAGGTTGCTCCTTACAGTTCTCCACTACTGTCAGTTCAGGGTGATGCCCATCAGGCCGCTATAGGCCCCCTTCCGCTCCTAGCGGCGCTCCGCGCATCCCCGTAGGACTCCGGCGatacc                                                                                                   >  2:364650/1‑147 (MQ=255)                                  cccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATATGGCCCATTGCGGTCCTGGCGGCGTTCGGCGGATACCGCTAGGACTCCGGCGCTATCCGGAACGTCGCTGCCGGCGTaa                                                                             >  1:307560/1‑149 (MQ=37)                                       gTAGGGTTGCTCATTACAGTTCTCCACTATTGCCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGTTATCCGGAACGTCCCTGCCGGTGTACTcaca                                                                        >  2:311968/1‑149 (MQ=255)                                                ctgaTTAGAGTTCTCCACTATTTCCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGGGGTCCTGCCGGCGTTCGGCGCATACCGCTAGGATTCCGGCGATATCCGGAAGGTCCCTGCTGGTGTAATGACAGACGCTAcc                                                               <  2:42773/146‑1 (MQ=25)                                                   aTTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGCTATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTg                                                            >  1:268910/1‑149 (MQ=255)                                                    ttACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCGCTGCCGGTGTAATGGCAGACGCTACCGTGc                                                           >  1:167378/1‑149 (MQ=255)                                                        aGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACCGGCGCTACGGTAccccc                                                       >  2:338602/1‑146 (MQ=25)                                                         tttCTCCACTGTTGTCAGTTCAGCGGGCTGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGTGGCGTTCGGCGTATACCGCTAGGTCTCGGGCGATATCCGGAACGGCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTg                                                      <  2:395362/148‑1 (MQ=21)                                                         gTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTg                                                      <  1:39502/149‑1 (MQ=255)                                                               ccatATTGGCAGTTCAGGGTGATGGCCATCGGGGAGGTAAAGGGCCCAGTGCGGTCCTGGCGGCGTTCTGCGTATACCGCTAGGACCCCGGCGATATCCGGAACGTCGCTGCAGGTGTAATGACAGACGCTACCGTGCCACTGGTAttt                                                <  2:99282/146‑1 (MQ=21)                                                                           ttCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCCGTa                                    >  1:40339/1‑149 (MQ=255)                                                                           ttCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTAGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTACCTACCGGTGTAATGACAGACGCTCACGTGCCACTGGTATTTGCCGGTGCAGTa                                    >  2:355043/1‑149 (MQ=21)                                                                            tCAGGGTGATCCTCATCAGGCAGGTATATGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAg                                   <  1:119838/149‑1 (MQ=255)                                                                                        tcatcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCgg                       <  1:432745/149‑1 (MQ=255)                                                                                        taatcaGGGAGGTTTCGGGCTCAGTGCGGTGCTGGGGGCGTTCGGCGTATACCGGTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGGGCAATGATAGACGCTACCGTGCCACTGGTATTTGCCGGTGTAGTAGCGAAAGATTCgg                       <  2:303209/147‑1 (MQ=255)                                                                                         catcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTGGCCGGTGCAGTGGCGAAAGATTCggg                      >  2:2572/1‑149 (MQ=255)                                                                                          atcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGa                     >  1:69171/1‑149 (MQ=255)                                                                                          atcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGGAGTGGCCGGTGCAGTAGCGAAAGATTCAGGa                     >  1:7058/1‑149 (MQ=21)                                                                                            caggcaggTATAGGGCCCATTGCGGTCCTTGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCGCAAGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACt                   <  1:653306/149‑1 (MQ=21)                                                                                             aggcaggTATAGGGCCCATTGCGGTTCTGGCGGCGTTCGCCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCgggccag                  <  2:246004/149‑5 (MQ=255)                                                                                                     ataGGGCCTATTGCGGTCCTGGCGGCGTTCGGGGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCt          <  2:247077/149‑1 (MQ=255)                                                                                                        gggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGc       <  1:638871/149‑1 (MQ=255)                                                                                                          gcaagttGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGTTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCgg     >  2:441801/6‑149 (MQ=255)                                                                                                           cccATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGt    >  1:302343/1‑149 (MQ=255)                                                                                                             cATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGtat  >  2:174491/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                   CTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTAT  >  NC_000913_3_hsa_tpiA/263956‑264212

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