Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	1,197,502	C→T	G11R (GGA→AGA)	ymfD ←	e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	1,197,502	0	C	T	100.0%	38.5 / NA	16	G11R (GGA→AGA)	ymfD	e14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/5);  total (11/5)

ATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGA  >  NC_000913_3_hsa_tpiA/1197360‑1197646                                                                                                                                               |                                                                                                                                                 aTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCat                                                                                                                                            >  2:11539/1‑149 (MQ=255)         atcaGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGAGGCCATTTTTCTCAATGCTTCTAATATACTCTTTGGGGAGGCTAGAGGGCATCGAATAAAGTTTTGCGGCATTCCCTGAGCGAATATTAACTCTGTGCATATTATAGt                                                                                                                                    >  2:151696/1‑149 (MQ=255)          tcaGTATATCTAAGGGTTCCGCAGCCAAAATCAAGGGGGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGCTCCTGTTGCAGTAACTGAGCGAATATTAACTCTGTGCATATTATAGtt                                                                                                                                   >  2:426341/1‑149 (MQ=255)                        ttttCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATAACTAGAAGGCATCGTATGAGGTTTTGCTGCATTATCTGAGCGAATATTAACTCTGTGCATATTATATTTCAAAGCAACTAGCa                                                                                                                     >  2:141126/1‑149 (MQ=255)                                                    ccATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGCCATCGTATGAGGTTTTGCTGCATTCTCTGAGCTAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGaaaaa                                                                                         <  1:202127/149‑1 (MQ=255)                                                                 ttcTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttatttta                                                                            <  1:151925/147‑1 (MQ=255)                                                                                          tCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATCTTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGCACAATTCTTCTCGTTTATCGTCCGCGATTTTTCCgacttgac                                                   >  2:231600/1‑149 (MQ=255)                                                                                             agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGCATTTTTCCGACTTCTCttt                                                >  2:166231/1‑149 (MQ=255)                                                                                                 ggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCTCCTTGTCTTTCAGt                                            >  2:155903/1‑149 (MQ=255)                                                                                                      cGTATGAGGTTTTGCTGACATCTCTGAGCGATTATTAACTCTGTGCATACTATAGTTCAAAGCCAGTACCCTTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTgg                                       <  1:426341/149‑1 (MQ=255)                                                                                                         aggAGGTTTTGCTGCATTCTCTGAGCGAATAATAACTCTGCGCATATTATAGTTCATCGCCAAAACAATTCCATTTATGTTATCTTTAGTAAACCATTCTTTTAATTTAACGTTCGGTCTTTTTTTTCCTTTCCTTGAAGTTCTGGGCt                                    >  2:419441/3‑149 (MQ=255)                                                                                                          acagatgttGGTGAATTCTCTGAGCGAATATTCACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCACTTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTATCTTTCAGTTCTCGGAtt                                   >  2:368224/8‑149 (MQ=255)                                                                                                                       cATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCga                                                         <  1:486660/114‑1 (MQ=255)                                                                                                                       cATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCga                                                         >  2:486660/1‑114 (MQ=255)                                                                                                                                       atTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCATTTCTGTGCTTAACTGTCAAAGCAAAAATAATACcactca      >  2:586015/1‑149 (MQ=255)                                                                                                                                           aCTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTga  <  1:297673/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                 ATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGA  >  NC_000913_3_hsa_tpiA/1197360‑1197646

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