Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	119,028	T→C	V99A (GTC→GCC)	ampD →	1,6‑anhydro‑N‑acetylmuramoyl‑L‑alanine amidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	119,028	0	T	C	100.0%	162.1 / NA	43	V99A (GTC→GCC)	ampD	1,6‑anhydro‑N‑acetylmuramoyl‑L‑alanine amidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/20);  total (23/20)

AATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGG  >  NC_000913/118981‑119074                                                |                                               aaTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcct                                              >  1:2109822/1‑50 (MQ=255)     gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctctc                                           <  1:439770/49‑1 (MQ=255)     gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                            >  1:3772834/1‑48 (MQ=255)      tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                         >  1:1153453/1‑50 (MQ=255)      tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                         <  1:4084756/50‑1 (MQ=255)      tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                         >  1:2372348/1‑50 (MQ=255)        cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGt                                        <  1:1328299/49‑1 (MQ=255)        cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTa                                       >  1:1684186/1‑50 (MQ=255)        cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTa                                       <  1:3504933/50‑1 (MQ=255)         aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                      >  1:1361888/1‑50 (MQ=255)            aTGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                   >  1:1079529/1‑50 (MQ=255)             tGTTCCTTTCTATAACCGTGCATGGCATGCGGGAGCCTCTCAGTATCatc                                  >  1:802305/1‑48 (MQ=255)             tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                   <  1:4051151/49‑1 (MQ=255)                 ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                               <  1:3592472/49‑1 (MQ=255)                 ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                               <  1:1209311/49‑1 (MQ=255)                  cTTTCGATAAACGTGCATGGCATGCGGGGGCCTCTCAGTATCAGGGgcg                              <  1:281714/49‑1 (MQ=255)                      cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                         >  1:200554/1‑50 (MQ=255)                      cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                         >  1:4215539/1‑50 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                         <  1:2223935/49‑1 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                         <  1:3489079/49‑1 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                        >  1:811211/1‑50 (MQ=255)                         tAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCt                      >  1:2557677/1‑50 (MQ=255)                         tAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCt                      <  1:949321/50‑1 (MQ=255)                         tAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCt                      >  1:3607110/1‑50 (MQ=255)                          aaaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                     >  1:1685564/1‑50 (MQ=255)                           aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                     <  1:2402437/49‑1 (MQ=255)                            aCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGc                    >  1:3863855/1‑49 (MQ=255)                              gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                 <  1:2488159/50‑1 (MQ=255)                               tGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                 >  1:812986/1‑49 (MQ=255)                                 catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGa               <  1:2369629/49‑1 (MQ=255)                                 catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGa               <  1:3986476/49‑1 (MQ=255)                                  atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt             >  1:98380/1‑50 (MQ=255)                                  atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt             >  1:2603831/1‑50 (MQ=255)                                     gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttttt          <  1:1770104/50‑1 (MQ=255)                                     gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttttt          <  1:846266/50‑1 (MQ=255)                                     gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttttt          >  1:1706547/1‑50 (MQ=255)                                     gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtttt           <  1:2457864/49‑1 (MQ=255)                                      catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc         >  1:2730281/1‑50 (MQ=255)                                       atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc         >  1:4309960/1‑49 (MQ=255)                                       atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc         >  1:772683/1‑49 (MQ=255)                                         gCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAt      >  1:2999032/1‑50 (MQ=255)                                           gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg    >  1:2525120/1‑50 (MQ=255)                                           gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg    >  1:917613/1‑50 (MQ=255)                                            ggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTgg   <  1:3924075/50‑1 (MQ=255)                                            ggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTgg   <  1:3126071/50‑1 (MQ=255)                                             gAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTggg  <  1:1629793/50‑1 (MQ=255)                                                |                                               AATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGG  >  NC_000913/118981‑119074

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: ‑