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

GATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAT  >  NC_000913/2793824‑2793913                                            |                                               gATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCgcgc                                          >  1:1920224/1‑50 (MQ=255) gATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCgcgc                                          >  1:1757189/1‑50 (MQ=255)  aTCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGcc                                         >  1:1383662/1‑50 (MQ=255)    cacaACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcg                                       >  1:2008998/1‑50 (MQ=255)     acaacaAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcg                                       <  1:699402/49‑1 (MQ=255)       aacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGa                                    >  1:2614082/1‑50 (MQ=255)       aacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGa                                    >  1:1610356/1‑50 (MQ=255)         caaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATa                                  <  1:832596/50‑1 (MQ=255)         caaGCCGGACGCAAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATa                                  <  1:1234848/50‑1 (MQ=39)          aaGCCGGACGCAAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATaa                                 <  1:1788455/50‑1 (MQ=39)           aGCCGGACGCCAACCCCACCGCCGAGATCGTGGATCGCGCCCGCGATaaa                                <  1:225141/50‑1 (MQ=255)           aGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATaaa                                <  1:2533302/50‑1 (MQ=255)            gCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAg                               >  1:3853/1‑50 (MQ=255)               ggACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGcc                            >  1:1160181/1‑50 (MQ=255)                gACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt                           <  1:1580175/50‑1 (MQ=255)                  cGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGa                         >  1:1643686/1‑50 (MQ=255)                  cGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGa                         <  1:252358/50‑1 (MQ=255)                   gCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAt                        >  1:607008/1‑50 (MQ=255)                    ccAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAtt                       <  1:2278178/50‑1 (MQ=255)                     cAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttc                      >  1:2800866/1‑50 (MQ=255)                     cAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttc                      >  1:2441079/1‑50 (MQ=255)                      aaaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttct                     <  1:280612/50‑1 (MQ=255)                        aCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttcttc                   <  1:1304827/50‑1 (MQ=255)                         cTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATtcttct                  <  1:2240158/50‑1 (MQ=255)                           cACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctct                >  1:2148840/1‑50 (MQ=255)                            aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctctc               <  1:861177/50‑1 (MQ=255)                            aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctctc               >  1:216950/1‑50 (MQ=255)                               gccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTg            <  1:2177305/50‑1 (MQ=255)                                 cgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCg          >  1:2020063/1‑50 (MQ=255)                                  gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg         >  1:2039610/1‑50 (MQ=255)                                  gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg         <  1:2021274/50‑1 (MQ=255)                                  gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg         >  1:1607868/1‑50 (MQ=255)                                   agaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGc        >  1:1542068/1‑50 (MQ=255)                                   agaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGc        >  1:1143127/1‑50 (MQ=255)                                    gaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGcc       <  1:1610295/50‑1 (MQ=255)                                     aTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc      <  1:1303242/50‑1 (MQ=255)                                       cGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg     >  1:2808460/1‑49 (MQ=255)                                        gTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGt    <  1:762067/49‑1 (MQ=255)                                        gTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTa   >  1:1529914/1‑50 (MQ=255)                                         tGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAt  >  1:2105455/1‑50 (MQ=255)                                         tGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAt  >  1:2006861/1‑50 (MQ=255)                                            |                                               GATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAT  >  NC_000913/2793824‑2793913

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