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

CGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGT  >  NC_000913/2793820‑2793911                                                |                                             cGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATc                                            <  1:1538361/50‑1 (MQ=255) cGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATc                                            >  1:3454695/1‑50 (MQ=255) cGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATc                                            >  1:2667061/1‑50 (MQ=255)  ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcg                                           <  1:3565373/50‑1 (MQ=255)   gCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcgc                                          <  1:25509/50‑1 (MQ=255)   acGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcgc                                          >  1:1436307/2‑50 (MQ=255)     gATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCgcgc                                        >  1:969788/1‑50 (MQ=255)     gATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCgcgc                                        >  1:1627924/1‑50 (MQ=255)     gATCACAACAAACCGGACGCCAAACTCACCGCCGAGATCGTGGATCgcgc                                        >  1:304187/1‑50 (MQ=255)      aTCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGcc                                       >  1:3453501/1‑50 (MQ=255)       tCACAACAAGCGGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGccc                                      <  1:463735/50‑1 (MQ=255)       tCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGccc                                      >  1:512672/1‑50 (MQ=255)       tCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGccc                                      >  1:2014271/1‑50 (MQ=255)         acaacaAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                    >  1:2979673/1‑50 (MQ=255)         acaacaAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                    <  1:402920/50‑1 (MQ=255)          caacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgcg                                   <  1:427646/50‑1 (MQ=255)          caacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgcg                                   <  1:1252193/50‑1 (MQ=255)          caacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                    >  1:1745786/1‑49 (MQ=255)           aacaaGCCGGACGCAAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGa                                  <  1:1209235/50‑1 (MQ=39)              aaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATa                                >  1:2286625/1‑49 (MQ=255)               aGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATaaa                              <  1:477614/50‑1 (MQ=255)                gCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAg                             <  1:499053/50‑1 (MQ=255)                gCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAg                             >  1:2858793/1‑50 (MQ=255)                gCCGGACCCAAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAg                             <  1:2785559/50‑1 (MQ=255)                 ccGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAgg                            >  1:2884909/1‑50 (MQ=255)                 ccGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAgg                            <  1:2101548/50‑1 (MQ=255)                  cGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc                           >  1:2007917/1‑50 (MQ=255)                  cGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc                           >  1:3599307/1‑50 (MQ=255)                   ggACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGcc                          >  1:629779/1‑50 (MQ=255)                   ggACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGcc                          >  1:838199/1‑50 (MQ=255)                    gACGCAAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt                         <  1:699874/50‑1 (MQ=255)                     aCGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt                         >  1:171965/1‑49 (MQ=255)                      cGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGa                       <  1:3623170/50‑1 (MQ=255)                      cGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGa                       >  1:753357/1‑50 (MQ=255)                        ccAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAtt                     >  1:3424775/1‑50 (MQ=255)                          aaaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttct                   >  1:1701985/1‑50 (MQ=255)                          aaaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttct                   >  1:3203147/1‑50 (MQ=255)                          aaaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttct                   >  1:2477297/1‑50 (MQ=255)                           aaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttctt                  <  1:3028230/50‑1 (MQ=255)                           aaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttctt                  <  1:963492/50‑1 (MQ=255)                           aaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttctt                  <  1:2507318/50‑1 (MQ=255)                            aCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttcttc                 <  1:2436409/50‑1 (MQ=255)                               cACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctct              >  1:255840/1‑50 (MQ=255)                                aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctct              <  1:749499/49‑1 (MQ=255)                                  cgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCt           <  1:1814962/50‑1 (MQ=255)                                   gccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTg          >  1:1834174/1‑50 (MQ=255)                                    ccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGc         <  1:3921128/50‑1 (MQ=255)                                     cgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCg        >  1:3155011/1‑50 (MQ=255)                                      gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg       >  1:3192773/1‑50 (MQ=255)                                      gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg       >  1:1906059/1‑50 (MQ=255)                                         aTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc    <  1:3442974/50‑1 (MQ=255)                                         aTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc    <  1:1050545/50‑1 (MQ=255)                                          tCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg   <  1:3170858/50‑1 (MQ=255)                                           cGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg   >  1:2168581/1‑49 (MQ=255)                                            gTGGATCGCGCCCGGGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGt  <  1:1262040/49‑1 (MQ=38)                                                |                                             CGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGT  >  NC_000913/2793820‑2793911

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