Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,047,030	A→G	V671A (GTG→GCG)	gcvP ←	glycine decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,047,030	0	A	G	100.0%	97.6 / NA	29	V671A (GTG→GCG)	gcvP	glycine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/14);  total (15/14)

ATGCAGAAAGTTTTATGTAGGTTAAGGTGTGAAACGTCCGCACCAATAAAGCCCGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTC  >  NC_000913/3046907‑3047146                                                                                                                            |                                                                                                                     aTGCAGAAAGTTTTATGTAGGTTAAGGTGTGAAACGTCCGCACCAATAAAGCCCGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCGCGACTTcac                                                                                                             <  1:101987/133‑1 (MQ=255)  ttcAGAAAGCTTTATCTAGGTTAAGGGGTGAACCGTCCGCAACAATAAAGCCCGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCGCGACTTcaca                                                                                                            <  2:245380/131‑1 (MQ=255)  tGCAGAAAGTTTTATGTAGGTTAAGGTGTGAAACGTCGGCACCAATAAAGCCCGGCGAGGTAACGCCAACCTGGGGGTTCATGTTCGCGCCATCAAGGTAAACCTAACCGCCGAACTGATGCGCGACTTcaca                                                                                                            <  2:328448/133‑1 (MQ=255)  tGCAGAAAGTTTTATATAGGTTAACGTGCGAATCGTCCGCACCAATAAAGCCCGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCGCGACTTcaca                                                                                                            <  2:199909/133‑1 (MQ=255)              tATGTAGGTTAAGGTGTGAAACGTCCGCACCAATAAAGCCCGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCGCGACTTCACACACTTCACGGAt                                                                                                >  1:444264/1‑133 (MQ=255)                          ggTGTGAAACGTCCGCACCAATAAAGCCCGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCGCGACTTCACACACTTCACGGATCGTTTCTTCATa                                                                                    >  1:271154/1‑133 (MQ=255)                          ggTGTGAAACGTCCGCACCAATAAAGCCCGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCGCGACTTCACACACTTCACGGATCGTTTCTTCATa                                                                                    >  2:2863/1‑133 (MQ=255)                              tgAAACGTCCGCACCAATAAAGCCCGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCGCGACTTCACACACTTCACGGATCGTTTCTTCATACACg                                                                                >  1:258212/1‑133 (MQ=255)                               gAAACGTCCGCACCAATAAAGCCCGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATTCGCGACTTTACAGACTACACGGATCATTTCTTCATACACGc                                                                               >  2:161397/1‑133 (MQ=255)                               gAAACGTCCGCACCAATAAAGCCCGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCGCGACTTCAGACACTTCGCGGATCATTTCTTTATACACGc                                                                               >  1:226541/1‑133 (MQ=255)                               gAAACGTCCGCACCAATAAAGCCCGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCGCGACTTCAAACACTTCACGGTTCCTTTCTTCATACACGc                                                                               >  1:179454/1‑133 (MQ=255)                                aaaCGTCCGCACCAATAAAGCCCGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCGCGACTTCACACACTTCACGGATCGTTTCTTCATACACGcc                                                                              <  1:321090/133‑1 (MQ=255)                                  aCGTCCGCCCAAATAAAGCCCGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCGCGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGt                                                                            <  1:313284/133‑1 (MQ=255)                                  aCGTCCGCACAAATAGAGCCCGGGGAGGTGATCCCAACCTGTGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCGCGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGt                                                                            <  2:226541/133‑1 (MQ=255)                                             aaTAAAGCCCGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCGCGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGa                                                                 <  1:323726/133‑1 (MQ=255)                                                  aGCCCGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCGCGACTTCACACACTTCACGGATTGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGt                                                            <  1:258597/133‑1 (MQ=255)                                                  aGCCCGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCGCGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGt                                                            <  1:81698/133‑1 (MQ=255)                                                     ccGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCGCGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCAc                                                         <  2:258212/133‑1 (MQ=255)                                                      cGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCGCGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCAcc                                                        >  1:340782/1‑133 (MQ=255)                                                      cGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCGCGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCAcc                                                        >  1:473417/1‑133 (MQ=255)                                                       ggCGAGGTAATGCCACCCAGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCGCGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCa                                                       <  1:161397/133‑1 (MQ=255)                                                        gCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCGCGACTTCACACACTTCACGGATCGTTTCTTCATGCACGCCGTGGGGAGTAGGATAGGTCACCAt                                                      >  1:155066/1‑133 (MQ=255)                                                              tGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCGCGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACa                                                >  1:366027/1‑133 (MQ=255)                                                                                                  ggTAAACCTGACCGCCGAACTGATGCGCGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTcgcgcg            <  1:297880/133‑1 (MQ=255)                                                                                                      acccTGACCGCCGAACTGATGCGCGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGa        <  2:366027/131‑1 (MQ=255)                                                                                                      aaCCTGACCGCCGAACTGATGCGCGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGa        >  2:398314/1‑133 (MQ=255)                                                                                                            aCCGCCGAACTGATGCGCGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGAtcagtc  >  1:81248/1‑133 (MQ=255)                                                                                                            aCCGCCGAACTGATGCGCGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGAtcagtc  >  2:309228/1‑133 (MQ=255)                                                                                                            aCCGCCGAACTGATGCGCGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGAtcagtc  >  1:287566/1‑133 (MQ=255)                                                                                                                            |                                                                                                                     ATGCAGAAAGTTTTATGTAGGTTAAGGTGTGAAACGTCCGCACCAATAAAGCCCGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTC  >  NC_000913/3046907‑3047146

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: ‑