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

GGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTA  >  NC_000913/2793821‑2793912                                               |                                              ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATc                                             >  1:1724149/1‑49 (MQ=255)  gCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcg                                            >  1:2061355/1‑49 (MQ=255)   cGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcgc                                           >  1:2626162/1‑49 (MQ=255)    gATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcgcg                                          >  1:1113624/1‑49 (MQ=255)      tCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGcc                                        >  1:2214082/1‑49 (MQ=255)         caacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                     >  1:136967/1‑49 (MQ=255)           acaaGCCGCACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGa                                   >  1:1161946/1‑49 (MQ=255)               gCCGGACGCCAACCTCACCGCCGAGATCGTGGATCGCGCCCGCGATaaa                               <  1:926471/49‑1 (MQ=39)               gCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATa                                 >  1:280976/1‑47 (MQ=255)                ccGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAg                              >  1:998901/1‑49 (MQ=255)                ccGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAg                              <  1:106312/49‑1 (MQ=255)                 cGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAgg                             <  1:469360/49‑1 (MQ=255)                 cGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAgg                             <  1:1321674/49‑1 (MQ=255)                 cGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAgg                             >  1:339638/1‑49 (MQ=255)                  ggACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc                            <  1:1276858/49‑1 (MQ=255)                  ggACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc                            >  1:2231917/1‑49 (MQ=255)                  ggACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc                            >  1:2602286/1‑49 (MQ=255)                     cGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTg                         <  1:1863637/49‑1 (MQ=255)                      gCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGa                        >  1:2727435/1‑49 (MQ=255)                       ccAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAt                       >  1:1423862/1‑49 (MQ=255)                        cAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAtt                      >  1:744524/1‑49 (MQ=255)                        aaaaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAtt                      <  1:2548799/48‑1 (MQ=255)                             tCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATtcttct                 >  1:112512/1‑49 (MQ=255)                               aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATtcttct                 <  1:1496550/47‑1 (MQ=255)                               aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctct               >  1:520169/1‑49 (MQ=255)                                   ccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTg           >  1:1235996/1‑49 (MQ=255)                                    cgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGc          <  1:2130408/49‑1 (MQ=255)                                    cgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGc          <  1:106527/49‑1 (MQ=255)                                     gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCg         >  1:2356714/1‑49 (MQ=255)                                     gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCg         >  1:1531573/1‑49 (MQ=255)                                     gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCg         >  1:1087883/1‑49 (MQ=255)                                        aTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGcc      >  1:1121872/1‑49 (MQ=255)                                         tCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc     <  1:67702/49‑1 (MQ=255)                                            tGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTa  >  1:1548582/1‑49 (MQ=255)                                               |                                              GGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTA  >  NC_000913/2793821‑2793912

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