Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	950,582	G→A	R167C (CGC→TGC)	pflA ←	pyruvate formate‑lyase activating enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	950,582	0	G	A	100.0%	134.0 / NA	39	R167C (CGC→TGC)	pflA	pyruvate formate‑lyase activating enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/21);  total (18/21)

GCCACCCATTTGTGTTTGCCCAGCTCGTGGTAGGGGAGAAGCTCGATTTTCTCAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACT  >  NC_000913/950446‑950715                                                                                                                                         |                                                                                                                                      gCCACCCATTTGTGTTTGCCCAGCTCGTGGTAGGGGAGAAGCTCGATTTTCTCAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGAt                                                                                                                                    >  2:544779/1‑140 (MQ=255)      ccATTTGTGTTTGCCCAGCTCGTGGTAGGGGAGAAGCTCGATTTTCTCAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCcaca                                                                                                                               >  1:355826/1‑140 (MQ=255)                ttGCCCAGCTCGTGGTAGGGGAGAAGCTCGATTTTCTCAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACAtt                                                                                                                     >  2:28855/1‑140 (MQ=255)                ttGCCCAGCTCGTGGTAGGGGAGAAGCTCGATTTTCTCAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACAtt                                                                                                                     >  2:789627/1‑140 (MQ=255)                     cAGCTCGTGGTAGGGGAGAAGCTCGATTTTCTCAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGt                                                                                                                <  1:612783/140‑1 (MQ=255)                            tGGTAGGGGAGAAGCTCGATTTTCTCAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAg                                                                                                         <  2:76848/140‑1 (MQ=255)                              gTAGGGGAGAAGCTCGATTTTCTCAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGat                                                                                                       <  1:704872/140‑1 (MQ=255)                                  gggAGAAGCTCGATTTTCTCAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATAtt                                                                                                    <  1:779144/139‑1 (MQ=255)                                  gggAGAAGCTCGATTTTCTCAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATAtt                                                                                                    >  2:779144/1‑139 (MQ=255)                                   ggAGAAGCTCGATTTTCTCAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTa                                                                                                  >  2:39666/1‑140 (MQ=255)                                       aaGCTCGATTTTCTCAACGTTGCCCATATCACGGGTAAATTCACCGAGACGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGa                                                                                              >  2:762188/1‑140 (MQ=255)                                        aGCTCGATTTTCTCAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGaa                                                                                             >  2:217849/1‑140 (MQ=255)                                          cTCGATTTTCTCAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACt                                                                                           <  1:306405/140‑1 (MQ=255)                                                tttCTCAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCg                                                                                     <  1:428971/140‑1 (MQ=255)                                                  tctcAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTg                                                                                   <  1:28855/140‑1 (MQ=255)                                                  tctcAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTg                                                                                   <  1:217849/140‑1 (MQ=255)                                                   ctcAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGc                                                                                  <  2:4844/140‑1 (MQ=255)                                                    tcAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCg                                                                                 <  2:95244/140‑1 (MQ=255)                                                            gCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTgg                                                                         <  1:789627/140‑1 (MQ=255)                                                             cccATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGa                                                                        >  1:519234/1‑140 (MQ=255)                                                                     aCGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCaa                                                                <  1:39666/140‑1 (MQ=255)                                                                       gggTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAAcc                                                              >  1:546268/1‑140 (MQ=255)                                                                            aaTTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGAtt                                                         <  2:519234/140‑1 (MQ=255)                                                                                   cGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTtggtgg                                                  >  1:202909/1‑140 (MQ=255)                                                                                        cGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTtggt                                                    >  1:74755/1‑133 (MQ=255)                                                                                        cGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTtggt                                                    <  2:74755/133‑1 (MQ=255)                                                                                         gATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCtcg                                            <  2:546268/140‑1 (MQ=255)                                                                                               cTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTc                                      <  1:544779/140‑1 (MQ=255)                                                                                                 gAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCAt                                    >  2:826702/1‑140 (MQ=255)                                                                                                       tcatcgtcaGACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGttt                              >  2:157818/1‑140 (MQ=255)                                                                                                            gtcaGTCCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGAt                         <  2:5474/140‑1 (MQ=255)                                                                                                                gACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAg                     <  2:103679/140‑1 (MQ=255)                                                                                                                gACCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAg                     <  1:324076/140‑1 (MQ=255)                                                                                                                 aCCAGCCTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGc                    >  1:518222/1‑140 (MQ=255)                                                                                                                      ccTGGGACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTAc               >  1:332187/1‑140 (MQ=255)                                                                                                                         gggACAACAACGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAg            <  2:518222/140‑1 (MQ=255)                                                                                                                                caacGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGtt     >  2:405514/1‑140 (MQ=255)                                                                                                                                caacGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGA‑TTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTa    >  2:70470/1‑140 (MQ=255)                                                                                                                                   cGTAGCAGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACt  <  1:157818/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      GCCACCCATTTGTGTTTGCCCAGCTCGTGGTAGGGGAGAAGCTCGATTTTCTCAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACT  >  NC_000913/950446‑950715

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