Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,547,856	C→A	G152G (GGC→GGA)	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,547,856	0	C	A	100.0%	88.4 / NA	26	G152G (GGC→GGA)	fdnG	formate dehydrogenase N subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (15/11);  total (15/11)

ACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATGCTGGCGGTAGACAACCAGGCGCG  >  NC_000913/1547729‑1547981                                                                                                                                |                                                                                                                              aCCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGAGTAACGGTAAAc                                                                                                                   <  2:440118/140‑1 (MQ=255)            cGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTc                                                                                                        >  1:493231/1‑140 (MQ=255)            cGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCGGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTc                                                                                                        <  2:69325/140‑1 (MQ=255)             gTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCt                                                                                                       >  1:137190/1‑140 (MQ=255)              tGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTa                                                                                                      <  1:253480/140‑1 (MQ=255)              tGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTa                                                                                                      >  2:378108/1‑140 (MQ=255)                 gcCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTACCg                                                                                                   >  1:258906/1‑140 (MQ=255)                 gcCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTACCg                                                                                                   >  2:349451/1‑140 (MQ=255)                      gTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTACCGGTATg                                                                                              <  2:258906/140‑1 (MQ=255)                                             aGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCaa                                                                       <  2:493231/140‑1 (MQ=255)                                              gCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAAc                                                                      >  1:429/1‑140 (MQ=255)                                                   gaagaaGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAAc                                                                 >  2:76930/1‑140 (MQ=255)                                                                         cGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAAttt                                           >  1:457973/1‑140 (MQ=255)                                                                         cGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAAttt                                           >  2:430627/1‑140 (MQ=255)                                                                           aaGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGc                                         >  1:432921/1‑140 (MQ=255)                                                                            aGCTGATGAAAGCTGACCGTGACACTAACTTTATTGAAAAGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGcc                                        <  1:544756/140‑1 (MQ=255)                                                                              cTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCg                                      <  2:457973/140‑1 (MQ=255)                                                                              cTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCg                                      >  2:536912/1‑140 (MQ=255)                                                                                   gaAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTccc                                 >  1:437913/1‑140 (MQ=255)                                                                                               tGACGCTAACTTTATTGAAAAGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATGCTgg                     <  1:336237/140‑1 (MQ=255)                                                                                                        cTTTATTGAAAAGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATGCTGGCGGTAGACa            >  1:468672/1‑140 (MQ=255)                                                                                                              tGAAAAGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATGCTGGCGGTAGACAACCAgg      <  2:137190/140‑1 (MQ=255)                                                                                                                aaaaGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATGCTGGCGGTAGACAACCAGgcg    <  1:349451/140‑1 (MQ=255)                                                                                                                 aaaGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATGCTGGCg                   <  2:398212/124‑1 (MQ=255)                                                                                                                 aaaGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATGCTGGCg                   >  1:398212/1‑124 (MQ=255)                                                                                                                  aaGAACGAGCAGGGAGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATGCTGGCGGTAGACAACCAGgcgcg  >  2:121598/1‑140 (MQ=255)                                                                                                                                |                                                                                                                              ACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATGCTGGCGGTAGACAACCAGGCGCG  >  NC_000913/1547729‑1547981

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