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%	307.0 / NA	76	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 (46/30);  total (46/30)

ATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGG  >  NC_000913/118982‑119073                                               |                                              aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                            <  1:4029089/50‑1 (MQ=255) aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                            >  1:6002342/1‑50 (MQ=255) aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                            >  1:2723772/1‑50 (MQ=255)  tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:810758/1‑50 (MQ=255)  tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           <  1:6312874/50‑1 (MQ=255)   aGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:5855945/1‑49 (MQ=255)   aGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:2217786/1‑49 (MQ=255)    gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctctc                                          <  1:5825430/49‑1 (MQ=255)    gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:3617767/1‑48 (MQ=255)    gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:1689747/1‑48 (MQ=255)    gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:5761605/1‑48 (MQ=255)    gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:4500022/1‑48 (MQ=255)     tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                        <  1:4710164/50‑1 (MQ=255)     tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                        >  1:4550143/1‑50 (MQ=255)     tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                        <  1:726512/50‑1 (MQ=255)      ccAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                        >  1:409544/1‑49 (MQ=255)       cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTa                                      >  1:4696973/1‑50 (MQ=255)       cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTa                                      <  1:4473259/50‑1 (MQ=255)       cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTa                                      <  1:3040476/50‑1 (MQ=255)        aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                     >  1:875322/1‑50 (MQ=255)         gTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                     >  1:2416350/1‑49 (MQ=255)         gTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                     >  1:5559102/1‑49 (MQ=255)          tATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                     >  1:4887646/1‑48 (MQ=255)          tATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCa                                   <  1:959773/50‑1 (MQ=255)           aTGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                  >  1:5579241/1‑50 (MQ=255)           aTGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                  >  1:4698711/1‑50 (MQ=255)           aTGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                  >  1:6793059/1‑50 (MQ=255)            tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                  <  1:1223326/49‑1 (MQ=255)             gTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                                >  1:4264239/1‑50 (MQ=255)              ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgggg                               >  1:1602626/1‑50 (MQ=255)              ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                                <  1:2752932/49‑1 (MQ=255)                ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcg                             >  1:1834814/1‑50 (MQ=255)                ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcg                             >  1:7657802/1‑50 (MQ=255)                ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                              <  1:1006955/49‑1 (MQ=255)                ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                              <  1:6982618/49‑1 (MQ=255)                ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                              <  1:7180818/49‑1 (MQ=255)                 cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                            >  1:7067838/1‑50 (MQ=255)                 cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                            <  1:6777811/50‑1 (MQ=255)                 cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                            <  1:5394646/50‑1 (MQ=255)                   ttCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGcgcg                           >  1:2558255/1‑49 (MQ=255)                   ttCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGcgcg                           >  1:7559614/1‑49 (MQ=255)                    tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGcgcg                           <  1:2157593/48‑1 (MQ=255)                    tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGa                          <  1:7804737/49‑1 (MQ=255)                     cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGaa                         <  1:7039765/49‑1 (MQ=255)                     cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGaa                         <  1:2137969/49‑1 (MQ=255)                     cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                        >  1:6091571/1‑50 (MQ=255)                     cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                        >  1:7809153/1‑50 (MQ=255)                      gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                        <  1:3458845/49‑1 (MQ=255)                      gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                        <  1:3200133/49‑1 (MQ=255)                      gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                        <  1:4094202/49‑1 (MQ=255)                      gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                       >  1:1454444/1‑50 (MQ=255)                      gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                       >  1:6537193/1‑50 (MQ=255)                      gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                       >  1:8024541/1‑50 (MQ=255)                       aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                      >  1:7929822/1‑50 (MQ=255)                       aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                      >  1:7028025/1‑50 (MQ=255)                       aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                      <  1:594712/50‑1 (MQ=255)                       aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                      >  1:1135571/1‑50 (MQ=255)                         aaaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                    >  1:4874122/1‑50 (MQ=255)                          aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                    >  1:2213898/1‑49 (MQ=255)                          aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                    <  1:4452354/49‑1 (MQ=255)                          aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGc                   <  1:7880033/50‑1 (MQ=255)                           aCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGc                   >  1:4458437/1‑49 (MQ=255)                           aCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGc                   >  1:2768106/1‑49 (MQ=255)                           aCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGc                   >  1:3100115/1‑49 (MQ=255)                            cGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCaa                 <  1:3596335/50‑1 (MQ=255)                            cGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCaa                 >  1:3186768/1‑50 (MQ=255)                            cGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCaa                 >  1:7103093/1‑50 (MQ=255)                             gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                <  1:7614929/50‑1 (MQ=255)                             gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                >  1:7471161/1‑50 (MQ=255)                             gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                >  1:7971659/1‑50 (MQ=255)                                catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt             <  1:1781902/50‑1 (MQ=255)                                catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt             >  1:2752401/1‑50 (MQ=255)                                catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt             <  1:1457057/50‑1 (MQ=255)                                 atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt            >  1:710059/1‑50 (MQ=255)                                  tggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttt           >  1:6467153/1‑50 (MQ=255)                                   ggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtttt          >  1:4973757/1‑50 (MQ=255)                                    gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtttt          <  1:7561966/49‑1 (MQ=255)                                     catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc        >  1:7490027/1‑50 (MQ=255)                                     catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc        >  1:7810353/1‑50 (MQ=255)                                     catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc        >  1:223510/1‑50 (MQ=255)                                      atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc        >  1:5275860/1‑49 (MQ=255)                                      atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCt       <  1:4703208/50‑1 (MQ=255)                                       tgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTa      <  1:6966786/50‑1 (MQ=255)                                          gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg   >  1:4056283/1‑50 (MQ=255)                                          gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg   >  1:6942761/1‑50 (MQ=255)                                          gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg   >  1:1997409/1‑50 (MQ=255)                                          gggAGCCTCTCAGTATCAGGGGCGCGAACACTGCAATGATTTTTCTATTg   >  1:3969429/1‑50 (MQ=39)                                           ggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTgg  <  1:6812261/50‑1 (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: ‑