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%	222.2 / NA	56	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 (25/31);  total (25/31)

AATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTG  >  NC_000913/118981‑119072                                                |                                             aaTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcct                                            >  1:2475432/1‑50 (MQ=255)  aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                           <  1:2630878/50‑1 (MQ=255)   tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                          <  1:993439/50‑1 (MQ=255)     gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctctc                                         <  1:2246471/49‑1 (MQ=255)      tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                       <  1:5062695/50‑1 (MQ=255)       ccAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                       >  1:2935919/1‑49 (MQ=255)        cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGt                                      <  1:125082/49‑1 (MQ=255)        cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGt                                      <  1:1347279/49‑1 (MQ=255)        cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGt                                      <  1:2532785/49‑1 (MQ=255)        cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTa                                     >  1:3264904/1‑50 (MQ=255)        cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTa                                     <  1:2706274/50‑1 (MQ=255)         aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                    <  1:1343292/50‑1 (MQ=255)          gTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                    >  1:675569/1‑49 (MQ=255)           tATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                    >  1:170568/1‑48 (MQ=255)           tATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCa                                  <  1:4284344/50‑1 (MQ=255)           tATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCa                                  <  1:2121635/50‑1 (MQ=255)           tATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCa                                  <  1:3072737/50‑1 (MQ=255)            aTGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                 >  1:4297812/1‑50 (MQ=255)             tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgg                                >  1:5087514/1‑50 (MQ=255)              gTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                               >  1:4920410/1‑50 (MQ=255)                 ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcg                            >  1:4713319/1‑50 (MQ=255)                 ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                             <  1:1597556/49‑1 (MQ=255)                     tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGcgcg                          <  1:964190/48‑1 (MQ=255)                     tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGcgcg                          <  1:4732421/48‑1 (MQ=255)                     tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGa                         <  1:4354735/49‑1 (MQ=255)                      cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGaa                        <  1:3695646/49‑1 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                       <  1:4227520/49‑1 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                      >  1:662071/1‑50 (MQ=255)                         tAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCt                    <  1:2360948/50‑1 (MQ=255)                          aaaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                   >  1:931513/1‑50 (MQ=255)                           aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                   >  1:776908/1‑49 (MQ=255)                           aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                   >  1:4931427/1‑49 (MQ=255)                            aCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGc                  >  1:4662971/1‑49 (MQ=255)                            aCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCa                 <  1:910533/50‑1 (MQ=255)                             cGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCaa                >  1:361892/1‑50 (MQ=255)                               tGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCAAACGCTGCAATg              <  1:4312046/50‑1 (MQ=255)                                gcatggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGa             >  1:2817774/1‑50 (MQ=255)                                 catggtatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt            <  1:2312098/50‑1 (MQ=39)                                 catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt            <  1:1889737/50‑1 (MQ=255)                                  atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt           <  1:2590735/50‑1 (MQ=255)                                   tggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttt          >  1:1436531/1‑50 (MQ=255)                                     gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttttt        <  1:3961252/50‑1 (MQ=255)                                     gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttttt        >  1:3554863/1‑50 (MQ=255)                                     gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttttt        >  1:3250786/1‑50 (MQ=255)                                     gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtttt         <  1:1917985/49‑1 (MQ=255)                                      catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc       >  1:2332760/1‑50 (MQ=255)                                      catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc       <  1:4408693/50‑1 (MQ=255)                                      catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc       <  1:177153/50‑1 (MQ=255)                                       atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc       >  1:326072/1‑49 (MQ=255)                                       atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc       >  1:2623294/1‑49 (MQ=255)                                       atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc       >  1:2779748/1‑49 (MQ=255)                                       atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc       >  1:4229237/1‑49 (MQ=255)                                       atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc       >  1:4553548/1‑49 (MQ=255)                                       atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCt      <  1:3835704/50‑1 (MQ=255)                                          cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAtt   <  1:3641116/50‑1 (MQ=255)                                          cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAtt   <  1:4043606/50‑1 (MQ=255)                                          cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAt    <  1:3955677/49‑1 (MQ=255)                                           gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg  <  1:275022/50‑1 (MQ=255)                                           gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg  <  1:949405/50‑1 (MQ=255)                                           gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg  >  1:2593114/1‑50 (MQ=255)                                                |                                             AATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTG  >  NC_000913/118981‑119072

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