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

GGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTA  >  NC_000913/2793821‑2793912                                               |                                              ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATc                                             <  1:1070337/49‑1 (MQ=255) ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATc                                             <  1:2078755/49‑1 (MQ=255)     aTCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCgcgc                                         <  1:554396/49‑1 (MQ=255)     aTCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCgcgc                                         >  1:1245845/1‑49 (MQ=255)     aTCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCgcgc                                         >  1:1264223/1‑49 (MQ=255)      tCACAACAAGCCGGACGCCAACCTCACCGCCGAGATCGTGGATCGCGcc                                        <  1:475348/49‑1 (MQ=255)       cacaACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGccc                                       >  1:490599/1‑49 (MQ=255)       cacaACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGccc                                       >  1:1596043/1‑49 (MQ=255)        acaacaAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcg                                      >  1:1150634/1‑49 (MQ=255)         caacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                     >  1:904299/1‑49 (MQ=255)          aacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                     >  1:659170/1‑48 (MQ=255)          aacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc                                     >  1:1847740/1‑48 (MQ=255)            caaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGa                                   >  1:157328/1‑48 (MQ=255)            caaGCCGGACGCAAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGAt                                  <  1:711789/49‑1 (MQ=38)             aaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATa                                 <  1:2100320/49‑1 (MQ=255)              aGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATa                                 >  1:1554921/1‑48 (MQ=255)                ccGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAg                              >  1:2338839/1‑49 (MQ=255)                 cGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAgg                             >  1:529402/1‑49 (MQ=255)                  ggACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc                            >  1:1917832/1‑49 (MQ=255)                   gACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGcc                           >  1:1818882/1‑49 (MQ=255)                   gACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGcc                           >  1:1262934/1‑49 (MQ=255)                     cGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt                          >  1:1544280/1‑48 (MQ=255)                     cGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTg                         <  1:1978812/49‑1 (MQ=255)                     cGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTg                         <  1:1332772/49‑1 (MQ=255)                      gCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGa                        >  1:228543/1‑49 (MQ=255)                        cAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAtt                      >  1:1225527/1‑49 (MQ=255)                         aaaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttc                     >  1:1046202/1‑49 (MQ=255)                         aaaCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttc                     >  1:2115374/1‑49 (MQ=255)                               aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctct               <  1:330009/49‑1 (MQ=255)                               aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctct               <  1:1200729/49‑1 (MQ=255)                                ccgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctctc              <  1:2209305/49‑1 (MQ=255)                                ccgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctctc              >  1:851850/1‑49 (MQ=255)                                  gccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCt            <  1:1087040/49‑1 (MQ=255)                                    cgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGc          >  1:755272/1‑49 (MQ=255)                                      agaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg        <  1:1264623/49‑1 (MQ=255)                                         tCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc     >  1:1928307/1‑49 (MQ=255)                                         tCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc     <  1:1668204/49‑1 (MQ=255)                                          cGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg    >  1:262653/1‑49 (MQ=255)                                           gTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg    <  1:649330/48‑1 (MQ=255)                                            tGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTa  <  1:1664805/49‑1 (MQ=255)                                               |                                              GGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTA  >  NC_000913/2793821‑2793912

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