Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	1,775,562	Δ1 bp	intergenic (+118/‑25)	aroD → / → ydiF	3‑dehydroquinate dehydratase/putative acetyl‑CoA:acetoacetyl‑CoA transferase: alpha subunit/beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	1,775,559	0	A	.	100.0%	117.1 / NA	26	intergenic (+115/‑28)	aroD/ydiF	3‑dehydroquinate dehydratase/putative acetyl‑CoA:acetoacetyl‑CoA transferase: alpha subunit/beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (10/16);  total (10/16)

ACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAACCACCTCGTATTAATGGACGAGTGCCGGTCCTGTCGGCACAGGAAGCGGTGAATTATATTCCCGACGAAGCAACACTTTGTGTGTTAGGCGCTGG  >  NC_000913_3_hsa_tpiA/1775424‑1775696                                                                                                                                        |                                                                                                                                          aCTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGG‑AAATCACATCTGaa                                                                                                                             >  1:129386/1‑149 (MQ=255) aCTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCCGTGCGCTAATTCTACATTTGTCGTAGGG‑AAATCACATCTGaa                                                                                                                             >  2:286020/1‑149 (MQ=255) aCTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCACCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGG‑AAATCACATCTGaa                                                                                                                             <  1:79481/149‑1 (MQ=255)          cacCAGGCAAAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCGCTGTGCGCTAATTCTCCATTTGGCGTAGGG‑AAATCACATCTGAATCAGGAAtt                                                                                                                    <  2:200995/149‑1 (MQ=255)           acCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGG‑AAATCACATCTGAATCAGGAATTa                                                                                                                   >  1:284986/1‑149 (MQ=255)                    aaGCAATAATATTTCGGGGGGAGCCCGCTCACCGGCGAACTAAAAAATATATTTAATCGTATTTAATAAAAATATTTCTTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGG‑AAATCACCTCTGAATCAGGATTTAACACTGaaa                                                                                                          <  2:365088/149‑1 (MQ=255)                            atatTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGG‑AAATCACATCTGAATCAGGAATTAACAATGAAACCTGTaaa                                                                                                  >  1:32042/1‑149 (MQ=255)                               tttCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGTCGTAGGG‑AAATCACATCTGAATCACGAATTAACAATGAAACCTCTAAAacc                                                                                               >  2:336012/1‑149 (MQ=255)                                       ggAATACCCTCCCTGCCTAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGG‑AAATCAGATCTGAATCAGGAAATAACAATGAAACCTGTAAAACCACCTCGTa                                                                                       <  2:198350/149‑1 (MQ=255)                                                tCCCCACCTAACTAAAAAATATACACAATCGGATTTAATAAAAATATTTCGAGAGTCTCTGTGCGCTAATTCTCCATTTGACGTACGG‑AAATCACATCTGAATCAACAACTAACAATTAAACCTGTAAAACCACCTCGTATTAATGGAc                                                                              <  2:65636/149‑1 (MQ=255)                                                        gAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGG‑AAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAACCACCTCGTATTAATGGACGAGTGCCg                                                                      <  1:314623/149‑1 (MQ=255)                                                                aaaTATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGG‑AAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAACCACCTCGTATTAATGGACGAGTGCCGGTCCTGTc                                                              <  1:111554/149‑1 (MQ=255)                                                                       ttCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGG‑AAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAACCACCTCGTATTAATGGACGAGTGCCGGTCCTGTCGGCACAg                                                       >  2:362736/1‑149 (MQ=255)                                                                             cGTATTTAATAAAAATATTTCGTGAGTCTCTGGGCGCTAATTCTCCATTTGGCGTAGGG‑AAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAACCACCTCGTATTAATGGACGAGTGCCGGTCCTGTCGGCACAGGAAGCg                                                 <  2:129386/149‑1 (MQ=255)                                                                                  ttAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGG‑AAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAACCACCTCGTATTAATGGACGAGTGCCGGTCCTGTCGGCACAGGAAGCGGTGaa                                            >  2:18728/1‑149 (MQ=255)                                                                                  ttAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGG‑AAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAACCACCTCGTATTAATGGACGAGTGCCGGTCCTGTCGGCACAGGAAGCGGTGaa                                            <  1:292993/149‑1 (MQ=255)                                                                                    aataaAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGG‑AAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAACCACCTCGTATTAATGGACGAGTGCCGGTCCTGTCGGCACAGGAAGCGGTGAAtt                                          <  2:103245/149‑1 (MQ=255)                                                                                      taaAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGG‑AAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAACCACCTCGTATTAATGGACGAGTGCCGGTCCTGTCGGCACAGGAAGCGGt                                               >  1:253197/1‑142 (MQ=255)                                                                                      taaAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGG‑AAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAACCACCTCGTATTAATGGACGAGTGCCGGTCCTGTCGGCACAGGAAGCGGt                                               <  2:253197/142‑1 (MQ=255)                                                                                                      gTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGG‑AAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAACCACCTCGTATTAATGGACGAGTGCCGGTCCTGTCGGCACAGGAAGCGGTGAATTATATTCCCGACGAAGCaa                        <  1:175478/149‑1 (MQ=255)                                                                                                           ggtgCGTTAATTCTCCATTTGGCGTAGGG‑AAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAACCACCTCGTATTAATGGACGAGTGCCGGTCCTGTCGGCACAGGAAGCGGTGAATTATATTCCCGACGAAGCAACACtt                   <  2:129438/148‑1 (MQ=255)                                                                                                             tgCGCTAATTCTCCATTTGGATTAGGG‑AAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAACCACCTCGTATTAATGGACGAGTGCTGGTCCTGTCGGCACAGGAAGCGGTGAATTATATTCCCGACGAAGCAACACTTtg                 <  1:41282/149‑1 (MQ=255)                                                                                                                 cTAATTCTCCATTTGGCGTAGGG‑AAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAACCACCTCGTATTAATGGACGAGTGCCGGTCCTGTCGGCACAGGAAGCGGTGAATTATATTCCCGACGAAGCAACACTTtgtgtg             <  2:32042/149‑1 (MQ=255)                                                                                                                        tcCATTTGGCGTAGGG‑AAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAACCACCTCGTATTAATGGACGAGTGCCGGTCCTGTCGGCACAGGAAGCGGTGAATTATATTCCCGACGAAGCAACACTTTGTGTGTTAGgcg      >  2:253931/1‑149 (MQ=255)                                                                                                                        tcCATTTGGCGTAGGG‑AAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAACCACCTCGTATTAATGGACGAGTGCCGGTCCTGTCGGCACAGGAAGCGGTGAATTATATTCCCGACGAAGCAACACTTTGTGTGTTAGgcg      >  1:129863/1‑149 (MQ=255)                                                                                                                            tttGGCGTAGGG‑AAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAACCACCTCGTATTAATGGACGAGTGCCGGTCCTGTCGGCACAGGAAGCGGTGAATTATATTCCCGACGAAGCAACACTTTGTGTGTTAGGCGCTgg  <  1:18728/149‑1 (MQ=255)                                                                                                                                        |                                                                                                                                          ACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATATTTCGTGAGTCTCTGTGCGCTAATTCTCCATTTGGCGTAGGGAAAATCACATCTGAATCAGGAATTAACAATGAAACCTGTAAAACCACCTCGTATTAATGGACGAGTGCCGGTCCTGTCGGCACAGGAAGCGGTGAATTATATTCCCGACGAAGCAACACTTTGTGTGTTAGGCGCTGG  >  NC_000913_3_hsa_tpiA/1775424‑1775696

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