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%	240.2 / NA	63	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 (39/24);  total (39/24)

AATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGG  >  NC_000913/118981‑119074                                                |                                               aaTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcct                                              >  1:3541492/1‑50 (MQ=255) aaTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcct                                              <  1:4544369/50‑1 (MQ=255)  aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                             >  1:2573074/1‑50 (MQ=255)  aTAGTCCAGTATGTTCCTTTCGATAAACGCGCATGGCATGCGGGAGcctc                                             >  1:546542/1‑50 (MQ=39)   tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGActctct                                            >  1:63424/1‑50 (MQ=255)   tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                            <  1:711264/50‑1 (MQ=255)   tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                            <  1:2493479/50‑1 (MQ=255)   tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                            >  1:2017731/1‑50 (MQ=255)   tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                            >  1:3740316/1‑50 (MQ=255)     gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctctc                                           <  1:2929540/49‑1 (MQ=255)     gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctctc                                           <  1:4844090/49‑1 (MQ=255)     gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                            >  1:1243059/1‑48 (MQ=255)      tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                         >  1:5632130/1‑50 (MQ=255)      tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                         >  1:5301915/1‑50 (MQ=255)      tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                         >  1:4589116/1‑50 (MQ=255)       ccAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                         >  1:394541/1‑49 (MQ=255)        cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTa                                       <  1:5600965/50‑1 (MQ=255)         aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                      >  1:1912844/1‑50 (MQ=255)            aTGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                   >  1:3516963/1‑50 (MQ=255)             tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgg                                  <  1:455702/50‑1 (MQ=255)             tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                   <  1:269723/49‑1 (MQ=255)             tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                   <  1:4377746/49‑1 (MQ=255)                 ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcg                              >  1:7062231/1‑50 (MQ=255)                 ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                               <  1:6580686/49‑1 (MQ=255)                  cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                             <  1:2812099/50‑1 (MQ=255)                   tttCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                             >  1:2569584/1‑49 (MQ=255)                    ttCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGcgcg                            >  1:6502856/1‑49 (MQ=255)                     tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGcgcg                            <  1:408664/48‑1 (MQ=255)                     tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGa                           <  1:6344149/49‑1 (MQ=255)                     tCGATAAACGTGCATGGCATACGGGAGCCTCTCAGTATCAGGGGCGCGaa                          >  1:3863151/1‑50 (MQ=39)                      cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGaa                          <  1:5107369/49‑1 (MQ=255)                      cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGaa                          <  1:124253/49‑1 (MQ=255)                      cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                         >  1:5605317/1‑50 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                         <  1:3852490/49‑1 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                         <  1:2903066/49‑1 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                        >  1:710924/1‑50 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                        >  1:5625607/1‑50 (MQ=255)                       gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                        <  1:2555019/50‑1 (MQ=255)                        aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                       >  1:1810926/1‑50 (MQ=255)                          aaaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                     >  1:3722848/1‑50 (MQ=255)                          aaaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                     >  1:2572883/1‑50 (MQ=255)                          aaaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                     >  1:4471151/1‑50 (MQ=255)                           aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                     >  1:5781707/1‑49 (MQ=255)                           aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGc                    <  1:3465963/50‑1 (MQ=255)                           aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGc                    <  1:1704986/50‑1 (MQ=255)                             cGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCaa                  >  1:2955594/1‑50 (MQ=255)                              gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                 >  1:3104080/1‑50 (MQ=255)                              gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                 >  1:6116150/1‑50 (MQ=255)                              gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                 >  1:6162657/1‑50 (MQ=255)                               tGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                 >  1:4566843/1‑49 (MQ=255)                               tGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATg                <  1:6769025/50‑1 (MQ=255)                               tGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATg                <  1:3464862/50‑1 (MQ=255)                                gcatggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGa               >  1:2089505/1‑50 (MQ=255)                                  atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt             >  1:4800781/1‑50 (MQ=255)                                  atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt             >  1:1516622/1‑50 (MQ=255)                                  atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt             >  1:5095566/1‑50 (MQ=255)                                  atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt             <  1:2544148/50‑1 (MQ=255)                                  atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt             >  1:603543/1‑50 (MQ=255)                                    ggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtttt           >  1:3505452/1‑50 (MQ=255)                                      catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc         >  1:2029794/1‑50 (MQ=255)                                      catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc         >  1:4955256/1‑50 (MQ=255)                                      catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc         <  1:3458651/50‑1 (MQ=255)                                      catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc         >  1:4545063/1‑50 (MQ=255)                                       atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc         >  1:5610480/1‑49 (MQ=255)                                       atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc         >  1:2432122/1‑49 (MQ=255)                                       atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCt        <  1:3119725/50‑1 (MQ=255)                                         gCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAt      >  1:3533657/1‑50 (MQ=255)                                          cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAtt     <  1:3801984/50‑1 (MQ=255)                                           gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg    >  1:5978184/1‑50 (MQ=255)                                           gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg    >  1:539179/1‑50 (MQ=255)                                           gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg    <  1:6628310/50‑1 (MQ=255)                                           gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg    <  1:1331928/50‑1 (MQ=255)                                           gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg    >  1:3711902/1‑50 (MQ=255)                                            ggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTgg   >  1:1231436/1‑50 (MQ=255)                                             gAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTggg  <  1:4635731/50‑1 (MQ=255)                                                |                                               AATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGG  >  NC_000913/118981‑119074

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