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_000913119,0280TC100.0% 324.8 / NA 82V99A (GTC→GCC) ampD1,6‑anhydro‑N‑acetylmuramoyl‑L‑alanine amidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (40/42);  total (40/42)

ATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGG  >  NC_000913/118982‑119074
                                              |                                              
aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                             >  1:3320161/1‑50 (MQ=255)
aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                             >  1:3276203/1‑50 (MQ=255)
aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                             >  1:2444083/1‑50 (MQ=255)
aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                             >  1:6785309/1‑50 (MQ=255)
 tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                            <  1:1001151/50‑1 (MQ=255)
 tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                            <  1:1396758/50‑1 (MQ=255)
 tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                            <  1:4553854/50‑1 (MQ=255)
 tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                            >  1:3570375/1‑50 (MQ=255)
  aGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                            >  1:5534806/1‑49 (MQ=255)
  aGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                            >  1:4100836/1‑49 (MQ=255)
  aGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                            >  1:5403321/1‑49 (MQ=255)
   gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctctc                                           <  1:6691432/49‑1 (MQ=255)
   gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctctc                                           <  1:4002443/49‑1 (MQ=255)
   gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                            >  1:3408287/1‑48 (MQ=255)
    tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                         <  1:4915147/50‑1 (MQ=255)
    tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                         >  1:1986004/1‑50 (MQ=255)
    tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                         >  1:1452908/1‑50 (MQ=255)
      cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGt                                        <  1:3487943/49‑1 (MQ=255)
       aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                      >  1:1619281/1‑50 (MQ=255)
       aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                      >  1:3938414/1‑50 (MQ=255)
       aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                      <  1:398905/50‑1 (MQ=255)
       aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                      <  1:4457322/50‑1 (MQ=255)
       aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                      >  1:6504823/1‑50 (MQ=255)
       aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                      <  1:2429926/50‑1 (MQ=255)
        gTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                      >  1:7023245/1‑49 (MQ=255)
        gTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATc                                     <  1:7535896/50‑1 (MQ=255)
         tATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                      >  1:4612304/1‑48 (MQ=255)
         tATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCa                                    <  1:6263549/50‑1 (MQ=255)
         tATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCa                                    <  1:4922246/50‑1 (MQ=255)
           tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgg                                  <  1:1683675/50‑1 (MQ=255)
           tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                   <  1:731078/49‑1 (MQ=255)
           tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                   <  1:6502414/49‑1 (MQ=255)
            gTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                                 >  1:7012080/1‑50 (MQ=255)
            gTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                                 >  1:6656976/1‑50 (MQ=255)
             ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgggg                                <  1:4519780/50‑1 (MQ=255)
             ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgggg                                >  1:4907790/1‑50 (MQ=255)
              tCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                               >  1:4108839/1‑50 (MQ=255)
              tCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                               >  1:2756871/1‑50 (MQ=255)
               ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                               <  1:3471253/49‑1 (MQ=255)
               ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                               <  1:2556487/49‑1 (MQ=255)
                cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                             <  1:3873036/50‑1 (MQ=255)
                cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                             >  1:5610656/1‑50 (MQ=255)
                cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                             <  1:2949481/50‑1 (MQ=255)
                cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                             <  1:4858297/50‑1 (MQ=255)
                cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                             <  1:5918372/50‑1 (MQ=255)
                 tttCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                             >  1:5006741/1‑49 (MQ=255)
                  ttCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGcgcg                            >  1:3464599/1‑49 (MQ=255)
                   tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGa                           <  1:6262135/49‑1 (MQ=255)
                    cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                         >  1:5184990/1‑50 (MQ=255)
                     gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                         <  1:913199/49‑1 (MQ=255)
                     gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                        <  1:2292220/50‑1 (MQ=255)
                      aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                       >  1:47405/1‑50 (MQ=255)
                       tAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCt                      <  1:5795376/50‑1 (MQ=255)
                       tAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCt                      <  1:4826425/50‑1 (MQ=255)
                        aaaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                     >  1:6045888/1‑50 (MQ=255)
                         aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                     >  1:2257579/1‑49 (MQ=255)
                         aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGc                    <  1:2067256/50‑1 (MQ=255)
                         aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGc                    <  1:3778228/50‑1 (MQ=255)
                          aCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGc                    >  1:6227449/1‑49 (MQ=255)
                          aCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGc                    >  1:2825690/1‑49 (MQ=255)
                          aCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCa                   <  1:4638098/50‑1 (MQ=255)
                           cGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCaa                  >  1:2476883/1‑50 (MQ=255)
                            gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                 >  1:2821290/1‑50 (MQ=255)
                            gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                 <  1:5153479/50‑1 (MQ=255)
                            gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                 >  1:5002463/1‑50 (MQ=255)
                             tGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                 >  1:5073225/1‑49 (MQ=255)
                             tGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATg                <  1:3182925/50‑1 (MQ=255)
                              gcatggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGa               >  1:2418775/1‑50 (MQ=255)
                               catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt              >  1:2797098/1‑50 (MQ=255)
                               catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt              >  1:304806/1‑50 (MQ=255)
                               catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt              >  1:3301426/1‑50 (MQ=255)
                               catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt              >  1:3413777/1‑50 (MQ=255)
                                atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt             <  1:5284895/50‑1 (MQ=255)
                                atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt             >  1:4810371/1‑50 (MQ=255)
                                 tggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttt            >  1:6458127/1‑50 (MQ=255)
                                   gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttttt          >  1:4089641/1‑50 (MQ=255)
                                   gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttttt          <  1:4382166/50‑1 (MQ=255)
                                   gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtttt           <  1:4431136/49‑1 (MQ=255)
                                    catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc         <  1:2032605/50‑1 (MQ=255)
                                    catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc         >  1:5701870/1‑50 (MQ=255)
                                    catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc         <  1:4372256/50‑1 (MQ=255)
                                     atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc         >  1:6492833/1‑49 (MQ=255)
                                     atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc         >  1:3628306/1‑49 (MQ=255)
                                     atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCt        <  1:1185151/50‑1 (MQ=255)
                                      tgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTa       <  1:5410480/50‑1 (MQ=255)
                                       gCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAt      >  1:4726192/1‑50 (MQ=255)
                                        cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAt      <  1:2636625/49‑1 (MQ=255)
                                        cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAt      <  1:3781380/49‑1 (MQ=255)
                                        cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAt      <  1:689068/49‑1 (MQ=255)
                                         gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg    >  1:513619/1‑50 (MQ=255)
                                         gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg    <  1:4978076/50‑1 (MQ=255)
                                          ggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTgg   >  1:2080809/1‑50 (MQ=255)
                                          ggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTgg   <  1:1767/50‑1 (MQ=255)
                                           gAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTggg  <  1:6971464/50‑1 (MQ=255)
                                           gAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTggg  <  1:7394064/50‑1 (MQ=255)
                                           gAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTgg   >  1:7468970/1‑49 (MQ=255)
                                              |                                              
ATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGG  >  NC_000913/118982‑119074

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: