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%	143.4 / NA	37	A378D (GCT→GAT)	gabT	4‑aminobutyrate aminotransferase GabT
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (19/18);  total (19/18)

CGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAT  >  NC_000913/2793820‑2793913                                                |                                               cGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATc                                              <  1:4616981/50‑1 (MQ=255)  ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcg                                             <  1:2064370/50‑1 (MQ=255)     gATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCgcgc                                          >  1:4214536/1‑50 (MQ=255)      aTCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGcc                                         >  1:2319891/1‑50 (MQ=255)         acaacaAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                      >  1:5812878/1‑50 (MQ=255)          caacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgcg                                     <  1:3463083/50‑1 (MQ=255)          caacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                      >  1:3057363/1‑49 (MQ=255)             caaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATa                                  <  1:2096214/50‑1 (MQ=255)              aaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATaa                                 <  1:1560479/50‑1 (MQ=255)               aGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATa                                  >  1:824655/1‑48 (MQ=255)                  cGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc                             <  1:1130468/50‑1 (MQ=255)                  cGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc                             >  1:1009750/1‑50 (MQ=255)                   ggACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGcc                            >  1:1932972/1‑50 (MQ=255)                   ggACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGcc                            >  1:3499003/1‑50 (MQ=255)                     aCGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt                           >  1:3134655/1‑49 (MQ=255)                     aCGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt                           >  1:5155875/1‑49 (MQ=255)                     aCGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTg                          <  1:2430753/50‑1 (MQ=255)                      cGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGa                         >  1:2678048/1‑50 (MQ=255)                       gCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAt                        <  1:2342668/50‑1 (MQ=255)                        ccAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAtt                       <  1:220268/50‑1 (MQ=255)                           aaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttctt                    >  1:5303738/1‑50 (MQ=255)                            aCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttcttc                   <  1:4058366/50‑1 (MQ=255)                                aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctct                <  1:1139270/49‑1 (MQ=255)                                  cgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCt             <  1:3759837/50‑1 (MQ=255)                                   gccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTg            <  1:1280069/50‑1 (MQ=255)                                   gccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTg            >  1:1118038/1‑50 (MQ=255)                                     cgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCg          <  1:3249796/50‑1 (MQ=255)                                     cgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCg          <  1:5408489/50‑1 (MQ=255)                                      gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg         >  1:4539261/1‑50 (MQ=255)                                      gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg         >  1:5981293/1‑50 (MQ=255)                                      gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg         >  1:1297885/1‑50 (MQ=255)                                       agaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGc        >  1:5421844/1‑50 (MQ=255)                                       agaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGc        >  1:827683/1‑50 (MQ=255)                                        gaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGcc       >  1:4749447/1‑50 (MQ=255)                                         aTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc      <  1:309253/50‑1 (MQ=255)                                             tGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAt  <  1:671021/50‑1 (MQ=255)                                             tGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAt  <  1:702487/50‑1 (MQ=255)                                                |                                               CGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAT  >  NC_000913/2793820‑2793913

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