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%	317.6 / NA	80	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 (49/31);  total (49/31)

AATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGG  >  NC_000913/118981‑119073                                                |                                              aaTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcct                                             <  1:4861747/50‑1 (MQ=255) aaTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcct                                             <  1:1467970/50‑1 (MQ=255)  aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                            <  1:1224111/50‑1 (MQ=255)  aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                            >  1:6609123/1‑50 (MQ=255)   tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:1626584/1‑50 (MQ=255)    aGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:2232440/1‑49 (MQ=255)    aGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:3949353/1‑49 (MQ=255)     gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctctc                                          <  1:2114823/49‑1 (MQ=255)     gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctctc                                          <  1:5160886/49‑1 (MQ=255)     gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctctc                                          <  1:6015830/49‑1 (MQ=255)     gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:3958664/1‑48 (MQ=255)     gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCa                                         <  1:5752620/50‑1 (MQ=255)     gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCa                                         <  1:700228/50‑1 (MQ=255)       ccAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                        >  1:5131302/1‑49 (MQ=255)        cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTa                                      >  1:5256897/1‑50 (MQ=255)        cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTa                                      >  1:3039820/1‑50 (MQ=255)         aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                     >  1:4995903/1‑50 (MQ=255)         aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                     <  1:1337530/50‑1 (MQ=255)          gTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                     >  1:5802880/1‑49 (MQ=255)          gTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATc                                    <  1:1159726/50‑1 (MQ=255)            aTGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                  >  1:1973473/1‑50 (MQ=255)            aTGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                  >  1:2272734/1‑50 (MQ=255)             tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                  <  1:3723288/49‑1 (MQ=255)              gTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                                >  1:938392/1‑50 (MQ=255)              gTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                                >  1:4222429/1‑50 (MQ=255)              gTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                                >  1:339954/1‑50 (MQ=255)               ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgggg                               >  1:1792753/1‑50 (MQ=255)               ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                                <  1:6175256/49‑1 (MQ=255)               ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                                <  1:1409522/49‑1 (MQ=255)               ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                                <  1:1862057/49‑1 (MQ=255)                 ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcg                             >  1:1888706/1‑50 (MQ=255)                 ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                              <  1:3370114/49‑1 (MQ=255)                  cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                            <  1:1835480/50‑1 (MQ=255)                  cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                            >  1:6345549/1‑50 (MQ=255)                     tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGcgcg                           <  1:194346/48‑1 (MQ=255)                     tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGcgcg                           <  1:5648330/48‑1 (MQ=255)                     tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGaa                         >  1:6404699/1‑50 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                        <  1:2394811/49‑1 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                        <  1:6028822/49‑1 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                       <  1:2334014/50‑1 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                       >  1:2890440/1‑50 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                       >  1:1480163/1‑50 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                       >  1:4729532/1‑50 (MQ=255)                        aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                      >  1:5211035/1‑50 (MQ=255)                        aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                      <  1:1346630/50‑1 (MQ=255)                         tAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCt                     >  1:936662/1‑50 (MQ=255)                         tAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCt                     >  1:5926874/1‑50 (MQ=255)                          aaaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                    >  1:2812678/1‑50 (MQ=255)                          aaaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                    >  1:2251838/1‑50 (MQ=255)                           aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                    <  1:2974731/49‑1 (MQ=255)                           aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                    >  1:1988494/1‑49 (MQ=255)                           aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                    >  1:6351010/1‑49 (MQ=255)                            aCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCa                  <  1:3790639/50‑1 (MQ=255)                             cGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCaa                 >  1:1074110/1‑50 (MQ=255)                              gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                >  1:2590323/1‑50 (MQ=255)                              gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                >  1:2109457/1‑50 (MQ=255)                               tGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                >  1:2348836/1‑49 (MQ=255)                               tGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATg               <  1:5806346/50‑1 (MQ=255)                               tGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATg               <  1:3289863/50‑1 (MQ=255)                                gcatggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGa              >  1:6399867/1‑50 (MQ=255)                                 catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt             <  1:5128195/50‑1 (MQ=255)                                 catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt             >  1:6736863/1‑50 (MQ=255)                                  atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt            <  1:1696096/50‑1 (MQ=255)                                  atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt            >  1:4819068/1‑50 (MQ=255)                                    ggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtttt          >  1:3471046/1‑50 (MQ=255)                                     gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttttt         >  1:4224350/1‑50 (MQ=255)                                     gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttttt         >  1:3497299/1‑50 (MQ=255)                                     gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttttt         >  1:3591751/1‑50 (MQ=255)                                     gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttttt         >  1:2007666/1‑50 (MQ=255)                                     gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttttt         <  1:3272275/50‑1 (MQ=255)                                     gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttttt         <  1:5375918/50‑1 (MQ=255)                                     gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtttt          <  1:4900523/49‑1 (MQ=255)                                      catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc        >  1:3967355/1‑50 (MQ=255)                                      catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc        >  1:6036989/1‑50 (MQ=255)                                       atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc        >  1:755287/1‑49 (MQ=255)                                       atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc        >  1:3099638/1‑49 (MQ=255)                                       atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc        >  1:1571134/1‑49 (MQ=255)                                        tgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTa      <  1:5843341/50‑1 (MQ=255)                                         gCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAt     >  1:4752916/1‑50 (MQ=255)                                         gCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAt     >  1:1755233/1‑50 (MQ=255)                                          cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAtt    <  1:6742441/50‑1 (MQ=255)                                          cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAtt    <  1:3848249/50‑1 (MQ=255)                                          cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAtt    <  1:6686450/50‑1 (MQ=255)                                          cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAtt    >  1:6130983/1‑50 (MQ=255)                                          cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAt     <  1:2193658/49‑1 (MQ=255)                                           gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg   >  1:4559237/1‑50 (MQ=255)                                           gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg   >  1:1198238/1‑50 (MQ=255)                                           gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg   <  1:148425/50‑1 (MQ=255)                                             gAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTgg  >  1:6419295/1‑49 (MQ=255)                                             gAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTgg  >  1:3638572/1‑49 (MQ=255)                                             gAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTgg  >  1:5631100/1‑49 (MQ=255)                                                |                                              AATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGG  >  NC_000913/118981‑119073

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: ‑