Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	2,204,477	G→T	intergenic (+188/‑119)	yehK → / → yehL	uncharacterized protein/putative hexameric AAA+ MoxR family ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	2,204,477	0	G	T	100.0%	116.9 / NA	34	intergenic (+188/‑119)	yehK/yehL	uncharacterized protein/putative hexameric AAA+ MoxR family ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (16/18);  total (16/18)

CAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCA  >  NC_000913_3_hsa_tpiA/2204334‑2204621                                                                                                                                                |                                                                                                                                                 cAATGGCCTGGATTATAAAGCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTtttttt                                                                                                                                            <  2:627297/150‑1 (MQ=255)  aaTGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGtttttttt                                                                                                                                           >  2:38328/1‑150 (MQ=255)             ttATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGa                                                                                                                                >  1:490789/1‑150 (MQ=255)               ataAATCCCATTATCTTAATTGCAACGGGGTCAAGCACAGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATc                                                                                                                              <  1:797672/150‑1 (MQ=255)                        aTTATCTTAATTGCAACGGGGTCCATCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCa                                                                                                                     <  1:314347/150‑1 (MQ=255)                           aTCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGat                                                                                                                  <  2:409115/150‑1 (MQ=255)                            tCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGata                                                                                                                 >  1:989469/1‑150 (MQ=255)                                 atggCACCGGCGTCCAGACGTGTGTATGATTCGTTCGTTGGAACACGATGTTGCACAGGCTTTGGTGTTAGGCCTTAAAATAGTAAATAAAAGAACCGGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGATAttttt                                                                                                            <  2:295877/147‑1 (MQ=255)                                 aTTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGATAttttt                                                                                                            <  1:38328/150‑1 (MQ=255)                                  ttGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGATAtttttt                                                                                                           >  2:699476/1‑149 (MQ=255)                                  ttGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTc                                                                                                           >  1:771530/1‑150 (MQ=255)                                    gCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATACAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTCCTGTGATCTGCTTTCCAGATATTTTTcta                                                                                                         >  2:255032/1‑148 (MQ=255)                                         ggggTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGata                                                                                                                 >  1:159388/1‑137 (MQ=255)                                         ggggTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGata                                                                                                                 <  2:159388/137‑1 (MQ=255)                                            gTCCAGCCGTGGTTATAATACGTTCGTTAAAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAaacaa                                                                                                 <  1:381121/150‑1 (MQ=255)                                             tCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAaacaac                                                                                                <  2:490789/150‑1 (MQ=255)                                                   cGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATgc                                                                                          <  2:458558/150‑1 (MQ=255)                                                    gTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATgcg                                                                                         >  1:412317/1‑150 (MQ=255)                                                     tGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATgcgc                                                                                        <  2:599912/150‑1 (MQ=255)                                                           aaaTGCGTTCGTAAGAAGAAGATAATGCACAGGCTGTGGTGTTAGGCGTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTCCTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACa                                                                                  <  2:962119/149‑1 (MQ=255)                                                               tcgttcgtTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTa                                                                              <  2:921606/150‑1 (MQ=255)                                                                     gttGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGCTATTTTTCGCTCAAACAACTAATGCGCCAAAAATTTATTgcgc                                                                        >  1:342390/1‑150 (MQ=255)                                                                                aTGTTGCACAGCTTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTACTTTCCAGATATTTCTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCg                                                             <  1:901770/150‑1 (MQ=255)                                                                                         aGGCTGTGGCGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTATTTTTATCTGCTTTCCAGATATTTTTCTCTCAATCATCTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCa                                                    >  2:609299/1‑150 (MQ=255)                                                                                          ggCCGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCTCGTAAAATATCGTTTATTTCAt                                                   >  2:994696/1‑150 (MQ=255)                                                                                                    ttAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTGACCAGCAGTTTTTTTTTTTTATTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAGTTATTGTGCGTAACATTTTTTGCGCGTAAATTATCGTTTTTTTTTTTTATACtttt                                         >  2:897857/1‑150 (MQ=255)                                                                                                        gCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAgg                                     >  1:980529/1‑150 (MQ=255)                                                                                                                   gTAAATAAAAGAACCTGCCTCTCCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAatatat                          <  1:612264/150‑1 (MQ=255)                                                                                                                    taaataaaAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTCTTTCATTCATACATTTCAGGGATGAatatatt                         >  2:861175/1‑149 (MQ=255)                                                                                                                    taaataaaAGAACCTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATg                         >  1:727790/1‑150 (MQ=255)                                                                                                                                 cTGCCTCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGaata            >  2:630176/1‑150 (MQ=255)                                                                                                                                    ccTCACCAGCATTTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAatc         <  2:412317/150‑1 (MQ=255)                                                                                                                                      tCACCAGCAGTTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAatcac       >  1:604108/1‑149 (MQ=255)                                                                                                                                           agcagTTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCa  <  2:272540/150‑1 (MQ=255)                                                                                                                                           agcagTTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCa  <  1:255032/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                 CAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCA  >  NC_000913_3_hsa_tpiA/2204334‑2204621

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