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%	277.5 / NA	70	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 (50/20);  total (50/20)

ATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGG  >  NC_000913/118982‑119073                                               |                                              aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                            >  1:1239651/1‑50 (MQ=255) aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                            >  1:4331454/1‑50 (MQ=255)  tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:7549271/1‑50 (MQ=255)   aGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:2243323/1‑49 (MQ=255)   aGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:183042/1‑49 (MQ=255)    gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctctc                                          <  1:7405797/49‑1 (MQ=255)    gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctctc                                          <  1:3180268/49‑1 (MQ=255)    gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctctc                                          <  1:7638063/49‑1 (MQ=255)    gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:6888265/1‑48 (MQ=255)    gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:7270385/1‑48 (MQ=255)    gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:544742/1‑48 (MQ=255)    gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:1284001/1‑48 (MQ=255)    gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:1365725/1‑48 (MQ=255)     tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                        >  1:17359/1‑50 (MQ=255)     tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                        >  1:7226622/1‑50 (MQ=255)     tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                        >  1:5631018/1‑50 (MQ=255)      ccAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                        >  1:5753590/1‑49 (MQ=255)      ccAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                        >  1:3759504/1‑49 (MQ=255)        aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                     >  1:7928045/1‑50 (MQ=255)        aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                     >  1:3311701/1‑50 (MQ=255)         gTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                     >  1:158090/1‑49 (MQ=255)          tATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                     >  1:1439746/1‑48 (MQ=255)          tATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCa                                   <  1:1582071/50‑1 (MQ=255)          tATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCa                                   <  1:6125217/50‑1 (MQ=255)           aTGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                  >  1:7465785/1‑50 (MQ=255)            tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgg                                 >  1:1613348/1‑50 (MQ=255)            tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgg                                 >  1:7342447/1‑50 (MQ=255)            tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                  <  1:5477379/49‑1 (MQ=255)            tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                  <  1:7661984/49‑1 (MQ=255)             gTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                                >  1:269697/1‑50 (MQ=255)              ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgggg                               >  1:8278989/1‑50 (MQ=255)               tCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                              >  1:1204483/1‑50 (MQ=255)               tCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                              >  1:6275540/1‑50 (MQ=255)               tCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                              >  1:7853610/1‑50 (MQ=255)               tCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                              >  1:188380/1‑50 (MQ=255)               tCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGCATCAGGGgc                              >  1:8224474/1‑50 (MQ=255)                ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcg                             >  1:632754/1‑50 (MQ=255)                ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                              <  1:5493024/49‑1 (MQ=255)                ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                              <  1:5894817/49‑1 (MQ=255)                ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                              <  1:6680209/49‑1 (MQ=255)                ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                              <  1:2387052/49‑1 (MQ=255)                 cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                            >  1:4117830/1‑50 (MQ=255)                 cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcg                             <  1:1651590/49‑1 (MQ=255)                  tttCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                            >  1:7884443/1‑49 (MQ=255)                   ttCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGcgcg                           >  1:4185736/1‑49 (MQ=255)                   ttCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGcgcg                           >  1:3148158/1‑49 (MQ=255)                   ttCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGcgcg                           >  1:8256908/1‑49 (MQ=255)                    tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGa                          <  1:4456851/49‑1 (MQ=255)                     cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                        >  1:4915156/1‑50 (MQ=255)                      gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                        <  1:6689175/49‑1 (MQ=255)                      gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                       >  1:5368671/1‑50 (MQ=255)                      gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                       <  1:2233299/50‑1 (MQ=255)                       aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                      >  1:198070/1‑50 (MQ=255)                       aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                      >  1:7748557/1‑50 (MQ=255)                         aaaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                    >  1:100533/1‑50 (MQ=255)                         aaaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                    >  1:4470544/1‑50 (MQ=255)                          aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                    >  1:7537694/1‑49 (MQ=255)                          aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                    >  1:2084718/1‑49 (MQ=255)                          aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                    >  1:3352457/1‑49 (MQ=255)                          aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGAGCGAACGCTg                    <  1:3474605/49‑1 (MQ=255)                           aCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGc                   >  1:8166651/1‑49 (MQ=255)                            cGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCaa                 >  1:6703459/1‑50 (MQ=255)                             gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                >  1:5390534/1‑50 (MQ=255)                             gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                <  1:812005/50‑1 (MQ=255)                             gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                >  1:3159528/1‑50 (MQ=255)                               gcatggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGa              >  1:2748502/1‑50 (MQ=255)                                catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGa              <  1:898882/49‑1 (MQ=255)                                catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt             >  1:7952367/1‑50 (MQ=255)                                 atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt            >  1:623822/1‑50 (MQ=255)                                 atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt            >  1:3019540/1‑50 (MQ=255)                                 atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt            <  1:5697304/50‑1 (MQ=255)                                    gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtttt          <  1:3985738/49‑1 (MQ=255)                                     catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc        >  1:5174221/1‑50 (MQ=255)                                     catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc        >  1:1950877/1‑50 (MQ=255)                                     catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc        <  1:1349896/50‑1 (MQ=255)                                     catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc        >  1:1277289/1‑50 (MQ=255)                                      atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc        >  1:3258119/1‑49 (MQ=255)                                       tgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTa      >  1:7301543/1‑50 (MQ=255)                                        gCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAt     >  1:5535049/1‑50 (MQ=255)                                         cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAtt    <  1:2930438/50‑1 (MQ=255)                                         cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAtt    <  1:9281/50‑1 (MQ=255)                                            gAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTgg  >  1:6637750/1‑49 (MQ=255)                                            gAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTgg  >  1:658509/1‑49 (MQ=255)                                               |                                              ATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGG  >  NC_000913/118982‑119073

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