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,0280TC93.9% 461.4 / 12.7 132V99A (GTC→GCC) ampD1,6‑anhydro‑N‑acetylmuramoyl‑L‑alanine amidase
Reads supporting (aligned to +/- strand):  ref base T (3/5);  new base C (65/59);  total (68/64)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 146 total aligned reads displayed.

AATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGAT  >  NC_000913/118981‑119076
                                               |                                                
 aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGtctc                                               <  1:2225/50‑1 (MQ=255)
 aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                               <  1:4961844/50‑1 (MQ=255)
 aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                               <  2:2134644/50‑1 (MQ=255)
 aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                               >  2:1781852/1‑50 (MQ=255)
 aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                               >  2:1578651/1‑50 (MQ=255)
 aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                               >  1:2349046/1‑50 (MQ=255)
  tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                              <  2:679314/50‑1 (MQ=255)
  tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                              <  1:3930493/50‑1 (MQ=255)
  tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                              <  1:935981/50‑1 (MQ=255)
    gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCa                                            <  2:1969490/50‑1 (MQ=255)
    gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCa                                            <  2:3966020/50‑1 (MQ=255)
    gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCa                                            <  2:112005/50‑1 (MQ=255)
     tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                           >  2:3248078/1‑50 (MQ=255)
     tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                           >  2:3604089/1‑50 (MQ=255)
     tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                           >  1:4903378/1‑50 (MQ=255)
     tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                           >  2:321633/1‑50 (MQ=255)
     tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                           >  2:2140887/1‑50 (MQ=255)
     tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTAAg                                           <  2:4015754/50‑1 (MQ=39)
      ccAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                           >  2:1774878/1‑49 (MQ=255)
      ccAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                           >  1:306055/1‑49 (MQ=255)
       cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGt                                          <  2:4099269/49‑1 (MQ=255)
       cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTa                                         <  1:4725052/50‑1 (MQ=255)
       cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTa                                         >  1:68303/1‑50 (MQ=255)
        aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                        >  1:3240325/1‑50 (MQ=255)
        aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                        <  1:1752788/50‑1 (MQ=255)
        aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                        >  2:5057666/1‑50 (MQ=255)
         gTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                        >  1:1224074/1‑49 (MQ=255)
         gTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATc                                       <  1:4907314/50‑1 (MQ=255)
          tATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                        >  2:2754615/1‑48 (MQ=255)
            tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgg                                    <  2:3331317/50‑1 (MQ=255)
            tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgg                                    <  1:68719/50‑1 (MQ=255)
            tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                     <  1:3730459/49‑1 (MQ=255)
             gTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                                   >  2:4717739/1‑50 (MQ=255)
              ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgggg                                  >  1:2191419/1‑50 (MQ=255)
              ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgggg                                  <  2:1158255/50‑1 (MQ=255)
               tCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                                 >  2:5025889/1‑50 (MQ=255)
                ccTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGgc                                 <  2:3619754/49‑1 (MQ=255)
                ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcg                                >  1:245872/1‑50 (MQ=255)
                ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                                 <  2:2396714/49‑1 (MQ=255)
                 cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                               <  1:3356879/50‑1 (MQ=255)
                 cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                               >  1:1772660/1‑50 (MQ=255)
                 cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                               <  2:4295277/50‑1 (MQ=255)
                  tttCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                               >  2:3815184/1‑49 (MQ=255)
                  tttCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                               >  1:5251329/1‑49 (MQ=255)
                   ttCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGcgcg                              >  1:2254949/1‑49 (MQ=255)
                   ttCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGcgcg                              >  1:1958871/1‑49 (MQ=255)
                    tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGaa                            >  2:733918/1‑50 (MQ=255)
                    tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGa                             <  1:2275590/49‑1 (MQ=255)
                    tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGa                             <  1:22162/49‑1 (MQ=255)
                    tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGa                             <  2:4859217/49‑1 (MQ=255)
                     cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                           >  2:4979550/1‑50 (MQ=255)
                     cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                           >  2:2827042/1‑50 (MQ=255)
                     cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                           >  1:862238/1‑50 (MQ=255)
                     cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                           >  2:5350511/1‑50 (MQ=255)
                      gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                           <  2:309759/49‑1 (MQ=255)
                      gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                           <  2:2081418/49‑1 (MQ=255)
                      gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                           <  2:802175/49‑1 (MQ=255)
                       aTAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGc                         >  1:1330045/1‑50 (MQ=255)
                       aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                         >  1:2925359/1‑50 (MQ=255)
                        tAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCt                        >  2:5026213/1‑50 (MQ=255)
                        tAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCt                        >  2:915377/1‑50 (MQ=255)
                        tAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCt                        <  1:5231140/50‑1 (MQ=255)
                        tAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCt                        <  2:1801622/50‑1 (MQ=255)
                        tAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCt                        >  1:5350171/1‑50 (MQ=255)
                         aaaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                       >  2:4449684/1‑50 (MQ=255)
                         aaaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                       >  2:4798730/1‑50 (MQ=255)
                          aaCGTGCATGGTATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                       >  2:4969102/1‑49 (MQ=38)
                          aaCGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTg                       <  2:294122/49‑1 (MQ=255)
                          aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                       <  1:5037994/49‑1 (MQ=255)
                          aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                       >  2:534987/1‑49 (MQ=255)
                           aCGTGCATGGCATGCTGGAGCCTCTCAGTATCAGGGGCACGAACGCTGCa                     <  2:1357217/50‑1 (MQ=255)
                           aCGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCa                     <  1:3392776/50‑1 (MQ=255)
                            cGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCaa                    <  1:2789574/50‑1 (MQ=255)
                             gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                   >  2:927724/1‑50 (MQ=255)
                             gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                   <  1:317555/50‑1 (MQ=255)
                             gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                   >  1:1436495/1‑50 (MQ=255)
                              tGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                   >  1:2900417/1‑49 (MQ=255)
                              tGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                   >  1:13634/1‑49 (MQ=255)
                              tGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATg                  <  2:1694614/50‑1 (MQ=255)
                               gcatggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGa                 >  1:1806193/1‑50 (MQ=255)
                               gcatggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGa                 >  1:2588172/1‑50 (MQ=255)
                                catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt                >  1:3831922/1‑50 (MQ=255)
                                catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt                <  2:3278104/50‑1 (MQ=255)
                                catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt                <  2:1430690/50‑1 (MQ=255)
                                catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt                <  2:2707844/50‑1 (MQ=255)
                                 atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt               >  1:4107703/1‑50 (MQ=255)
                                 atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt               >  2:5256811/1‑50 (MQ=255)
                                  tggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttt              <  1:2780547/50‑1 (MQ=255)
                                    gcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttttt            >  2:4333278/1‑50 (MQ=255)
                                     catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc           >  2:4640105/1‑50 (MQ=255)
                                     catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc           <  2:5160687/50‑1 (MQ=255)
                                     catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc           <  2:693270/50‑1 (MQ=255)
                                     catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc           <  2:2994274/50‑1 (MQ=255)
                                     catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc           <  2:255201/50‑1 (MQ=255)
                                      atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCt          <  2:2196508/50‑1 (MQ=255)
                                       tgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTa         >  2:2383396/1‑50 (MQ=255)
                                       tgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTa         <  1:1704299/50‑1 (MQ=255)
                                        gCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAt        >  2:3437400/1‑50 (MQ=255)
                                         cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAtt       <  1:3729354/50‑1 (MQ=255)
                                         cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAtt       >  1:1338198/1‑50 (MQ=255)
                                               |                                                
AATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGAT  >  NC_000913/118981‑119076

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