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,6280AC87.1% 66.9 / 1.7 31T137P (ACG→CCG) ampD1,6‑anhydro‑N‑acetylmuramyl‑L‑alanine amidase AmpD
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base C (17/10);  total (20/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01

AATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATC  >  NZ_CP009273/115544‑115712
                                                                                    |                                                                                     
aaTGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCgg                                                                                  >  2:71180/1‑90 (MQ=255)
aaTGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCgg                                                                                  <  1:783650/90‑1 (MQ=255)
  tGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGc                                                                                <  2:836679/90‑1 (MQ=255)
  tGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGc                                                                                <  1:312266/90‑1 (MQ=255)
    aTTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCcgggcc                                                                              >  2:605164/1‑87 (MQ=255)
             ttGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGGACTgattgatt                                                                     >  2:822351/1‑90 (MQ=255)
             ttGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTgattgatt                                                                     <  1:279956/90‑1 (MQ=255)
               gggATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGc                                                                   <  1:523485/90‑1 (MQ=255)
                          tGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGattgattg                                                                    <  1:507940/78‑1 (MQ=255)
                          tGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGattgattg                                                                    >  2:507940/1‑78 (MQ=255)
                                cACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCAc                                                                              <  1:83672/62‑1 (MQ=255)
                                cACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCAc                                                                              >  2:83672/1‑62 (MQ=255)
                                          cTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCCCGGGCGCTGATTGAATGCTATCCGGATATCTCTAAAAAAATGACg                                        >  2:46351/1‑90 (MQ=255)
                                            ggCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGGTCCGCGGGCACTTATTGATTGCTATCCGGGTATCGCTAAAAACATGACggg                                      >  2:398522/1‑90 (MQ=255)
                                              cGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGcc                                    <  1:71180/90‑1 (MQ=255)
                                                tataCCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCAt                                  <  1:605164/90‑1 (MQ=255)
                                                  taCCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCAt                                  <  1:540121/88‑1 (MQ=255)
                                                  taCCGATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCAt                                  >  2:540121/1‑88 (MQ=255)
                                                      gATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGat                            >  1:129731/1‑90 (MQ=255)
                                                      gATGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAAAATGACGGGCCATTGTGat                            >  1:617526/1‑90 (MQ=255)
                                                      gATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAAGACGGGGCCTTGTGat                            >  2:164319/1‑90 (MQ=255)
                                                       aTGCGCAGTATCAACAGCTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGata                           >  1:296703/1‑90 (MQ=255)
                                                                        cTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGGCCTTGTGGTATTTCGCCGGGTTGGaaa         >  2:145520/1‑90 (MQ=255)
                                                                        cTTGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGaaa         >  1:64918/1‑90 (MQ=255)
                                                                        cTTGCGGCGGTTCCCCGGGGCCTGATTGATTGCTATCCGGGTATCGCTAAAAACATGAGGGGGCATTGTGATATTTCGCCGGGTCGGgaa         >  2:709579/1‑90 (MQ=255)
                                                                        cTTGCGGCGGTTACGCGGGGACTGATTGATTGCTATCCGGATATCGCTAAAAAAATGACGGGCCATTGTGATATTGCGCCGGATCGGaaa         >  1:628553/1‑90 (MQ=255)
                                                                        cTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                            <  2:137739/56‑1 (MQ=255)
                                                                        cTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                            >  1:137739/1‑56 (MQ=255)
                                                                          tGCGGCGGTTCCGCGGGCACTGATTGATTTCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGGTCGGAAAAc       >  2:533707/1‑90 (MQ=255)
                                                                          tGCGGCGGTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTTCGCCGGATCGGaaaa      >  1:825711/1‑89 (MQ=255)
                                                                                ggTTCCGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAAAATGACGGGCCATTGTGGTATTTCGCCGGATCGGGAAAAccaacc  >  2:607016/1‑86 (MQ=255)
                                                                                    |                                                                                     
AATGATTTTTCTATTGGGATTGAGCTTGAAGGCACCGATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATC  >  NZ_CP009273/115544‑115712

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