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	100.0%	102.8 / NA	34	L50L (CTC→CTA)	ykfH	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (17/17);  total (17/17)

CCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATTT  >  NC_000913_3_hsa_tpiA/263959‑264243                                                                                                                                            |                                                                                                                                                  ccAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCg                                                                                                                                          >  1:275851/1‑149 (MQ=255)  cAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCgg                                                                                                                                         <  1:516067/149‑1 (MQ=255)  cAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGAGGTATTGTCTTCGGTCGGCGTATACCGCTAGGACTCCgg                                                                                                                                         <  1:593219/149‑1 (MQ=255)   gggCGCAGTGTGTCCTCCCGCTGCGGGCCCCTGGTAGGTTTGCTCATTACTGTTCTCCACTATTGTCAGTTCAGGGTGATGCGCATCAGTCAGGTATAGGGCCCTTTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGc                                                                                                                                        <  2:115821/148‑1 (MQ=255)     gcgcAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATGACAGTGCTCCACTATTGTAAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGa                                                                                                                                      <  2:26722/149‑1 (MQ=255)           tttAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCg                                                                                                                                <  2:311057/147‑1 (MQ=255)                tctcCCGCTCCAGGCCCCTGGTAGGGTTGCCCATTCCAGTTCTCCACTATTGTCAGTTCAGGGTGCTGCTCATCAGGCAGGTATAGGGCCCGTTGCGGTCCTGCCGGCGTCCGCCGTATCCCGCTTGGACTCCGGCGCTCTCCGGGACg                                                                                                                           >  2:137003/1‑149 (MQ=255)                 ctcCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGt                                                                                                                          >  1:596151/1‑149 (MQ=255)                  tcCCGCTGCAGGCCCCTGGTAGGGTAGCACATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTAGGCGTATACCGCTAGGACTCCGCCGATATCCGGAACGTc                                                                                                                         >  2:491059/1‑149 (MQ=255)                    ccGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTccc                                                                                                                       <  2:219623/149‑1 (MQ=255)                     cGCTGCAGGCCCCTGGTAGGGCTGCCCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCt                                                                                                                      <  1:270543/149‑1 (MQ=255)                       cTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGc                                                                                                                    <  2:596151/149‑1 (MQ=255)                        tGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGcc                                                                                                                   <  1:167397/149‑1 (MQ=255)                          cAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGcc                                                                                                                   >  1:207493/1‑147 (MQ=255)                          cAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGcc                                                                                                                   <  2:207493/147‑1 (MQ=255)                               cccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGAGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCGCTGCCggaggca                                                                                                            >  1:278787/1‑144 (MQ=255)                                                     aGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCAGCGTATACCGCTAGGACTCCGGCGATATCCGCAACGTCCCTGCCGGTGTAATGACAGACCCTACCGTGCCACt                                                                                      >  2:498893/1‑149 (MQ=255)                                                                 ttGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCgg                                                                          >  2:182378/1‑149 (MQ=255)                                                                    tCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGc                                                                       >  2:213063/1‑149 (MQ=255)                                                                        ttCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTa                                                                   <  1:351254/149‑1 (MQ=255)                                                                                gATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGa                                                           <  2:89546/149‑1 (MQ=255)                                                                                     tgatcaGGCGGGTATGGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCgg                                                      <  2:275851/147‑1 (MQ=255)                                                                                       atcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGa                                                    >  1:500517/1‑149 (MQ=255)                                                                                              aggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGCCAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGCAGCGAAAGATTCGGGGCTGAGGa                                             >  2:437996/1‑149 (MQ=255)                                                                                              aggTATAGGGCCCATTGCGGTCCTGGCGGCGTCCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGGCGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGCCGCAggg                                             >  1:27429/1‑148 (MQ=18)                                                                                                gTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATg                                           >  2:393295/1‑149 (MQ=255)                                                                                                   taGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTg                                        >  2:586959/1‑149 (MQ=255)                                                                                                       gCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGa                                             <  1:502611/140‑1 (MQ=255)                                                                                                       gCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGa                                             >  2:502611/1‑140 (MQ=255)                                                                                                         ctatggCGGTCCTGGCGGCCTTCGGCGCTGACCGCTAGGCCTCCGGCAATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGTTATTTGCCGGTGCAGTAGCGAAAGATTCGGGCCTCAGGATGCTGGCGGta                                  <  1:154926/144‑1 (MQ=14)                                                                                                         ccATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGta                                  <  1:426822/149‑1 (MQ=255)                                                                                                                  tCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGGACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGCCTCAGGATGCTGGCGGTATATCGTCAt                         >  1:102148/1‑149 (MQ=255)                                                                                                                             ttCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGtt              <  1:182378/149‑1 (MQ=255)                                                                                                                                    tataCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGAt       >  1:311486/1‑149 (MQ=255)                                                                                                                                         cgctAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAttc  <  1:393295/149‑2 (MQ=255)                                                                                                                                            |                                                                                                                                                  CCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATTT  >  NC_000913_3_hsa_tpiA/263959‑264243

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