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%	308.7 / NA	78	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 (47/31);  total (47/31)

ATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGG  >  NC_000913/118982‑119074                                               |                                               aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                             <  1:773436/50‑1 (MQ=255) aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                             <  1:7090638/50‑1 (MQ=255) aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                             >  1:2353025/1‑50 (MQ=255) aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                             >  1:5636849/1‑50 (MQ=255)  tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                            <  1:7249291/50‑1 (MQ=255)  tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                            <  1:4610336/50‑1 (MQ=255)  tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                            >  1:6815890/1‑50 (MQ=255)   aGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                            >  1:7413538/1‑49 (MQ=255)   aGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                            >  1:6010838/1‑49 (MQ=255)    gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctctc                                           <  1:3258848/49‑1 (MQ=255)    gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                            >  1:1054538/1‑48 (MQ=255)     tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                         >  1:66696/1‑50 (MQ=255)     tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                         >  1:6845573/1‑50 (MQ=255)      ccAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                         >  1:1925703/1‑49 (MQ=255)      ccAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                         >  1:3768747/1‑49 (MQ=255)       cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGt                                        <  1:4832863/49‑1 (MQ=255)       cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGt                                        <  1:4540422/49‑1 (MQ=255)       cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTa                                       >  1:6617599/1‑50 (MQ=255)       cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTa                                       >  1:2238463/1‑50 (MQ=255)       cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTa                                       >  1:6897598/1‑50 (MQ=255)        aGTATGTTCCTTTCGATCAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                      >  1:863748/1‑50 (MQ=255)        aGTATGTTCCTTTCGATAAACGTGCGTGGCATGCGGGAGCCTCTCAGTAt                                      <  1:1097793/50‑1 (MQ=39)        aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                      >  1:7185152/1‑50 (MQ=255)          tATGTTCCTTTCGATAAACGTGCATGGCATGCGGGATCCTCTCAGTATcc                                    >  1:4149076/1‑49 (MQ=38)          tATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCa                                    <  1:1280794/50‑1 (MQ=255)           aTGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                   >  1:2614376/1‑50 (MQ=255)            tGTTCCTTTCGATAAACGTGCATGGCATGCGTGAGCCTCTCAGTATCAgg                                  >  1:5236700/1‑50 (MQ=39)            tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgg                                  >  1:2568240/1‑50 (MQ=255)              ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgggg                                >  1:3848095/1‑50 (MQ=255)              ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgggg                                <  1:2993228/50‑1 (MQ=255)              ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgggg                                >  1:7401738/1‑50 (MQ=255)              ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                                 <  1:2501052/49‑1 (MQ=255)               tCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                               >  1:4888340/1‑50 (MQ=255)                 cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                             >  1:1877952/1‑50 (MQ=255)                 cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                             >  1:6231794/1‑50 (MQ=255)                 cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                             >  1:1901328/1‑50 (MQ=255)                 cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                             <  1:5073666/50‑1 (MQ=255)                 cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                             <  1:5479953/50‑1 (MQ=255)                 cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                             >  1:2474870/1‑50 (MQ=255)                    tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGcgcg                            <  1:2637190/48‑1 (MQ=255)                    tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGaa                          >  1:4316106/1‑50 (MQ=255)                    tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGaa                          >  1:1063782/1‑50 (MQ=255)                     cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                         >  1:1594766/1‑50 (MQ=255)                      gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                         <  1:664067/49‑1 (MQ=255)                      gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                         <  1:5881952/49‑1 (MQ=255)                      gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                        >  1:1695864/1‑50 (MQ=255)                      gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                        >  1:1597577/1‑50 (MQ=255)                      gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                        <  1:7254829/50‑1 (MQ=255)                       aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                       >  1:1703239/1‑50 (MQ=255)                       aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                       <  1:6064085/50‑1 (MQ=255)                       aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                       <  1:5026750/50‑1 (MQ=255)                         aaaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                     >  1:2648490/1‑50 (MQ=255)                         aaaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                     >  1:4306195/1‑50 (MQ=255)                          aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                     >  1:3710936/1‑49 (MQ=255)                           aCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGc                    >  1:6511882/1‑49 (MQ=255)                           aCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGc                    >  1:2366279/1‑49 (MQ=255)                           aCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCa                   <  1:2604615/50‑1 (MQ=255)                            cGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCaa                  <  1:3713631/50‑1 (MQ=255)                            cGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCaa                  >  1:2301509/1‑50 (MQ=255)                            cGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCaa                  <  1:3595899/50‑1 (MQ=255)                             gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                 <  1:4105428/50‑1 (MQ=255)                             gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                 <  1:3276286/50‑1 (MQ=255)                              tGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                 >  1:6672877/1‑49 (MQ=255)                              tGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATg                <  1:1945585/50‑1 (MQ=255)                               gcatggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGa               >  1:6731875/1‑50 (MQ=255)                               gcatggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGa               >  1:4220320/1‑50 (MQ=255)                               gcatggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGa               >  1:7173560/1‑50 (MQ=255)                                catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt              >  1:5612449/1‑50 (MQ=255)                                catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt              >  1:3653050/1‑50 (MQ=255)                                 atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt             <  1:412235/50‑1 (MQ=255)                                 atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt             >  1:472921/1‑50 (MQ=255)                                 atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt             <  1:2160074/50‑1 (MQ=255)                                  tggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttt            >  1:5601508/1‑50 (MQ=255)                                   ggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtttt           >  1:522522/1‑50 (MQ=255)                                     catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc         <  1:5508295/50‑1 (MQ=255)                                      atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc         >  1:3029067/1‑49 (MQ=255)                                      atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc         >  1:5283079/1‑49 (MQ=255)                                      atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCt        <  1:6610519/50‑1 (MQ=255)                                         cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAtt     >  1:6658526/1‑50 (MQ=255)                                         cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAtt     <  1:5716024/50‑1 (MQ=255)                                         cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAt      <  1:5962128/49‑1 (MQ=255)                                          gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg    <  1:4657433/50‑1 (MQ=255)                                          gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg    <  1:1107187/50‑1 (MQ=255)                                          gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg    >  1:2530546/1‑50 (MQ=255)                                          gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg    >  1:5653133/1‑50 (MQ=255)                                           ggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTgg   <  1:7035420/50‑1 (MQ=255)                                           ggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTgg   <  1:1066641/50‑1 (MQ=255)                                            gAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTggg  <  1:2944051/50‑1 (MQ=255)                                            gAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTgg   >  1:6492067/1‑49 (MQ=255)                                               |                                               ATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGG  >  NC_000913/118982‑119074

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