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

CGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAT  >  NC_000913/2793820‑2793913                                                |                                               cGGGGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATc                                              <  2:690558/50‑1 (MQ=255) cGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATc                                              >  2:3075744/1‑50 (MQ=255)  ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcg                                             >  2:2806460/1‑50 (MQ=255)  ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcg                                             >  2:1718068/1‑50 (MQ=255)  ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcg                                             >  2:1582686/1‑50 (MQ=255)  ggCGATCACAACAAGCCGAACGCCAACCTCACCGACGAGATCGTGGATca                                             <  2:3719082/50‑2 (MQ=255)   gCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcgc                                            >  2:1934749/1‑50 (MQ=255)     gATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCgcgc                                          >  2:5645793/1‑50 (MQ=255)     gATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCgcgc                                          >  1:3886169/1‑50 (MQ=255)     gATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCgcgc                                          <  1:3517701/50‑1 (MQ=255)     gATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCgcgc                                          >  1:3339875/1‑50 (MQ=255)      aTCACAACAAGCCGGACGCCGAACTCACCGCCGAGATCGTGGATCGCGcc                                         >  1:366065/1‑50 (MQ=38)      aTCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGcc                                         >  2:327206/1‑50 (MQ=255)      aTCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGcc                                         >  1:3704052/1‑50 (MQ=255)      aTCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGcc                                         >  2:22934/1‑50 (MQ=255)       tCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGccc                                        <  1:50502/50‑1 (MQ=255)        cacaACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcg                                       >  2:2749241/1‑50 (MQ=255)          caacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                      >  1:5491697/1‑49 (MQ=255)              aaGCCGGACGCCAAGCTCACCGCCGAGATCGTGGATCGCGCCCGCGATaa                                 <  2:2928221/50‑1 (MQ=255)              aaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATa                                  >  1:394340/1‑49 (MQ=255)              aaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATa                                  >  2:5427235/1‑49 (MQ=255)               aGCCGAACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATaaa                                <  1:4080494/50‑1 (MQ=255)                gCCGGACGCAAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAg                               <  1:3804407/50‑1 (MQ=255)                 ccGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAgg                              >  1:5016786/1‑50 (MQ=255)                  cGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc                             >  1:3454458/1‑50 (MQ=255)                  cGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc                             >  1:3953555/1‑50 (MQ=255)                  cGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc                             >  2:3121713/1‑50 (MQ=255)                  cGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc                             <  1:3896413/50‑1 (MQ=255)                   ggACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGcc                            >  1:4732122/1‑50 (MQ=255)                   ggACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGcc                            >  1:463702/1‑50 (MQ=255)                    gACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt                           >  1:118316/1‑50 (MQ=255)                    gACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt                           >  1:4990301/1‑50 (MQ=255)                    gACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt                           >  1:1965194/1‑50 (MQ=255)                    gACGCAAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt                           <  2:3185023/50‑1 (MQ=255)                     aCGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt                           >  2:2818916/1‑49 (MQ=255)                      cGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGa                         <  2:526611/50‑1 (MQ=255)                      cGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGa                         >  1:2070210/1‑50 (MQ=255)                       gCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAt                        >  2:5399246/1‑50 (MQ=255)                         cAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttc                      >  2:418325/1‑50 (MQ=255)                         cAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttc                      >  1:596423/1‑50 (MQ=255)                         cAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttc                      >  2:4105957/1‑50 (MQ=255)                         cAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttc                      >  1:4673910/1‑50 (MQ=255)                          aaaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttct                     >  1:2514536/1‑50 (MQ=255)                          aaaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttct                     >  1:5125634/1‑50 (MQ=255)                             cTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATtcttct                  <  1:4953608/50‑1 (MQ=255)                             cTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATtcttct                  >  1:1350322/1‑50 (MQ=255)                             cTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATtcttct                  <  1:2433512/50‑1 (MQ=255)                               cACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctct                >  2:1729291/1‑50 (MQ=255)                               cACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctct                <  1:2204085/50‑1 (MQ=255)                                aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctctc               >  1:1649138/1‑50 (MQ=255)                                aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctct                <  1:4818840/49‑1 (MQ=255)                                aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctc                 <  2:1198390/48‑1 (MQ=255)                                 ccgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTcc              <  1:2549334/50‑1 (MQ=255)                                  cgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCt             <  2:1103308/50‑1 (MQ=255)                                  cgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCt             <  1:2039438/50‑1 (MQ=255)                                  cgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCt             <  2:2869390/50‑1 (MQ=255)                                  cgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCt             <  2:5106534/50‑1 (MQ=255)                                   gccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTg            >  1:3532401/1‑50 (MQ=255)                                   gccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTg            <  2:4783670/50‑1 (MQ=255)                                   gccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTg            >  2:5077428/1‑50 (MQ=255)                                   gccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTg            >  2:1422588/1‑50 (MQ=255)                                   gccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTg            >  1:925490/1‑50 (MQ=255)                                   gccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTg            <  1:2962052/50‑1 (MQ=255)                                    ccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGc           >  1:4064104/1‑50 (MQ=255)                                     cgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCg          >  1:1475568/1‑50 (MQ=255)                                      gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg         <  1:2182214/50‑1 (MQ=255)                                      gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg         >  1:2057599/1‑50 (MQ=255)                                      gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg         >  2:4112559/1‑50 (MQ=255)                                      gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg         >  2:5593530/1‑50 (MQ=255)                                      gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg         <  1:3619261/50‑1 (MQ=255)                                      gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg         >  1:3260652/1‑50 (MQ=255)                                       agaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGc        >  1:4017308/1‑50 (MQ=255)                                        gaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGcc       >  1:1623435/1‑50 (MQ=255)                                         aTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc      <  1:3170201/50‑1 (MQ=255)                                         aTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc      >  2:3182507/1‑50 (MQ=255)                                         aTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc      <  1:3347643/50‑1 (MQ=255)                                         aTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc      <  1:2406268/50‑1 (MQ=255)                                          tCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg     >  2:4101472/1‑50 (MQ=255)                                          tCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg     >  1:4403055/1‑50 (MQ=255)                                          tCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg     >  2:4852038/1‑50 (MQ=255)                                          tCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg     >  1:2234725/1‑50 (MQ=255)                                           cGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg     >  2:2315965/1‑49 (MQ=255)                                            gTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGt    <  2:40527/49‑1 (MQ=255)                                            gTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGt    <  1:1042074/49‑1 (MQ=255)                                            gTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTa   <  1:2457690/50‑1 (MQ=255)                                            gTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTa   <  1:4575184/50‑1 (MQ=255)                                             tGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAt  >  2:2184159/1‑50 (MQ=255)                                             tGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAt  >  2:960757/1‑50 (MQ=255)                                                |                                               CGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAT  >  NC_000913/2793820‑2793913

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