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%	132.1 / NA	41	G11R (GGA→AGA)	ymfD	e14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/20);  total (21/20)

GATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACT  >  NC_000913_3_hsa_tpiA/1197359‑1197644                                                                                                                                                |                                                                                                                                               gATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCat                                                                                                                                          <  1:845042/150‑1 (MQ=255)          atcaGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGtt                                                                                                                                 <  2:439600/150‑1 (MQ=255)           tcaGAATATCTAAGTTTTCCGCAGCCAAAATCAAGAGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGtt                                                                                                                                 >  2:190465/1‑149 (MQ=255)           tcaGAATATCTAAGTTTTCCGCAGCCAAAATCAAGAGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGtt                                                                                                                                 <  1:190465/149‑1 (MQ=255)              gAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCTGCTCTTCTAAGATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCaaa                                                                                                                             <  1:29255/150‑1 (MQ=255)                  atCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATGTTTCTCAATGCTTCTAATATACTCGCAAAGATATCGAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCaa                                                                                                                         <  2:626714/150‑1 (MQ=255)                             ccGAAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATCTTCCTAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCaa                                                                                                              <  1:237443/150‑1 (MQ=255)                                  gCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTAACGGAGCGAATATTAACGCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt                                                                                                         <  2:117286/150‑1 (MQ=255)                                    cAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGt                                                                                                       >  1:24559/1‑150 (MQ=255)                                             gggCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAg                                                                                              >  1:79450/1‑150 (MQ=255)                                                    gCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGCACCTTAACACTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGaaaaa                                                                                       >  1:1005566/1‑150 (MQ=255)                                                    gCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGaaaaa                                                                                       >  2:364960/1‑150 (MQ=255)                                                     ccATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATTTAGAGGGCATCGTATGAGGTTTTGCTTCATTCTCTGACCGAATATTAACTCTGTACATATTATAGTTCAAAAACAGCAGCAATTAATAAATCTTATTTTTAGTGaataat                                                                                      >  2:983002/1‑150 (MQ=255)                                                     ccATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCaaa                                                                                                                             <  1:72648/111‑1 (MQ=255)                                                     ccATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCaaa                                                                                                                             >  2:72648/1‑111 (MQ=255)                                                        tttttCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGTATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCt                                                                                   <  1:255328/150‑1 (MQ=255)                                                                aaTGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttatttt                                                                           >  1:54992/1‑150 (MQ=255)                                                                aaTGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttatttt                                                                           >  2:769831/1‑150 (MQ=255)                                                                     ttCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGt                                                                      <  2:495574/150‑1 (MQ=255)                                                                        tAATATACTCGCAAAGATATCTAGAGGGCATCGTATGCGGTTTTGCTTCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCg                                                                   <  1:191773/150‑1 (MQ=255)                                                                        tAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCg                                                                   >  2:527994/1‑150 (MQ=255)                                                                            ataCTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGAt                                                               >  2:230309/1‑150 (MQ=255)                                                                              aCTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCGTTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGAttt                                                             <  2:1021192/150‑1 (MQ=255)                                                                                tCGCAAAGATATCTAGGGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGAttttt                                                           >  1:895395/1‑150 (MQ=255)                                                                                          atCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCCTTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCCGCATTTTTCCgactttat                                                 >  2:539252/1‑147 (MQ=255)                                                                                          atCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgacttgac                                                 <  2:206051/150‑1 (MQ=255)                                                                                                   gCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTc                                        >  1:323096/1‑150 (MQ=255)                                                                                                     aTCGTATGGGGTTTGGCTGCATTCCCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTg                                      <  2:228528/150‑1 (MQ=255)                                                                                                       cGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTggg                                    <  2:413879/150‑1 (MQ=255)                                                                                                         tATGAGGTTTTGCTGGATTCTCTGTGCGATGTTTAAGTGTGTGTATATTGTAGTTGAAAGTAGGTAGCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTTCGACTTGACTTTTAGTTCTGGGCt                                  >  2:1013239/1‑150 (MQ=255)                                                                                                           tGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTa                                >  1:801966/1‑150 (MQ=255)                                                                                                            gAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTaa                               <  1:923222/150‑1 (MQ=255)                                                                                                                ttttGCTGCATTCTCTGAGCTAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATCAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGt                           <  1:672970/150‑1 (MQ=255)                                                                                                                     cagaaTTCTGTGAGCGAATATTAACGCTGTTCATATTAAAGTTAAAAGCAAGAAGCATTGAATTAATGTTATTTTTATTTAAAATTTTTTTTATTTTGTCGTCCGGGATTTTTCCGCCTTGACTTTCAGTTCTGGGCTTAACTGTcaaag                      <  2:216978/146‑1 (MQ=255)                                                                                                                         aTTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCTACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagcaaa                  >  2:383052/1‑150 (MQ=255)                                                                                                                         aTTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTACTGGCCTTAACTGTcaaagcaaa                  >  2:141323/1‑150 (MQ=255)                                                                                                                           tctcTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGa                <  1:527994/150‑1 (MQ=255)                                                                                                                             tctGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTTGGCTTAACTGTCAAAGCAAAGata              >  2:870812/1‑150 (MQ=255)                                                                                                                                     aaTATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcact      >  1:721260/1‑150 (MQ=255)                                                                                                                                       tatTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactca    <  1:769831/150‑1 (MQ=255)                                                                                                                                         ttAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactcact  >  1:140395/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                               GATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACT  >  NC_000913_3_hsa_tpiA/1197359‑1197644

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