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

ACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACT  >  NC_000913_3_hsa_tpiA/1197356‑1197640                                                                                                                                                   |                                                                                                                                           aCTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGTGTGGACATTTTTTACATTGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                                                                                                          <  1:226369/149‑1 (MQ=255)  cTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCtgtg                                                                                                                                         <  2:22367/149‑1 (MQ=255)          ttCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCCAATGAGGTTTTGCTGCATTCTCTGACCGAATATTAACTCTTTGCAtattat                                                                                                                                 >  2:603635/1‑149 (MQ=255)                tgAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCa                                                                                                                           <  2:802273/148‑1 (MQ=255)                            ttttCCGCAGCCAAAAACAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGAGATCTAGAGGGCATCGTATAAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTATGTGCATATTATAGTTCAAAGCATGTACCa                                                                                                               >  1:837552/1‑149 (MQ=255)                                      ccAAAATCAAGGGCGGGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt                                                                                                     <  2:250497/149‑1 (MQ=255)                                       cAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATg                                                                                                    >  1:624806/1‑149 (MQ=255)                                        aaaaTCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCtgtg                                                                                                                                         <  1:435281/111‑1 (MQ=255)                                        aaaaTCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCtgtg                                                                                                                                         >  2:435281/1‑111 (MQ=255)                                                             tttCTCAATGATTCTAATAGACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTAATATTAATAAAAAATTCtt                                                                              <  2:269136/149‑1 (MQ=255)                                                              ttCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttt                                                                             <  2:829552/149‑1 (MQ=255)                                                               tctcAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCtttt                                                                            <  1:558250/149‑1 (MQ=255)                                                                   aaTGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttattt                                                                        >  2:403585/1‑149 (MQ=255)                                                                           tAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTcc                                                                >  2:825225/1‑149 (MQ=255)                                                                                   tCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGAtttt                                                        >  2:305721/1‑149 (MQ=255)                                                                                     gCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTc                                                      <  2:353585/149‑1 (MQ=255)                                                                                      cAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTACGAGATTTTTcc                                                     >  1:300398/1‑149 (MQ=255)                                                                                         aGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgac                                                  <  2:624806/149‑1 (MQ=255)                                                                                           atatCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgactt                                                >  1:779017/1‑149 (MQ=255)                                                                                              tCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgacttgac                                             >  2:802526/1‑149 (MQ=255)                                                                                               cTAGAGGGCATCGTATGAGGTTTTGCTGCATCCTCTGAGCGAATATTAACTCTGTGCATATtata                                                                                                                                <  2:579767/65‑1 (MQ=255)                                                                                               cTAGAGGGCATCGTATGAGGTTTTGCTGCATCCTCTGAGCGAATATTAACTCTGTGCATATtata                                                                                                                                >  1:579767/1‑65 (MQ=255)                                                                                                tAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgacttgact                                            >  1:803061/1‑148 (MQ=255)                                                                                                tAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgacttgact                                            <  2:803061/148‑1 (MQ=255)                                                                                                  gaaggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTc                                         <  2:531383/146‑1 (MQ=255)                                                                                                              tGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCtt                             >  2:130382/1‑149 (MQ=255)                                                                                                              tGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCtt                             >  1:507529/1‑149 (MQ=255)                                                                                                               gAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTa                            >  1:368587/1‑149 (MQ=255)                                                                                                                   gtttGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAGAGCAAGTACCATTCAATTAATGTTATTTATAGTCAAAAATTCTATTATTTTATCGTTCGGGATTCTTCCAACTTGACTTTCATTGTTGCACTTAACTg                        >  2:730726/2‑149 (MQ=255)                                                                                                                         tgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTTTTTTTTTTTTTTCAGTTCTGAGCTTAACTGTtaaag                  >  2:717413/1‑149 (MQ=255)                                                                                                                         tgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGCCTTGACTTTCAGTTCTGGCCTTATCTGTcaacg                  >  2:587707/1‑147 (MQ=255)                                                                                                                             ttCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagcaaa              <  1:802526/149‑1 (MQ=255)                                                                                                                               ctctGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGa            >  1:310942/1‑149 (MQ=255)                                                                                                                                    aGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGataaaa       >  2:508144/1‑147 (MQ=255)                                                                                                                                     gCGAATATTAACTCTGTGCATATTACAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCTACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATAc      <  1:867293/149‑1 (MQ=255)                                                                                                                                        gatatTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGTCTTAACTGTCAAAGAAAAGATATTACcat   >  2:392397/2‑148 (MQ=255)                                                                                                                                         atatTAACTCTGTGCATATTATAGTTCAAAGCAAGTACTATTTAATTAATGTTATTTTTAGTGAAAAATTATTTTATTTTATAGTTTGGGTTTTTTTTGAATTGAATTTTAGTTGTGGGCTTAAATGTTAAAGTAAAGataatagtagt  >  2:897314/1‑144 (MQ=255)                                                                                                                                                   |                                                                                                                                           ACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACT  >  NC_000913_3_hsa_tpiA/1197356‑1197640

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