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,4880TC96.0% 150.9 / ‑1.3 50F90S (TTC→TCC) ampD1,6‑anhydro‑N‑acetylmuramyl‑L‑alanine amidase AmpD
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (14/34);  total (15/35)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTGAA  >  NZ_CP009273/115405‑115573
                                                                                   |                                                                                     
tttctttGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATaa                                                                                 <  1:1574466/90‑1 (MQ=255)
tttctttGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATaa                                                                                 <  2:1497451/90‑1 (MQ=255)
tttctttGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATaa                                                                                 <  2:193750/90‑1 (MQ=255)
tttctttGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATaa                                                                                 <  2:576185/90‑1 (MQ=255)
  tctttGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAAc                                                                               <  2:18002/90‑1 (MQ=255)
     ttGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTg                                                                            >  2:47944/1‑90 (MQ=255)
       gCTGAGATCGCCCATTTGCGCGTCTCCGCTCGCTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTgca                                                                          >  2:785556/1‑90 (MQ=255)
       gCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTgca                                                                          >  2:1055614/1‑90 (MQ=255)
        cTGAGATCGCCCATTTGCGCGTCTCCGCTCCCTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTgcat                                                                         <  1:690603/90‑1 (MQ=255)
        cTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTgcat                                                                         <  1:1215902/90‑1 (MQ=255)
        cTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTgcat                                                                         <  1:639155/90‑1 (MQ=255)
        cTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTgcat                                                                         <  2:366447/90‑1 (MQ=255)
        cTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTgcat                                                                         <  2:604223/90‑1 (MQ=255)
        cTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTgcat                                                                         <  1:1468795/90‑1 (MQ=255)
          gagaTCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTgcatgg                                                                       <  2:529538/90‑1 (MQ=255)
          gagaTCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTgcatgg                                                                       <  1:1236043/90‑1 (MQ=255)
          gagaTCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTgcatgg                                                                       >  1:1320343/1‑90 (MQ=255)
            gaTCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTgcatggca                                                                     <  1:899239/90‑1 (MQ=255)
            gaTCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTgcatgg                                                                       >  1:3265/1‑88 (MQ=255)
            gaTCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTgcatgg                                                                       <  2:3265/88‑1 (MQ=255)
                gCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCg                                                                 <  2:347955/90‑1 (MQ=255)
                  ccATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCggg                                                               <  1:1282760/90‑1 (MQ=255)
                    aTTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAg                                                             >  1:1425315/1‑90 (MQ=255)
                     tttGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGt                                                            <  2:1212018/90‑1 (MQ=255)
                      ttGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGACATGCGGGAGtc                                                           <  1:600684/90‑1 (MQ=255)
                        gcgcGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGtctc                                                         >  2:1142299/1‑90 (MQ=255)
                          gcgTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGtctctc                                                       <  1:730873/90‑1 (MQ=255)
                          gcgTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGtctctc                                                       <  1:1341395/90‑1 (MQ=255)
                           cgTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCa                                                      <  2:1158410/90‑1 (MQ=255)
                           cgTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCa                                                      <  2:754484/90‑1 (MQ=255)
                           cgTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCa                                                      <  2:1129977/90‑1 (MQ=255)
                             tctcCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGt                                                    <  2:1317850/90‑1 (MQ=255)
                                       cTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGgcg                                          >  1:1543963/1‑90 (MQ=255)
                                            tGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAAc                                     <  1:432781/90‑1 (MQ=255)
                                            tGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAAc                                     <  1:453275/90‑1 (MQ=255)
                                                   ccGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCaa                              >  2:956953/1‑90 (MQ=255)
                                                   ccGTGATGGTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCaa                              <  1:68944/90‑1 (MQ=255)
                                                          ggTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAtttt                       <  2:427304/90‑1 (MQ=255)
                                                           gTGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGAttttt                      <  1:943837/90‑1 (MQ=255)
                                                            tGAAATAGTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTc                     >  2:87600/1‑90 (MQ=255)
                                                                   gTCCAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTggg              >  1:579188/1‑90 (MQ=255)
                                                                     ccAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGAt            >  2:74709/1‑90 (MQ=255)
                                                                     ccAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGAt            <  2:565398/90‑1 (MQ=255)
                                                                      cAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGAtt           <  1:1055614/90‑1 (MQ=255)
                                                                      cAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGAtt           <  2:926244/90‑1 (MQ=255)
                                                                      cAGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGAtt           <  2:222788/90‑1 (MQ=255)
                                                                       aGTATGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTg          >  1:1275595/1‑90 (MQ=255)
                                                                          aTGTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGc       >  2:680632/1‑90 (MQ=255)
                                                                            gTTCCTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCtt     >  2:162850/1‑90 (MQ=255)
                                                                             ttccttccGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGc                                   >  1:474777/1‑59 (MQ=255)
                                                                             ttccttccGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGc                                   <  2:474777/59‑1 (MQ=255)
                                                                             ttccttccGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTg    >  2:557959/1‑90 (MQ=255)
                                                                             ttccttccGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTg    >  2:684384/1‑90 (MQ=255)
                                                                             ttccttccGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTg    >  1:5233/1‑90 (MQ=255)
                                                                             ttccttccGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTg    >  2:800236/1‑90 (MQ=255)
                                                                               ccTTCCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTGaa  <  2:1085509/90‑1 (MQ=255)
                                                                                   |                                                                                     
TTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCGAACGCTGCAATGATTTTTCTATTGGGATTGAGCTTGAA  >  NZ_CP009273/115405‑115573

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