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%	169.9 / 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)

AATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATT  >  NC_000913/118981‑119071                                                |                                            aaTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcct                                           <  1:1947891/50‑1 (MQ=255)  aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                          >  1:4561039/1‑50 (MQ=255)   tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                         >  1:3120338/1‑50 (MQ=255)    aGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                         >  1:3451222/1‑49 (MQ=255)    aGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                         >  1:5659752/1‑49 (MQ=255)    aGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                         >  1:3035261/1‑49 (MQ=255)    aGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                         >  1:4974482/1‑49 (MQ=255)     gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctctc                                        <  1:1231110/49‑1 (MQ=255)     gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCa                                       <  1:4914899/50‑1 (MQ=255)      tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                      <  1:3990538/50‑1 (MQ=255)        cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTa                                    >  1:5058268/1‑50 (MQ=255)           tATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCa                                 <  1:2511474/50‑1 (MQ=255)            aTGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                >  1:2794167/1‑50 (MQ=255)             tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgg                               <  1:2178350/50‑1 (MQ=255)              gTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                              >  1:2021289/1‑50 (MQ=255)              gTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                              >  1:2735531/1‑50 (MQ=255)               ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                              <  1:2373803/49‑1 (MQ=255)                 ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                            <  1:382167/49‑1 (MQ=255)                 ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                            <  1:5301103/49‑1 (MQ=255)                   tttCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                          >  1:3719312/1‑49 (MQ=255)                   tttCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                          >  1:964716/1‑49 (MQ=255)                     tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGaa                       >  1:2559655/1‑50 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                      <  1:2550827/49‑1 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                     <  1:1705076/50‑1 (MQ=255)                        aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                    >  1:1975473/1‑50 (MQ=255)                        aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                    >  1:6036985/1‑50 (MQ=255)                        aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                    >  1:5361407/1‑50 (MQ=255)                         tAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCt                   <  1:1169227/50‑1 (MQ=255)                         tAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCt                   <  1:3225585/50‑1 (MQ=255)                           aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                  >  1:5652108/1‑49 (MQ=255)                            aCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGc                 >  1:28015/1‑49 (MQ=255)                            aCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCa                <  1:5678197/50‑1 (MQ=255)                             cGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCaa               >  1:1650354/1‑50 (MQ=255)                             cGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCaa               >  1:4475662/1‑50 (MQ=255)                             cGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCaa               >  1:2800019/1‑50 (MQ=255)                             cGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCaa               >  1:3703487/1‑50 (MQ=255)                              gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt              >  1:4665484/1‑50 (MQ=255)                              gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt              >  1:954416/1‑50 (MQ=255)                               tGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCaagg             <  1:4462440/50‑3 (MQ=39)                                gcatggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGa            >  1:2064116/1‑50 (MQ=255)                                 catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGa            <  1:4466712/49‑1 (MQ=255)                                 catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt           >  1:4701134/1‑50 (MQ=255)                                 catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt           <  1:5669411/50‑1 (MQ=255)                                 catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt           <  1:4069043/50‑1 (MQ=255)                                  atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt          <  1:615436/50‑1 (MQ=255)                                     gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttttt       >  1:59765/1‑50 (MQ=255)                                      catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc      <  1:5095100/50‑1 (MQ=255)                                      catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc      >  1:1386983/1‑50 (MQ=255)                                       atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc      >  1:315365/1‑49 (MQ=255)                                          cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAtt  <  1:4202282/50‑1 (MQ=255)                                          cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAtt  <  1:3428535/50‑1 (MQ=255)                                                |                                            AATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATT  >  NC_000913/118981‑119071

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