Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 115,634 Δ1 bp coding (415/552 nt) ampD → 1,6‑anhydro‑N‑acetylmuramyl‑L‑alanine amidase AmpD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273115,6320G.88.9% 90.0 / 4.7 27coding (413/552 nt)ampD1,6‑anhydro‑N‑acetylmuramyl‑L‑alanine amidase AmpD
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base . (12/12);  total (15/12)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAA  >  NZ_CP009273/115505‑115758
                                                                                                                               |                                                                                                                               
catgCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGACtgattga                                                                                                                     >  1:113336/1‑139 (MQ=255)
 atgCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACtgattgat                                                                                                                    <  1:311030/139‑1 (MQ=255)
    cGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGattgattg                                                                                                                  <  1:153643/138‑1 (MQ=255)
          tctctcAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGGTTGATTGCTCTCCg                                                                                                           >  2:211410/1‑139 (MQ=255)
          tctctcAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGGTAAGCGGCACTGGATGATTGCTATccc                                                                                                           >  1:178618/1‑138 (MQ=255)
             ctcAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGa                                                                                                         >  1:272645/1‑138 (MQ=255)
                 gTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGAATGATTTCTTTCCGGATATCg                                                                                                    >  1:392839/1‑139 (MQ=255)
                      aGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGGTTTGTTTCTATTCGGGTATCGCTaa                                                                                                >  2:189839/1‑138 (MQ=255)
                       ggggCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGATATCGCTaaaa                                                                                              <  1:138967/139‑1 (MQ=255)
                          gcgcGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGATATCGCTaaaaaaa                                                                                           >  1:287606/1‑137 (MQ=255)
                                                          ttGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCgg                                                           <  1:444881/139‑1 (MQ=255)
                                                                          cGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTc                                           <  2:272645/139‑1 (MQ=255)
                                                                          cGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTc                                           <  2:436539/139‑1 (MQ=255)
                                                                                           ccGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGc                           <  2:146933/138‑1 (MQ=255)
                                                                                             gATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGCCGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACg                       <  2:392839/139‑1 (MQ=255)
                                                                                             gATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACg                       <  2:345410/139‑1 (MQ=255)
                                                                                              aTGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTTCGCCGGATCGGGAAACCGATCCCGGTCCTGCATTTGATTGGGGACgg                      >  2:326763/1‑139 (MQ=255)
                                                                                              aTGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTTCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACgg                      >  1:351084/1‑139 (MQ=255)
                                                                                                 cgcAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGGTTGATTGCTATCCGGATATCGCTAAAAACATGACGGGGCATTGTGATATTTCGCCGGATCGGGAAAACGATCCCGGTCCTGCATTTGATTGGGGACGGttt                   >  2:426657/1‑139 (MQ=255)
                                                                                                 cgcAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGtt                    <  1:296202/138‑1 (MQ=255)
                                                                                                   cAGTATCAACAGCTTGCGGCGGTTAAGCGGGACTGATTGATTGCTATCCGGGTATCGCCAAAAAAATGACGGGCCCTTGTGATATTTCGCCGGAGCGGGAAAACGAACCCCGGCCCGCCTTTGATTTGGCACCGGtttg                  >  2:203649/1‑138 (MQ=255)
                                                                                                           aCAGCTTGCGGCGGTTACGCGGCACTGATTTATTTCTATCCGGGTATCGCTAAAAACATGACGGGCCATTGTGATATTTCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGGTTCGTGTGCTgg         >  1:392685/1‑139 (MQ=255)
                                                                                                                ttGCGGCGGTTACGCGGCCACGATTTATTGCTATCCGGATATCGCTAAAAAAATGACGGGCCCTTGTGATATTTCGCCGGATCCGGAAACCGATCCCGGTCCTTCCATTGATTGGGCACGGATTTGTGTGATGGTTAgg    >  1:275137/1‑138 (MQ=255)
                                                                                                                ttGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACg                       <  2:333603/120‑1 (MQ=255)
                                                                                                                ttGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAg     >  2:434636/1‑138 (MQ=255)
                                                                                                                ttGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGATATCGCCAAAAACATGACGGGCCATTGTGATATTTCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACg                       >  1:333603/1‑120 (MQ=255)
                                                                                                                   cggcggTTACGCGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCaa  <  2:275137/138‑1 (MQ=255)
                                                                                                                               |                                                                                                                               
CATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAA  >  NZ_CP009273/115505‑115758

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