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% 246.4 / NA 65V99A (GTC→GCC) ampD1,6‑anhydro‑N‑acetylmuramoyl‑L‑alanine amidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/29);  total (36/29)

AATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGG  >  NC_000913/118981‑119073
                                               |                                             
aaTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcct                                             >  1:5060284/1‑50 (MQ=255)
aaTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcct                                             <  1:36354/50‑1 (MQ=255)
aaTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcct                                             <  1:5435716/50‑1 (MQ=255)
 aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                            >  1:786733/1‑50 (MQ=255)
 aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                            >  1:4189952/1‑50 (MQ=255)
 aTAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctc                                            >  1:6053588/1‑50 (MQ=255)
  tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           <  1:4496644/50‑1 (MQ=255)
  tAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:793532/1‑50 (MQ=255)
    gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctctc                                          <  1:3344231/49‑1 (MQ=255)
    gTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGcctct                                           >  1:2032617/1‑48 (MQ=255)
     tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                        >  1:5322557/1‑50 (MQ=255)
     tCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                        <  1:2570244/50‑1 (MQ=255)
      ccAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAg                                        >  1:1946434/1‑49 (MQ=255)
       cAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGt                                       <  1:2577843/49‑1 (MQ=255)
        aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                     <  1:602207/50‑1 (MQ=255)
        aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                     <  1:5104417/50‑1 (MQ=255)
        aGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                     <  1:4585383/50‑1 (MQ=255)
         gTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                     >  1:4756174/1‑49 (MQ=255)
         gTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                     >  1:1589388/1‑49 (MQ=255)
         gTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATc                                    <  1:5330290/50‑1 (MQ=255)
          tATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTAt                                     >  1:1905325/1‑48 (MQ=255)
           aTGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAg                                  >  1:6322722/1‑50 (MQ=255)
            tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgg                                 >  1:3373188/1‑50 (MQ=255)
            tGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgg                                 <  1:3800123/50‑1 (MQ=255)
              ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgggg                               <  1:3281238/50‑1 (MQ=255)
              ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgggg                               >  1:6049759/1‑50 (MQ=255)
              ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAgggg                               >  1:3980818/1‑50 (MQ=255)
              ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                                <  1:5640800/49‑1 (MQ=255)
              ttCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAggg                                <  1:6070989/49‑1 (MQ=255)
               tCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                              >  1:5610866/1‑50 (MQ=255)
               tCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                              >  1:293109/1‑50 (MQ=255)
               tCCTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                              >  1:894916/1‑50 (MQ=255)
                ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcg                             >  1:4751339/1‑50 (MQ=255)
                ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                              <  1:478193/49‑1 (MQ=255)
                ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                              <  1:2275157/49‑1 (MQ=255)
                ccTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgc                              <  1:3292251/49‑1 (MQ=255)
                 cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                            >  1:3184321/1‑50 (MQ=255)
                 cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                            >  1:900868/1‑50 (MQ=255)
                 cTTTCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGgcgc                            <  1:1693370/50‑1 (MQ=255)
                  gttCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGcgcg                           <  1:5431125/49‑1 (MQ=255)
                    tCGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGaa                         >  1:1175339/1‑50 (MQ=255)
                     cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                        >  1:6006460/1‑50 (MQ=255)
                     cGATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                        >  1:625093/1‑50 (MQ=255)
                      gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAAc                        <  1:2435254/49‑1 (MQ=255)
                      gATAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACg                       <  1:5130470/50‑1 (MQ=255)
                       aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                      >  1:4461272/1‑50 (MQ=255)
                       aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                      <  1:6166465/50‑1 (MQ=255)
                       aTAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGc                      >  1:2039758/1‑50 (MQ=255)
                        tAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCt                     >  1:406853/1‑50 (MQ=255)
                        tAAACGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCt                     <  1:3632320/50‑1 (MQ=255)
                         aaaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                    >  1:91466/1‑50 (MQ=255)
                          aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                    >  1:765962/1‑49 (MQ=255)
                          aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                    >  1:1489021/1‑49 (MQ=255)
                          aaCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTg                    >  1:2480509/1‑49 (MQ=255)
                           aCGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCa                  <  1:5113033/50‑1 (MQ=255)
                            cGTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCaa                 >  1:4154972/1‑50 (MQ=255)
                             gTGCATGGCATGCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAAt                >  1:5827887/1‑50 (MQ=255)
                               gcatggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGa              >  1:4164692/1‑50 (MQ=255)
                               gcatggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGa              >  1:768998/1‑50 (MQ=255)
                                catggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt             >  1:4017641/1‑50 (MQ=255)
                                 atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt            <  1:5536350/50‑1 (MQ=255)
                                 atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt            <  1:500520/50‑1 (MQ=255)
                                 atggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtt            >  1:5298695/1‑50 (MQ=255)
                                  tggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttt           >  1:974049/1‑50 (MQ=255)
                                  tggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttt           <  1:2080648/50‑1 (MQ=255)
                                  tggcatgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttt           <  1:3953651/50‑1 (MQ=255)
                                     catgCGGGGGCCTCTCAGTAGCAGGGGCGCGAACGCTGCAATGATTTTTc        <  1:2266040/50‑1 (MQ=255)
                                     catgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc        >  1:2766077/1‑50 (MQ=255)
                                      atgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc        >  1:1553847/1‑49 (MQ=255)
                                       tgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTa      <  1:4168138/50‑1 (MQ=255)
                                       tgCGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTa      >  1:172144/1‑50 (MQ=255)
                                         cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAtt    <  1:3704945/50‑1 (MQ=255)
                                         cGGGAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAt     <  1:5769696/49‑1 (MQ=255)
                                          gggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTg   >  1:4918774/1‑50 (MQ=255)
                                           ggAGCCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTgg  <  1:2908658/50‑1 (MQ=255)
                                               |                                             
AATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGG  >  NC_000913/118981‑119073

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