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%	87.2 / NA	25	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 (10/15);  total (10/15)

AATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGG  >  NC_000913/118981‑119074                                                |                                               aaTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcct                                              <  1:1095129/50‑1 (MQ=255) aaTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcct                                              >  1:1236585/1‑50 (MQ=255)   tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                            <  1:1945197/50‑1 (MQ=255)    aGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGCAGcctct                                            >  1:1405743/1‑49 (MQ=38)      tCCAGTATGTTCCTTTCGAGAAACGTGCATGGCATGCGGGAGCCTCTCAg                                         <  1:1256040/50‑1 (MQ=255)        cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGt                                        <  1:766220/49‑1 (MQ=255)        cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTa                                       >  1:547308/1‑50 (MQ=255)         aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                      >  1:1281251/1‑50 (MQ=255)          gTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                      >  1:613465/1‑49 (MQ=255)             tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                   <  1:211106/49‑1 (MQ=255)             tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                   <  1:1442635/49‑1 (MQ=255)               ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgggg                                >  1:907701/1‑50 (MQ=255)               ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                                 <  1:1256927/49‑1 (MQ=255)                     tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGaa                          >  1:1454144/1‑50 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                        >  1:410710/1‑50 (MQ=255)                         tAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCt                      <  1:202258/50‑1 (MQ=255)                         tAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCt                      <  1:1064752/50‑1 (MQ=255)                          aaaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                     >  1:132506/1‑50 (MQ=255)                           aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                     >  1:458845/1‑49 (MQ=255)                            aCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCa                   <  1:1846871/50‑1 (MQ=255)                                gcatggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGa               >  1:331530/1‑50 (MQ=255)                                 catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt              <  1:2005045/50‑1 (MQ=255)                                  atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt             >  1:2070146/1‑50 (MQ=255)                                     gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttttt          <  1:407311/50‑1 (MQ=255)                                     gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtttt           <  1:988630/49‑1 (MQ=255)                                          cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAt      <  1:2231211/49‑1 (MQ=255)                                            ggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTgg   <  1:1906666/50‑1 (MQ=255)                                             gAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTggg  <  1:2051894/50‑1 (MQ=255)                                             gAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTggg  <  1:183614/50‑1 (MQ=255)                                                |                                               AATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGG  >  NC_000913/118981‑119074

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