Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,793,867	C→A	A378D (GCT→GAT)	gabT →	4‑aminobutyrate aminotransferase GabT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,793,867	0	C	A	100.0%	128.0 / NA	34	A378D (GCT→GAT)	gabT	4‑aminobutyrate aminotransferase GabT
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (15/19);  total (15/19)

CGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAT  >  NC_000913/2793820‑2793913                                                |                                               cGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATc                                              <  1:1198434/50‑1 (MQ=255)  ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcg                                             >  1:1246481/1‑50 (MQ=255)  ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcg                                             >  1:629683/1‑50 (MQ=255)  ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcg                                             <  1:1936714/50‑1 (MQ=255)    cGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcgcg                                           >  1:334833/1‑50 (MQ=255)     gATCACAACAAGCCGGGCGCCAAACTCACCGCCGAGATCGTGGATCgcgc                                          >  1:544627/1‑50 (MQ=255)      aTCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGcc                                         <  1:2005956/50‑1 (MQ=255)       tCAAAACAAGCCGGACGCCAAACTCGCCGCCGAGATCGTGGATCGCGccc                                        <  1:1981841/50‑1 (MQ=255)         acaacaAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                      >  1:616063/1‑50 (MQ=255)         acaacaAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                      <  1:929029/50‑1 (MQ=255)         acaacaAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                      <  1:849771/50‑1 (MQ=255)            acaaGCCGGACGCCAAACCAACCGCCGAGATCGTGGATCGCGCCCGCGAt                                   <  1:2156385/50‑1 (MQ=255)              aaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGCTa                                  >  1:1150639/1‑49 (MQ=255)              aaGCCGGACGAAAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATaa                                 <  1:472426/50‑1 (MQ=255)               aGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATaaa                                <  1:986895/50‑1 (MQ=255)                gCCGAACGCCAAACACACCGCCGAGATCGTGGATCGCGCCCGCGATAAAg                               <  1:2004417/50‑1 (MQ=255)                 ccGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAgg                              >  1:165745/1‑50 (MQ=255)                 ccGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAgg                              >  1:1500519/1‑50 (MQ=255)                   ggACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGcc                            >  1:365548/1‑50 (MQ=255)                    gACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt                           >  1:1365544/1‑50 (MQ=255)                     aCGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt                           >  1:135780/1‑49 (MQ=255)                         cAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttc                      <  1:1993490/50‑1 (MQ=255)                         aaaaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttc                      <  1:83753/49‑1 (MQ=255)                          aaaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttct                     <  1:989833/50‑1 (MQ=255)                           aaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttctt                    <  1:288049/50‑1 (MQ=255)                           aaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttctt                    <  1:2232915/50‑1 (MQ=255)                                aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctct                <  1:463933/49‑1 (MQ=255)                                      gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg         >  1:2008491/1‑50 (MQ=255)                                           cGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg     >  1:1414186/1‑49 (MQ=255)                                            gTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGt    <  1:719580/49‑1 (MQ=255)                                            gTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTa   >  1:263057/1‑50 (MQ=255)                                            gTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTa   >  1:234251/1‑50 (MQ=255)                                             tGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAt  <  1:654778/50‑1 (MQ=255)                                             tGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAt  <  1:781115/50‑1 (MQ=255)                                                |                                               CGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAT  >  NC_000913/2793820‑2793913

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