Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,917,231	C→T	G49S (GGC→AGC)	atpD ←	ATP synthase F1 complex subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,917,231	0	C	T	100.0%	34.9 / NA	15	G49S (GGC→AGC)	atpD	ATP synthase F1 complex subunit beta
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (5/10);  total (5/10)

CTTTCATGTCGACCGGTTCACCCAGTACGTTCATGATACGGCCCAGAGTCGCTTTACCTACCGGGACTTCAATCGGGTGTTCGAGGTCTTTTACATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGC  >  NC_000913/3917078‑3917413                                                                                                                                                          |                                                                                                                                                                                       cTTTCATGTCGACCGGTTCACCCAGTACGTTCATGATTCGGCCCAGAGTCGCTTTACCCACAGGGACTTCAATCGGATGTTCGAGGTCTTTTGCATCCTGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGCTACCGCTGCTGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCACGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCg                                                                                       >  1:464672/1‑251 (MQ=255)                    aCCCAGTACGTTCATGATACGGCCCAGAGTCGCTTTACCTACCGGGACTTCAATCGGGTGTTCGAGGTCTTTTACATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCTGCCGAGCTGCTGCTGAACTTCCAGCACCAGAc                                                                                                                                                        >  1:33623/1‑167 (MQ=255)                         gTACGTTCATGATACGGCCCAGAGTCGCTTTACCTACCGGGACTTCAATCGGGTGTTCGAGGTCTTTTACATCCAGACCGCGACGCAGACCGTCGGAGGAACTCAGTGCGATGGTACGTACGATACCGCTGCCGAGCTGCGGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGc                                                                                                       <  1:229249/211‑1 (MQ=255)                                       cGGCCCAGAGTCGCTTTACCTACCGGGACTTCAATCGGGTGTTCGAGGTCTTTTACATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCTGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTg                                                                                                 <  1:16983/203‑1 (MQ=255)                                                      ttACCTACCGGGACTTCAATCGGGTGTTCGAGGTCTTTTACATCCAGACCGCGACGCAGACCGTCGTAGGACCCCACTGCGATGGTACGTACGATACCGCTGCCGAGCTGCTGCTGAACTTCCAgn                                                                                                                                                               <  1:89507/126‑2 (MQ=255)                                                                                       tCTTTTACATCCAGAACGCGGCGCAGACCGTCGGGGGAACCCATTGCGCTGGACCGTACGATACGGCTGACGAGGTGCTGCTGCACTTCCAGGACCAGACGCTCATTACCATTTTGCCCCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTgg                                                            >  1:223853/1‑192 (MQ=255)                                                                                          nttACATNNAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCTGCCGAGCTGCTGCTGAACTTCCAGCACCAGAc                                                                                                                                                        <  2:33623/96‑1 (MQ=255)                                                                                                        gcgACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCTGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCGCAAGAGCATCGTACACTCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAAc                  <  1:267731/217‑1 (MQ=255)                                                                                                         cgACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCTGCCGAGCTGCTGCTGAACTGCCAGCACCAGACGCTCATTACCATTTTGCCCCTCAAGAGCATCGTACAcg                                                                                                                  <  1:186227/120‑1 (MQ=255)                                                                                                                             aCCCATGGCGATTGTACGTACGATACCGCTGCCAAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCTTTTCGCACCTCAAGATCATCGTACACGCGCGGTACg                                                                                                         <  1:173208/109‑1 (MQ=255)                                                                                                                              cccATTGCGATGGTACTTACGATACCGCTGCCGAGCTGCTGCGGAACTTCCAGCACCAGACGCTCATTACCAGTTGGCACCTCAAGAGCATCGTACACGCGCgg                                                                                                             <  1:320699/104‑1 (MQ=255)                                                                                                                                           tacgtacgATACCGCTGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACGATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAAtt                                                                                         >  1:294632/1‑111 (MQ=255)                                                                                                                                                 cgATACCGCTGCCGAGCTGCTGCTGANCTTCCANNACNAGCCGCTCATTACCATTTTNc                                                                                                                                       >  2:114117/1‑59 (MQ=255)                                                                                                                                                 cgATACCGCTGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGc                                                                                                                                       <  1:114117/59‑1 (MQ=255)                                                                                                                                                   aTACCGCTGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCgcggn  <  1:406456/190‑2 (MQ=255)                                                                                                                                                          |                                                                                                                                                                                       CTTTCATGTCGACCGGTTCACCCAGTACGTTCATGATACGGCCCAGAGTCGCTTTACCTACCGGGACTTCAATCGGGTGTTCGAGGTCTTTTACATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGC  >  NC_000913/3917078‑3917413

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