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%	315.4 / NA	79	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 (43/36);  total (43/36)

AATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGG  >  NC_000913/118981‑119073                                                |                                              aaTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcct                                             >  1:647428/1‑50 (MQ=255) aaTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcct                                             <  1:6275779/50‑1 (MQ=255)  aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                            >  1:2177318/1‑50 (MQ=255)  aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                            >  1:5359533/1‑50 (MQ=255)   tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:2497271/1‑50 (MQ=255)    aGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:988796/1‑49 (MQ=255)    aGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:5541907/1‑49 (MQ=255)     gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCa                                         <  1:2708964/50‑1 (MQ=255)     gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCa                                         <  1:4022287/50‑1 (MQ=255)      tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                        <  1:6351762/50‑1 (MQ=255)        cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGt                                       <  1:5684327/49‑1 (MQ=255)        cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGt                                       <  1:6675252/49‑1 (MQ=255)        cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTa                                      >  1:4786365/1‑50 (MQ=255)         aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                     >  1:4597267/1‑50 (MQ=255)          gTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                     >  1:175859/1‑49 (MQ=255)           tATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                     >  1:6182692/1‑48 (MQ=255)           tATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCa                                   <  1:5057176/50‑1 (MQ=255)           tATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCa                                   <  1:1038651/50‑1 (MQ=255)           tATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCa                                   <  1:951409/50‑1 (MQ=255)             tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgg                                 <  1:6017945/50‑1 (MQ=255)             tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgg                                 >  1:4108097/1‑50 (MQ=255)             tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgg                                 <  1:3617710/50‑1 (MQ=255)             tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                  <  1:44175/49‑1 (MQ=255)              gTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                                >  1:620188/1‑50 (MQ=255)              gTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                                >  1:4082475/1‑50 (MQ=255)              gTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                                >  1:1080698/1‑50 (MQ=255)               ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgggg                               <  1:6651277/50‑1 (MQ=255)               ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgggg                               >  1:1818603/1‑50 (MQ=255)               ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                                <  1:4125382/49‑1 (MQ=255)               ctCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgggg                               >  1:1441181/2‑50 (MQ=255)                   tttCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                            >  1:3905269/1‑49 (MQ=255)                   tttCGATAAACGTGCATGGCATGCGGGAGCCTCCCAGTATCAGGGgcgc                            >  1:2292191/1‑49 (MQ=255)                    ttCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGcgcg                           >  1:1738052/1‑49 (MQ=255)                    ttCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGcgcg                           >  1:2312904/1‑49 (MQ=255)                    ttCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGcgcg                           >  1:3726006/1‑49 (MQ=255)                     tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGaa                         >  1:5038649/1‑50 (MQ=255)                     tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGaa                         >  1:3342362/1‑50 (MQ=255)                     tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGa                          <  1:4145717/49‑1 (MQ=255)                      cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                        >  1:6367260/1‑50 (MQ=255)                      cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                        >  1:6035869/1‑50 (MQ=255)                      cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                        >  1:3126639/1‑50 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                        <  1:407709/49‑1 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                        <  1:1553916/49‑1 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                        <  1:4642360/49‑1 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                       >  1:1214668/1‑50 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                       <  1:931849/50‑1 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                       >  1:1321034/1‑50 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                       <  1:1345383/50‑1 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                       >  1:6113183/1‑50 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                       >  1:1755579/1‑50 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                       <  1:5487511/50‑1 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGAGCGCGAAc                        <  1:2832423/49‑1 (MQ=255)                        aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                      <  1:6540483/50‑1 (MQ=255)                        aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                      <  1:2205675/50‑1 (MQ=255)                        aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                      >  1:5280596/1‑50 (MQ=255)                         tAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCt                     >  1:3382462/1‑50 (MQ=255)                         tAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCt                     >  1:2091740/1‑50 (MQ=255)                           aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                    >  1:3669747/1‑49 (MQ=255)                           aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                    >  1:2940830/1‑49 (MQ=255)                           aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                    >  1:6323291/1‑49 (MQ=255)                           aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                    <  1:437816/49‑1 (MQ=255)                           aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGc                   <  1:5140751/50‑1 (MQ=255)                              gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                <  1:3985047/50‑1 (MQ=255)                              gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                <  1:6122134/50‑1 (MQ=255)                              gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                <  1:2519512/50‑1 (MQ=255)                               tGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                >  1:2429107/1‑49 (MQ=255)                                gcatggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGa              >  1:4723677/1‑50 (MQ=255)                                gcatggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGa              >  1:3726461/1‑50 (MQ=255)                                gcatggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGa              >  1:1439415/1‑50 (MQ=255)                                 catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt             >  1:6392104/1‑50 (MQ=255)                                  atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt            <  1:5902764/50‑1 (MQ=255)                                  atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt            >  1:686942/1‑50 (MQ=255)                                  atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt            >  1:6472484/1‑50 (MQ=255)                                    ggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtttt          >  1:1836748/1‑50 (MQ=255)                                     gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttttt         <  1:3165629/50‑1 (MQ=255)                                     gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttttt         <  1:1870828/50‑1 (MQ=255)                                      catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc        >  1:4039350/1‑50 (MQ=255)                                      catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc        <  1:6277576/50‑1 (MQ=255)                                      catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc        <  1:5795708/50‑1 (MQ=255)                                      catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc        <  1:1748205/50‑1 (MQ=255)                                       atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc        >  1:4859453/1‑49 (MQ=255)                                        tgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTa      <  1:1290512/50‑1 (MQ=255)                                          cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAt     <  1:3081141/49‑1 (MQ=255)                                           gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg   >  1:215279/1‑50 (MQ=255)                                           gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg   >  1:4314190/1‑50 (MQ=255)                                             gAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTgg  >  1:3002623/1‑49 (MQ=255)                                                |                                              AATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGG  >  NC_000913/118981‑119073

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