Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,548,895	A→G	S499G (AGC→GGC)	fdnG →	formate dehydrogenase N subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,548,895	0	A	G	100.0%	48.5 / NA	16	S499G (AGC→GGC)	fdnG	formate dehydrogenase N subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (8/8);  total (8/8)

ACTTAGGCCTGCTCTCTACCAGCCTGCCAGGTTATCTGACGCTGCCGTCAGAAAAACAGGTTGATTTGCAGTCGTATCTGGAAGCGAACACGCCGAAAGCGACGCTGGCTGATCAGGTGAACTACTGGAGCAACTATCCGAAGTTCTTCGTTAGCCTGATGAAATCTTTCTATGGCGATGCCGCGCAGAAAGAGAACAACTGGGGCTATGACTGGCTGCCGAAGTGGGACCAGACCTACGACGTCATCAAGTATTTCA  >  NC_000913/1548767‑1549024                                                                                                                                 |                                                                                                                                  actTAGGCCTGCTCTCTACCAGCCTGCCAGGTTATCTGACGCTGCCGTCAGAAAAACAGGTTGATTTGCAGTCGTATCTGGAAGCGAACACGCCGAAAGCGACGCTGGCTGATCAGGTGAACTACTGGGGCaa                                                                                                                               >  2:505988/1‑133 (MQ=255)            ctctctGCCAGCCTGCCAGGTTATCTGACGCTGCCGTCAGAAAAACAGGTTGATTTGCAGTCGTATCTGGAAGCGAACACGCCGAAAGCGACGCTGGCTGATCAGGTGAACTACTGGGGCAACt                                                                                                                             <  1:57011/124‑1 (MQ=255)            ctctctACCAGCCTGCCAGGTTATCTGACGCTGCCGTCAGAAAAACAGGTTGATTTGCAGTCGTATCTGGAAGCGAACACGCCGAAAGCGACGCTGGCTGATCAGGTGAACTACTGGGGCAACt                                                                                                                             >  2:57011/1‑124 (MQ=255)                 tACCAGCCTGCCAGGTTATCTGACGCTGCCGTCAGAAAAACAGGTTGATTTGCAGTCGTATCTGGAAGCGAACACGCCGAAAGCGACGCTGGCTGATCAGGTGAACTACTGGGGCAACTATCCGAAGttcttc                                                                                                               >  2:94067/1‑133 (MQ=255)                  aCCAGCCTGCCAGGTTATCTGACGCTGCCGTCAGAAAAACAGGTTGATTTGCAGTCGTATCTGGAAGCGAACACGCCGAAAGCGACGCTGGCTGATCAGGTGAACTACTGGGGCAACTATCCGAAGTTCTTCg                                                                                                              <  1:322138/133‑1 (MQ=255)                       ccTGCCAGGTTATCTGACGCTGCCGTCAGAAAAACAGGTTGATTTGCAGTCGTATCTGGAAGCGAACACGCCGAAAGCGACGCTGGCTGATCAGGTGAACTACTGGGGCAACTATCCGAAGTTCTTCGTTAGc                                                                                                         <  1:18071/133‑1 (MQ=255)                         tGCCAGGTTATCTGACGCTGCCGTCAGAAAAACAGGTTGATTTGCAGTCGTATCTGGAAGCGAACACGCCGAAAGCGACGCTGGCTGATCAGGTGAACTACTGGGGCAACTATCCGAAGTTCTTCGTTAGCCt                                                                                                       >  2:119805/1‑133 (MQ=255)                           ccAGGTTATGTGACGCTGCCGTCAGAAAAACAGGTTGATTTGCAGTCGTATCTGGAAGAGAACACGCCGAAAGCGACGCTGGCTGATCAGGTGAACTACTGGGGCAACTATCCGAAGATCTTCGTTAGCCtga                                                                                                     >  2:317088/1‑133 (MQ=255)                                          cTGCCGTCAGAAAAACAGGTTGATTTGCAGTCGTATCTGGAAGCGAACACGCCGAAAGCGACGCTGGCTGATCAGGTGAACTACTGGGGCAACTATCCGAAGTTCTTCGTTAGCCTGATGAAATCTTTCTATg                                                                                      >  2:340070/1‑133 (MQ=255)                                                      aaaCAGGTTGATTTGCAGTCGTATCTGGAAGCGAACACGCCGAAAGCGACGCTGGCTGATCAGGTGAACTACTGGGGCAACTATCCGAAGTTCTTCGTTAGCCTGATGAAATCTTTCTATGGCGATGCCgcgc                                                                          <  1:390151/133‑1 (MQ=255)                                                      aaaCAGGTTGATTTGCAGTCGTATCTGGAAGCGAACACGCCGAAAGCGACGCTGGCTGATCAGGTGAACTACTGGGGCAACTATCCGAAGTTCTTCGTTAGCCTGATGAAATCTTTCTATGGCGATGCCgcgc                                                                          <  1:358764/133‑1 (MQ=255)                                                                                  aaGCGAACACGCCGAAAGCGACGCTGGCTGATCAGGTGAACTACTGGGGCAACTATCCGAAGTTCTTCGTTAGCCTGATGAAATCTTTCTATGGCGATGCCGCGCAGAAAGAGAACAACTGGGGCTATGACTg                                              >  2:46897/1‑133 (MQ=255)                                                                                      gAACACGCCGAAAGCGACGCTGGCTGATCAGGTGAACTACTGGGGCAACTATCCGAAGTTCTTCGTTAGCCTGATGAAATCTTTCTATGGCGATGCCGCGCAGAAAGAGAACAACTGGGGCTATGACTGGCTg                                          >  2:467402/1‑133 (MQ=255)                                                                                       aaCACGCCGAAAGCGACGCTGGCTGATCAGGTGAACTACTGGGGCAACTATCCGAAGTTCTTCGTTAGCCTGATGAAATCTTTCTATGGCGTTGCCGCGCAGAAAGAGTACAACTGGAGCTACGACTGGCTGc                                         >  2:205305/1‑133 (MQ=255)                                                                                        acacGCCGAAAGCGACGCTGGCTGATCAGGTGAACTACTGGGGCAACTATCCGAAGTTCTTCGTTAGCCTGATGAAATCTTTCTATGGCGATGCCGCGCAGATAGAGAACGACTGGGTCTATGAGTGGCTGcc                                        >  1:181793/1‑133 (MQ=255)                                                                                              cGAAAGCGACGCTGGCTGATCAGGTGAACTACTGGGGCAACTATCCGAAGTTCTTCGTTAGCCTGATGAAATCTTTCTATGGCGATGCCGCGCAGAAAGAGAACAACTGGGGCTATGACTGGCTGCCGAAGTg                                  >  2:284789/1‑133 (MQ=255)                                                                                                            gctgATCAGGTGAACTACTGGGGCAACTATCCGAAGTTCTTCGTTAGCCTGATGAAATCTTTCTATGGCGATGCCGCGCAGAAAGAGAACAACTGGGGCTATGACTGGCTGCCGAAGTGGGACCAGACCTacg                    <  1:441404/133‑1 (MQ=255)                                                                                                               gATCAGGTAAACTACTGGGGCACCTATCCGGAGTTCTTCGTTAGCCTGATGAAATCTTTCTATGGCGATGCCGCGCAGAAAGAGAACAACTGGGGCTATGACTGGCTGCCGAAGTGGGACCAGACCTacgacg                 <  1:238293/133‑1 (MQ=255)                                                                                                                              tGGGGCAACTATCCGAAGTTCTTCGTTAGCCTGATGAAATCTTTCTATGGCGATGCCGCGCAGAAAGAGAACAACTGGGGCTATGACTGGCTGCCGAAGTGGGACCAGACCTACGACGTCATCAAGTATTTCa  <  1:317088/133‑1 (MQ=255)                                                                                                                                 |                                                                                                                                  ACTTAGGCCTGCTCTCTACCAGCCTGCCAGGTTATCTGACGCTGCCGTCAGAAAAACAGGTTGATTTGCAGTCGTATCTGGAAGCGAACACGCCGAAAGCGACGCTGGCTGATCAGGTGAACTACTGGAGCAACTATCCGAAGTTCTTCGTTAGCCTGATGAAATCTTTCTATGGCGATGCCGCGCAGAAAGAGAACAACTGGGGCTATGACTGGCTGCCGAAGTGGGACCAGACCTACGACGTCATCAAGTATTTCA  >  NC_000913/1548767‑1549024

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