Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,144,522	G→A	M176I (ATG→ATA)	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,522	0	G	A	100.0%	134.4 / NA	37	M176I (ATG→ATA)	pflD	putative formate acetyltransferase 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/19);  total (18/19)

ATAAAGGCCAGGGGCACATTATTATTGATTACCCACGCCTGTTGAATCACGGGCTGGGTGAGCTGGTGGCACAGATGCAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCCGCGTTGTTACTGCTGGAAGCCTCGCAGAAAC  >  NC_000913/4144446‑4144598                                                                             |                                                                             aTAAAGGCCAGGGGCACATTATTATTGATTACCCACGCCTGTTGAATCACGGGCTGGGTGAGCTGGTGGCACAGATacagc                                                                          <  1:820583/81‑1 (MQ=255)    aaGGCCAGGGGCACATTATTATTGATTACCCACGCCTGTTGAATCACGGGCTGGGTGAGCTGGTGGCACAGATACAGCaa                                                                        <  2:820584/80‑1 (MQ=255)    aaGGCCAGGGGCACATTATTATTGATTACCCACGCCTGTTGAATCACGGGCTGGGTGAGCTGGTGGCACAGATACAGCaa                                                                        >  2:820585/1‑80 (MQ=255)      ggCCAGGGGCACATTATTATTGATTACCCACGCCTGTTGAATCACGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACAt                                                                     >  2:820623/1‑81 (MQ=255)         cAGGGGCACATTATTATTGATTACCCACGCCTGTTGAATCACGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGt                                                                  <  1:820586/81‑1 (MQ=255)           ggggCACATTATTATTGATTACCCACGCCTGTTGAATCACGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGt                                                                  >  2:820614/1‑79 (MQ=255)                acaTTATTATTGATTACCCACGCCTGTTGAATCACGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCaa                                                            >  2:820617/1‑80 (MQ=255)                   ttattattGATTACCCACGCCTGTTGAATCACGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGc                                                        >  1:820627/1‑81 (MQ=255)                            aTTACCCACGCCTGTTGAATCACGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAAt                                                <  2:820587/80‑1 (MQ=255)                                cccACGCCTGTTGAATCGCGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACt                                            <  1:820588/80‑1 (MQ=255)                                cccACGCCTGTTGAATCACGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACtt                                           >  2:820631/1‑81 (MQ=255)                                   aCGCCTGTTGAATCACGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTa                                        <  2:820589/81‑1 (MQ=255)                                    cGCCTGTTGAATCACGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTAt                                       <  2:820590/81‑1 (MQ=255)                                       cTGTTGAATCACGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAg                                    >  1:820628/1‑81 (MQ=255)                                         gTTGAATCACGGGCTGGGTGAGCTGGCGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGc                                  >  2:820621/1‑81 (MQ=255)                                           tGAATCACGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAACCACTTTTATCAGGCcg                                >  2:820645/1‑81 (MQ=255)                                            gAATCACGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCcgc                               >  1:820638/1‑81 (MQ=255)                                            gAATCACGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCcg                                <  1:820591/80‑1 (MQ=255)                                             aaTCACGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCcgcg                              >  2:820641/1‑81 (MQ=255)                                             aaTCACGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCcgcg                              <  1:820592/81‑1 (MQ=255)                                              aTCACGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCATCAGCCGGAGAATCACTTTTATCAGGCCGCgt                             >  2:820644/1‑81 (MQ=255)                                              aTCACGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCCGCgt                             <  2:820594/81‑1 (MQ=255)                                              aTCACGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCCGCgt                             <  1:820593/81‑1 (MQ=255)                                                cACGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCCGCgttg                           >  2:820624/1‑81 (MQ=255)                                                 aCGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCCGCgttgt                          <  2:820595/81‑1 (MQ=255)                                                 aCGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCCGCgttgt                          >  1:820634/1‑81 (MQ=255)                                                  cGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCCGCgttgtt                         <  2:820596/81‑1 (MQ=255)                                                  cGGGCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCCGCgttgtt                         <  2:820597/81‑1 (MQ=255)                                                    ggCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCCGCGTTGTTAc                       <  1:820598/81‑1 (MQ=255)                                                     gCTGGGTGAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCCGCGTTGTTAct                      <  2:820599/81‑1 (MQ=255)                                                            gAGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCCGCGTTGTTACTGCTgg                 >  2:820630/1‑79 (MQ=255)                                                             aGCTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCCGCGTTGTTACTGCTGGAAg              <  1:820600/81‑1 (MQ=255)                                                               cTGGTGGCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCCGCGTTGTTACTGCTGGAAGcc            <  2:820601/81‑1 (MQ=255)                                                                 ggtggCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCCGCGTTGTTACTGCTGGAAGCCTc          >  2:820635/1‑81 (MQ=255)                                                                     gCACAGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCCGCGTTGTTACTGCTGGAAGCCTCGCa       <  2:820602/80‑1 (MQ=255)                                                                       acaGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCCGCGTTGTTACTGCTGGAAGCCTCGCAGa     >  1:820629/1‑80 (MQ=255)                                                                       acaGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCCGCGTTGTTACTGCTGGAAGCCTCGCAGa     >  1:820652/1‑80 (MQ=255)                                                                         aGATACAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCCGCGTTGTTACTGCTGGAAGCCTCGCAGAAac  <  2:820603/81‑1 (MQ=255)                                                                             |                                                                             ATAAAGGCCAGGGGCACATTATTATTGATTACCCACGCCTGTTGAATCACGGGCTGGGTGAGCTGGTGGCACAGATGCAGCAACATTGTCAGCAACAGCCGGAGAATCACTTTTATCAGGCCGCGTTGTTACTGCTGGAAGCCTCGCAGAAAC  >  NC_000913/4144446‑4144598

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