Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	941,442	T→G	S162A (TCT→GCT)	dmsA →	dimethyl sulfoxide reductase subunit A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	941,442	0	T	G	100.0%	104.5 / NA	30	S162A (TCT→GCT)	dmsA	dimethyl sulfoxide reductase subunit A
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (16/14);  total (16/14)

GCATTAGCTGGGAAGAAGCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGTCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTC  >  NC_000913/941368‑941519                                                                           |                                                                              gcATTAGCTGGGAAGAAGCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGgctatc                                                                          >  1:176696/1‑80 (MQ=255)       gCTGGGAAGAAGCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTg                                                                    >  2:176732/1‑80 (MQ=255)           ggAAGAAGCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACt                                                                <  2:176664/80‑1 (MQ=255)             aagaagCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATg                                                             >  2:176707/1‑81 (MQ=255)             aagaagCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATg                                                             >  1:176695/1‑81 (MQ=255)               gaagCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATgg                                                            >  2:176713/1‑80 (MQ=255)               gaagCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATgg                                                            >  2:176708/1‑80 (MQ=255)                aagCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCa                                                          <  1:176665/81‑1 (MQ=255)                 agCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCAc                                                         >  1:176706/1‑81 (MQ=255)                  gCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCAcc                                                        <  1:176666/81‑1 (MQ=255)                    cTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCgg                                                      >  1:176698/1‑81 (MQ=255)                      aCGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTa                                                    >  1:176694/1‑81 (MQ=255)                       cGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGt                                                     <  2:176667/79‑1 (MQ=255)                        gACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACg                                                  >  2:176701/1‑81 (MQ=255)                         aCATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGc                                                 >  1:176697/1‑81 (MQ=255)                              atcGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGggc                                            <  2:176668/81‑1 (MQ=255)                                 gcgACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGggcgg                                          >  2:176710/1‑80 (MQ=255)                                 gcgACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGggcagc                                         <  2:176669/81‑1 (MQ=255)                                  cgACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCa                                        <  1:176670/81‑1 (MQ=255)                                                     gATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTCCTGgcc                     >  1:176720/1‑81 (MQ=255)                                                      aTCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTCCTGgc                      <  1:176672/79‑1 (MQ=255)                                                      aTCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTCCTGgc                      <  1:176671/79‑1 (MQ=255)                                                         aaaGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTCCTGgccgcc                  <  1:176673/80‑1 (MQ=255)                                                         aaaGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTCCTGGccgccg                 >  1:176702/1‑81 (MQ=255)                                                            gagTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGa              <  2:176674/81‑1 (MQ=255)                                                                 cGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATAc         <  1:176675/81‑1 (MQ=255)                                                                  ggCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATAcc        >  2:176723/1‑81 (MQ=255)                                                                     aaCGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTg     >  2:176705/1‑81 (MQ=255)                                                                      aCGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTgg    >  1:176721/1‑81 (MQ=255)                                                                      aCGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTgg    <  1:176676/81‑1 (MQ=255)                                                                        gAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTc  <  1:176677/81‑1 (MQ=255)                                                                           |                                                                              GCATTAGCTGGGAAGAAGCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGTCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTC  >  NC_000913/941368‑941519

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