Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,225,294	A→G	G489G (GGA→GGG)	metH →	homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,225,294	0	A	G	100.0%	122.1 / NA	35	G489G (GGA→GGG)	metH	homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/17);  total (18/17)

TATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAG  >  NC_000913/4225168‑4225427                                                                                                                               |                                                                                                                                      tATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACt                                                                                                                           >  1:465767/1‑139 (MQ=255)         gCAAAGGCATTGTTAACGCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACg                                                                                                                   <  2:522150/139‑1 (MQ=255)           aaaGGCATTGTTAACTCTATCTTGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGa                                                                                                                 >  2:1043751/1‑139 (MQ=255)                     ttAACTCTATCTTGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGAt                                                                                                       <  1:1043751/139‑1 (MQ=255)                             aTCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGt                                                                                                >  1:1059052/1‑138 (MQ=255)                             aTCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGt                                                                                                <  2:1059052/138‑1 (MQ=255)                              tctcGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCg                                                                                              >  2:932393/1‑139 (MQ=255)                               ctcGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCgg                                                                                             <  1:106453/139‑1 (MQ=255)                                   aTGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGt                                                                                         >  1:404862/1‑139 (MQ=255)                                   aTGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGt                                                                                         >  2:1020160/1‑139 (MQ=255)                                      aaaGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACa                                                                                      <  1:932393/139‑1 (MQ=255)                                             gCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCt                                                                               <  1:1020160/139‑1 (MQ=255)                                              cGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTc                                                                              >  2:946741/1‑139 (MQ=255)                                              cGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTc                                                                              >  2:484395/1‑139 (MQ=255)                                                tCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAg                                                                                                                                    <  1:876926/83‑1 (MQ=255)                                                tCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAg                                                                                                                                    >  2:876926/1‑83 (MQ=255)                                                 cGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCAcc                                                                           >  1:545869/1‑139 (MQ=255)                                                           aTCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTg                                                                 >  2:756809/1‑139 (MQ=255)                                                                 cACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGc                                                                                           >  1:916102/1‑107 (MQ=255)                                                                 cACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGc                                                                                           <  2:916102/107‑1 (MQ=255)                                                                 cACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTcc                                                           <  2:404862/139‑1 (MQ=255)                                                                 cACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTcc                                                           <  1:946741/139‑1 (MQ=255)                                                                  aCGCGAAATTGTTGTGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCgg                                                                                             >  2:399194/1‑104 (MQ=255)                                                                  aCGCGAAATTGTTGTGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCgg                                                                                             <  1:399194/104‑1 (MQ=255)                                                                        aaTTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCg                                                         >  1:590321/1‑134 (MQ=255)                                                                        aaTTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCg                                                         <  2:590321/134‑1 (MQ=255)                                                                           tgttgCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGTTTCCCGCCagaaga                                                 >  1:464512/1‑139 (MQ=255)                                                                                                cAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCagca                                                   <  2:645462/116‑3 (MQ=255)                                                                                                cAGTGGTGGTAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCagaa                                                   >  1:645462/1‑116 (MQ=255)                                                                                                     gtggtAATGGCCTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTcg                       <  1:756809/139‑1 (MQ=255)                                                                                                              ggcTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCaa              >  2:644317/3‑139 (MQ=255)                                                                                                               ccTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAAc             <  1:484395/139‑1 (MQ=255)                                                                                                               ccTTTGACGAACAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAAc             >  2:61920/1‑139 (MQ=255)                                                                                                                          cAGGGGCAGGCCGATACTCGCGCACGGATAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAag  <  1:644317/139‑1 (MQ=255)                                                                                                                          cAGGGGCAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGTTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAag  <  2:464512/139‑1 (MQ=255)                                                                                                                               |                                                                                                                                      TATCCAGGGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAG  >  NC_000913/4225168‑4225427

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑