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

CGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAT  >  NC_000913/2793820‑2793913                                                |                                               cGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATc                                              <  1:3332689/50‑1 (MQ=255)  ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcg                                             >  1:2885844/1‑50 (MQ=255)  ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcg                                             >  1:3112841/1‑50 (MQ=255)  ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcg                                             >  1:1892844/1‑50 (MQ=255)  ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcg                                             >  1:3190889/1‑50 (MQ=255)   gCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcgc                                            <  1:1210246/50‑1 (MQ=255)   gCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcgc                                            <  1:4359843/50‑1 (MQ=255)   gCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcgc                                            <  1:2884606/50‑1 (MQ=255)   gCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcgc                                            >  1:2551342/1‑50 (MQ=255)     gATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCgcgc                                          >  1:554297/1‑50 (MQ=255)     gATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCgcgc                                          <  1:3603243/50‑1 (MQ=255)     gATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCgcgc                                          >  1:3296936/1‑50 (MQ=255)      aTCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGcc                                         >  1:1999396/1‑50 (MQ=255)         acaacaAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                      >  1:1888860/1‑50 (MQ=255)         acaacaAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                      <  1:1606012/50‑1 (MQ=255)         aaaacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                      <  1:3919166/48‑1 (MQ=255)          caacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                      >  1:2976365/1‑49 (MQ=255)            acaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGa                                    >  1:4265958/1‑49 (MQ=255)              aaGCCGCACGCAAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATaa                                 <  1:2527082/50‑1 (MQ=255)               aGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATaaa                                <  1:585213/50‑1 (MQ=255)               aGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATa                                  >  1:4043925/1‑48 (MQ=255)                 ccGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAgg                              >  1:3665781/1‑50 (MQ=255)                 ccGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAgg                              <  1:1472650/50‑1 (MQ=255)                  cGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc                             >  1:1959229/1‑50 (MQ=255)                     aCGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTg                          <  1:570989/50‑1 (MQ=255)                     aCGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTg                          <  1:3356122/50‑1 (MQ=255)                       gCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAt                        <  1:34847/50‑1 (MQ=255)                          aaaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttct                     >  1:3046578/1‑50 (MQ=255)                          aaaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttct                     >  1:2618505/1‑50 (MQ=255)                          aaaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttct                     >  1:2396226/1‑50 (MQ=255)                          aaaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttct                     >  1:1506719/1‑50 (MQ=255)                              tCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctc                 >  1:750203/1‑50 (MQ=255)                               cACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctct                >  1:4370883/1‑50 (MQ=255)                                aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctctc               <  1:1814750/50‑1 (MQ=255)                                aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctct                <  1:144789/49‑1 (MQ=255)                                aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctc                 <  1:4056512/48‑1 (MQ=255)                                   gccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTg            <  1:4241030/50‑1 (MQ=255)                                   gccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTg            <  1:542082/50‑1 (MQ=255)                                      gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg         <  1:4289671/50‑1 (MQ=255)                                        gaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGcc       <  1:3619964/50‑1 (MQ=255)                                         aTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc      <  1:2325635/50‑1 (MQ=255)                                           cGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg     >  1:3075717/1‑49 (MQ=255)                                           cGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg     >  1:1759556/1‑49 (MQ=255)                                            gTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTa   <  1:2275605/50‑1 (MQ=255)                                             tGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAt  >  1:1433805/1‑50 (MQ=255)                                                |                                               CGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAT  >  NC_000913/2793820‑2793913

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