Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,144,838	A→G	T282A (ACA→GCA)	pflD →	putative formate acetyltransferase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,144,838	0	A	G	100.0%	60.2 / NA	19	T282A (ACA→GCA)	pflD	putative formate acetyltransferase 2
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/10);  total (9/10)

AAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGCTCGAATCTTTATGGGTGAAATGCAACGACATCGTGCTGTTGCGCTCCACCAGTAGCGCGCGTTATTTCGCCGGTTTCCCG  >  NC_000913/4144712‑4144963                                                                                                                               |                                                                                                                              aaGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGGcatcat                                                                                                                         <  1:364427/133‑1 (MQ=255)                 ggCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCTGGCATCATTAACCCAGGGCGAAGa                                                                                                         <  1:417644/133‑1 (MQ=255)                     ggCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGGCATCATTAACCCAGGGCGAAGATgcg                                                                                                     <  2:491065/133‑1 (MQ=255)                     ggCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGGCATCATTAACCCAGGGCGAAGATgcg                                                                                                     <  1:411936/133‑1 (MQ=255)                                    cTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGGCATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGaa                                                                                      <  2:66491/133‑1 (MQ=255)                                                 atcatTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGGCATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGCTCGAATCtt                                                                         >  1:435104/1‑133 (MQ=255)                                                      tCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGGCATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGCTCGAATCTTTATgg                                                                    <  2:435104/133‑1 (MQ=255)                                                             tACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTACCAGGCATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGCTCGAATCTTTATGGGTg                                                                 <  1:303276/129‑1 (MQ=255)                                                             tACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTACCAGGCATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGCTCGAATCTTTATGGGTg                                                                 >  2:303276/1‑129 (MQ=255)                                                               cGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGGCATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGCTCGAATCTTTATGGGTGAAATGc                                                           >  2:460203/1‑133 (MQ=255)                                                                  aTCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGGCATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGCTCGAATCTTTATGGGTGAAATGCaa                                                         >  1:350110/1‑132 (MQ=255)                                                                  aTCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGGCATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGCTCGAATCTTTATGGGTGAAATGCaa                                                         <  2:350110/132‑1 (MQ=255)                                                                         gCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGGCATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGCTCGAATCTTTATGGGTGAAATCCAACGACATCg                                                 >  2:346548/1‑133 (MQ=255)                                                                                    aTCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGGCATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGCTCGAATCTTTATGGGTGAAATGCAACGACATCGTGCTGTTgcgc                                      <  1:500594/133‑1 (MQ=255)                                                                                                                ttGCCGTTCTATCAGGCATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGCTCGAATCTTTATGGGTGAAATGCAACGACATCGTGCTGTTGCGCTCCACCAGTAGCGCGCGTTATTTCGCCg          >  1:362837/1‑133 (MQ=255)                                                                                                                   ccGTTCTATCAGGCATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGCTCGAATCTTTATGGGTGAAATGCAACGACATCGTGCTGTTGCGCTCCACCAGTAGCGCGCGTTATTTCGCCGGtt       >  1:102841/1‑133 (MQ=255)                                                                                                                    cGTTCTATCAGGCATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGCTCGAATCTTTATGGGTGAAATGCAACGACATCGTGCTGTTGCGCTCCACCAGTAGCGCGCGTTATTTCGCCGGttt      >  2:378666/1‑133 (MQ=255)                                                                                                                    cGTTCTATCAGGCATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGCTCGAATCTTTATGGGTGAAATGCAACGACATCGTGCTGTTGCGCTCCACCAGTAGCGAGCATTATTTAGCCGGttt      >  1:242659/1‑133 (MQ=255)                                                                                                                        cTATCAGGCATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGCTCGAATCTTTATGGGTGAAATGCAACGACATCGTGCTGTTGCGCTCCACCAGTAGCGCGCGTTATTTCGCCGGTTTCCCg  <  2:102841/133‑1 (MQ=255)                                                                                                                               |                                                                                                                              AAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAACTGCTCGAATCTTTATGGGTGAAATGCAACGACATCGTGCTGTTGCGCTCCACCAGTAGCGCGCGTTATTTCGCCGGTTTCCCG  >  NC_000913/4144712‑4144963

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