Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 115,628 A→C T137P (ACG→CCG)  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,6280AC91.4% 82.3 / ‑5.6 35T137P (ACG→CCG) ampD1,6‑anhydro‑N‑acetylmuramyl‑L‑alanine amidase AmpD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (18/14);  minor base T (3/0);  total (21/14)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCC  >  NZ_CP009273/115542‑115714
                                                                                      |                                                                                      
gCAATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCcgc                                                                                     <  2:6529/90‑1 (MQ=255)
    tGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGc                                                                                 <  1:395961/90‑1 (MQ=255)
       tttttCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACt                                                                              <  2:309790/90‑1 (MQ=255)
          ttCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGGACtgat                                                                           >  2:75651/1‑90 (MQ=255)
                            tGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCCCTGATTGATTGCTATCCGGatat                                                         >  1:161734/1‑90 (MQ=255)
                            tGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCCCGGGGCCTGATTGATTTCTATCCGGatat                                                         >  1:719497/1‑90 (MQ=255)
                             gAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTTATTGCTATCCGGATATc                                                        >  1:28802/1‑90 (MQ=255)
                                 gCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTa                                                    <  2:407052/90‑1 (MQ=255)
                                  cACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTaa                                                   >  2:120931/1‑90 (MQ=255)
                                   aCCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTaaa                                                  >  2:163457/1‑90 (MQ=255)
                                          cGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAAGa                                           >  2:643426/1‑90 (MQ=255)
                                              ggCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACggg                                       >  2:82131/1‑90 (MQ=255)
                                                cGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGcc                                     <  2:651592/90‑1 (MQ=255)
                                                cGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGcc                                     <  1:229378/90‑1 (MQ=255)
                                                cGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGcc                                     <  2:183254/90‑1 (MQ=255)
                                                cGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGcc                                     <  2:200891/90‑1 (MQ=255)
                                                        gATGCGCAGTATCAACAGCTTGCGGCGGTTTCGCGGGCCCTGATTGATTTCTATCCGGATATCCCTAAAAAAAAGACGGGCCCTTGTGat                             >  1:631404/1‑90 (MQ=255)
                                                        gATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGGACTGATTGATTTCTATCCGGATATCGCTAAAAACATGACGGGCCCTTGTGat                             >  2:664382/1‑90 (MQ=255)
                                                        gATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGGACTGATTGATTGCTATCCGGATATCGCTAAAAAAAATACGGGCCATTGTGat                             >  1:618724/1‑90 (MQ=255)
                                                        gATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGat                             >  2:296495/1‑90 (MQ=255)
                                                        gATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGat                             >  2:651669/1‑90 (MQ=255)
                                                              cAGTATCAACAGCTTGCGGCGGTTTCGCGGGCACTGATTGATTGCTATCCCGATATTGCTAAAAAAAAGACGGGCCCTTTTGATAtttgg                       >  1:664677/1‑87 (MQ=255)
                                                               aGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGGTATCGCTAAAAACATGACGGGCCATTGTGATATTtcgc                      >  2:537827/1‑90 (MQ=255)
                                                                gTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGcc                     >  1:53785/1‑90 (MQ=255)
                                                                    cAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGAt                 <  1:105317/90‑1 (MQ=255)
                                                                        aGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGa             <  1:680338/90‑1 (MQ=255)
                                                                         gCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGaa            <  2:53785/90‑1 (MQ=255)
                                                                          cTTGCGGCGGTTTCGCGGGCACTGGTTGGGTGCTATCCGGGGATCGCTAAAAACAAGACGGGCCCTTTGGATATTGCGCCGGGTCGGaaa           >  1:600926/1‑90 (MQ=255)
                                                                          cTTGCGGCGGTTCCGCGGGCACTGATTGATTTCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGaaa           >  2:279166/1‑90 (MQ=255)
                                                                          cTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTTCGCCGGGTCGGaaa           >  2:573757/1‑90 (MQ=255)
                                                                          cTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGgaa           >  2:369078/1‑90 (MQ=255)
                                                                            tGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTTCGCCGGATCGGAAAAc         >  1:510559/1‑90 (MQ=255)
                                                                             gcggcgGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAAcc        <  1:82131/90‑1 (MQ=255)
                                                                              cggcggTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCg       <  2:212570/90‑1 (MQ=255)
                                                                                   gTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATccc  <  2:129074/90‑1 (MQ=255)
                                                                                      |                                                                                      
GCAATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCC  >  NZ_CP009273/115542‑115714

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