Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,791,708	C→T	P146S (CCG→TCG)	gabD →	succinate‑semialdehyde dehydrogenase (NADP(+)) GabD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,791,708	0	C	T	100.0%	27.9 / NA	11	P146S (CCG→TCG)	gabD	succinate‑semialdehyde dehydrogenase (NADP(+)) GabD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (5/6);  total (5/6)

GCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTC  >  NC_000913/2791502‑2791920                                                                                                                                                                                                               |                                                                                                                                                                                                                     gcAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGTCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCcggcn                                                                                                                                                                          <  1:64072/251‑2 (MQ=255)                          gCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGTCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCgatgat                                                                                                                                                                  >  1:38980/1‑234 (MQ=255)                                                   atgaCCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGTCGATTGGCGTCan                                                                                                                                                                                                         <  1:67873/170‑2 (MQ=255)                                                                                     cGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGTCGATTGGCGGCACCGCGGCTCTCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCgn                                                                                                                                                           <  1:48201/182‑2 (MQ=255)                                                                                                aTCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGTCGATTGGCGTCACCGCGGCTATCACGCCGGGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAg                                                                                                              <  1:131634/216‑1 (MQ=255)                                                                                                  cAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGTCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTAACCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGcgctggcgctgg                                                                         >  1:61817/1‑251 (MQ=255)                                                                                                             cctcctTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCCTTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGTCGCTTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCg                                                              >  1:94881/1‑251 (MQ=255)                                                                                                                         aGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGTCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCGGGCGATCCGCGCgggn                                                  <  1:68332/251‑2 (MQ=255)                                                                                                                                       gaagGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGTCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTg                                                                                                                 >  1:119462/1‑174 (MQ=255)                                                                                                                                                                                            ccTGATTGTTATCAAGCAGTCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCACTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGgggtg  >  1:7674/1‑231 (MQ=255)                                                                                                                                                                                                aTTGTTATCAAGCAGTCAATTGGCGTCCCCGCCGCTATCACGCCGTGGAACTGCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTgg                                                                <  1:248789/166‑1 (MQ=255)                                                                                                                                                                                                               |                                                                                                                                                                                                                     GCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTC  >  NC_000913/2791502‑2791920

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