Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 115,488 T→C F90S (TTC→TCC)  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,4880TC90.6% 82.1 / 1.8 32F90S (TTC→TCC) ampD1,6‑anhydro‑N‑acetylmuramyl‑L‑alanine amidase AmpD
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base C (9/20);  total (9/23)
Fisher's exact test for biased strand distribution p-value = 5.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

CTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTG  >  NZ_CP009273/115413‑115571
                                                                           |                                                                                   
cTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTgcat                                                                       >  1:464287/1‑90 (MQ=255)
  gagaTCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTgcatgg                                                                     <  2:27630/90‑1 (MQ=255)
   agaTCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTgcatggc                                                                    >  1:752281/1‑90 (MQ=255)
            aTTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAg                                                           <  1:426204/90‑1 (MQ=255)
             tttGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGt                                                          <  2:97979/90‑1 (MQ=255)
             tttGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGt                                                          <  1:359523/90‑1 (MQ=255)
              ttGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGtc                                                         <  2:685370/90‑1 (MQ=255)
              ttGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGtc                                                         <  1:15230/90‑1 (MQ=255)
              ttGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGtc                                                         <  1:563095/90‑1 (MQ=255)
              ttGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGtc                                                         <  1:629090/90‑1 (MQ=255)
                gcgcGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGtctc                                                       <  2:731187/90‑1 (MQ=255)
                 cgcgTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGtctct                                                      <  1:775349/90‑1 (MQ=255)
                 cgcgTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGtctct                                                      <  2:464287/90‑1 (MQ=255)
                 cgcgTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGtctct                                                      <  2:41480/90‑1 (MQ=255)
                  gcgTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGtctctc                                                     <  1:384445/90‑1 (MQ=255)
                   cgTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCa                                                    <  1:331992/90‑1 (MQ=255)
                          gCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAg                                             <  2:231606/90‑1 (MQ=255)
                          gCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAg                                             <  1:308703/90‑1 (MQ=255)
                                    tGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAAc                                   >  2:103863/1‑90 (MQ=255)
                                      aTTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGc                                 <  2:95251/90‑1 (MQ=255)
                                         cGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGc                              >  1:702440/1‑90 (MQ=255)
                                         cGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGc                              >  1:322894/1‑90 (MQ=255)
                                           ccGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGtctctc                                                     >  1:291930/1‑65 (MQ=255)
                                           ccGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGtctctc                                                     <  2:291930/65‑1 (MQ=255)
                                           ccGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGtctct                                                      >  1:221689/1‑64 (MQ=255)
                                           ccGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGtctct                                                      <  2:221689/64‑1 (MQ=255)
                                               gATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAt                        >  2:432285/1‑90 (MQ=255)
                                                        aTAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTAtt               <  2:72573/90‑1 (MQ=255)
                                                           gTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTGTATTggg            >  2:384629/1‑90 (MQ=255)
                                                              cAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGAtt         <  1:607580/90‑1 (MQ=255)
                                                                  aTGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGc     <  2:440915/90‑1 (MQ=255)
                                                                  aTGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGGTGCAATGATTTTTCTATTGGGATTGAGc     <  2:792161/90‑1 (MQ=255)
                                                                     ttccttccGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTg  >  1:44616/1‑90 (MQ=255)
                                                                           |                                                                                   
CTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTG  >  NZ_CP009273/115413‑115571

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