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.92.9% 47.7 / ‑2.6 14coding (413/552 nt)ampD1,6‑anhydro‑N‑acetylmuramyl‑L‑alanine amidase AmpD
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (9/4);  total (10/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAG  >  NZ_CP009273/115505‑115755
                                                                                                                               |                                                                                                                           
catgCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACttattga                                                                                                                 >  2:192947/1‑139 (MQ=255)
                aGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGata                                                                                                   <  2:101627/137‑1 (MQ=255)
                            gcgAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGGGGGACTGTTTTATTTCTATCCGGATATCGCTaaaaaaaa                                                                                      >  2:361752/1‑135 (MQ=255)
                                                                       cACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGATTGATTGGTATCCGGAGATCGCTAAAAAAAAGACGGGGCAATGTGGTATTGCGCCGGAGCGGGAAACCGATTCCg                                          >  2:182580/1‑139 (MQ=255)
                                                                          cGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTc                                       <  2:1828/139‑1 (MQ=255)
                                                                                       tataCCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTTTGATATTGCGCCGGATCGGAAAAACGATCCCGGTCCTGCATTTGATTg                          >  2:126251/1‑139 (MQ=255)
                                                                                              aTGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGATTGATTTCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACgg                   >  1:370025/1‑139 (MQ=255)
                                                                                              aTGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGGTATCGCTAAAAACATGACGGGCCATTGTGATATTTCGCCGGATCGGAAAAACGATCCCGGTCCTGCATTTGATTGGGCACgg                   >  1:71869/1‑139 (MQ=255)
                                                                                              aTGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTc                                       <  1:62441/119‑1 (MQ=255)
                                                                                              aTGCGCAGTATCAACAGCTTGCGGCGGTTAAGCGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTc                                       >  2:62441/1‑119 (MQ=255)
                                                                                                                ttGCGGCGGTTACGCGGCACTGGTTGATTTCTATCCGGATATCGCTAAAAACATGACGGGCCATTTTGATATTTCGCCGGGTCGGGAAACCGAACCCGGGCCTGCCTTTGATTTGGCACGGGTTTGTGTGCTGGTCAg  >  2:129439/1‑138 (MQ=255)
                                                                                                                ttGCGGCGGTTACGCGGCACTGATTGATTTCTATCCGGATATCGCt                                                                                              >  2:454066/1‑46 (MQ=255)
                                                                                                                ttGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGATATCGCt                                                                                              <  1:454066/46‑1 (MQ=255)
                                                                                                                ttGCGGCGGTTACGCGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAg  >  1:227976/1‑138 (MQ=255)
                                                                                                                               |                                                                                                                           
CATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAG  >  NZ_CP009273/115505‑115755

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