Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,416,934 A→C 9.6% H61P (CAC→CCC)  PP_1239 → metallo‑beta‑lactamase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,416,9340AC 9.6% 144.7 / 6.5 52H61P (CAC→CCC) PP_1239metallo‑beta‑lactamase family protein
Reads supporting (aligned to +/- strand):  ref base A (25/22);  new base C (3/2);  total (28/24)
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.03e-01

GCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGCCTGGAAG  >  NC_002947/1416797‑1417077
                                                                                                                                         |                                                                                                                                               
gCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATgggg                                                                                                                                      >  1:311376/1‑149 (MQ=255)
  tGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTc                                                                                                                                    >  5:333727/1‑149 (MQ=255)
     tCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGc                                                                                                                                 <  5:7293/149‑1 (MQ=255)
      cGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCt                                                                                                                                >  6:198476/1‑149 (MQ=255)
       gCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCtt                                                                                                                               <  5:198476/149‑1 (MQ=255)
          aGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACAAGCCGACCCCATGCACGTGGTCGGCATgt                                                                                                                             >  3:204579/1‑147 (MQ=255)
          aGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACAAGCCGACCCCATGCACGTGGTCGGCATgt                                                                                                                             <  4:204579/148‑2 (MQ=255)
                          tCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAAt                                                                                                            >  2:262080/1‑149 (MQ=255)
                             tGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTa                                                                                                         >  7:315308/1‑149 (MQ=255)
                               gTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTAc                                                                                                        <  6:333727/148‑1 (MQ=255)
                                   aTTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTc                                                                                                   >  5:246781/1‑149 (MQ=255)
                                            tttCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGc                                                                                          >  3:120232/1‑149 (MQ=255)
                                               cccTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGAAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCaa                                                                                       >  2:297481/1‑149 (MQ=255)
                                                ccTGCGTGAAACCGAGCGTCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCaa                                                                                       <  3:67240/148‑1 (MQ=255)
                                                ccTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAg                                                                                      <  7:150353/149‑1 (MQ=255)
                                                ccTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAg                                                                                      >  6:32504/1‑149 (MQ=255)
                                                ccTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAg                                                                                      <  3:69706/149‑1 (MQ=255)
                                                    cGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGAc                                                                                  <  8:316039/149‑1 (MQ=255)
                                                    cGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGAc                                                                                  <  1:297481/149‑1 (MQ=255)
                                                     gTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGgc                                                                                                                    <  6:183585/114‑1 (MQ=255)
                                                       gAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACttt                                                                               >  1:45546/1‑149 (MQ=255)
                                                        aaaCCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTg                                                                              >  6:94360/1‑149 (MQ=255)
                                                          aCCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGgc                                                                                                                    >  4:233447/1‑109 (MQ=255)
                                                          aCCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGgc                                                                                                                    <  3:233447/109‑1 (MQ=255)
                                                           ccGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCg                                                                                                      <  3:321196/122‑1 (MQ=255)
                                                           ccGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCg                                                                                                      >  4:321196/1‑122 (MQ=255)
                                                              aGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGc                                                                        <  2:311376/149‑1 (MQ=255)
                                                                      cTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATgcgc                                                                  >  5:37996/1‑147 (MQ=255)
                                                                         gcgcTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCaa                                                                <  3:183512/146‑1 (MQ=255)
                                                                         gcgcTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCaa                                                                >  4:183512/1‑146 (MQ=255)
                                                                           gctgctCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCgg                                                           <  5:94360/149‑1 (MQ=255)
                                                                                 cGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAgg                                                     >  2:311496/1‑149 (MQ=255)
                                                                                  ggggTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGt                                                    >  6:330567/1‑149 (MQ=255)
                                                                                       ttCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCg                                               >  7:148183/1‑149 (MQ=255)
                                                                                            cagctcagctCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGtttt                                          >  4:133184/1‑149 (MQ=255)
                                                                                             agctcagctcAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTtc                                         <  5:236532/149‑1 (MQ=255)
                                                                                                       agtgcagtgcTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTgg                               <  5:32504/149‑1 (MQ=255)
                                                                                                         tgcagtgcTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGc                              >  7:158800/1‑148 (MQ=255)
                                                                                                              tgcTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTtc                                         >  2:20046/1‑132 (MQ=255)
                                                                                                                cTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCt                      <  5:330567/149‑1 (MQ=255)
                                                                                                                 tGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTg                     >  8:167501/1‑149 (MQ=255)
                                                                                                                   gTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGc                    <  6:246781/148‑1 (MQ=255)
                                                                                                                       cccACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATc               <  8:315308/149‑1 (MQ=255)
                                                                                                                         cACGTACATGCGGACCACGGTCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCgg             <  7:208498/149‑1 (MQ=255)
                                                                                                                         cACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCgg             >  5:223538/1‑149 (MQ=255)
                                                                                                                           cGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCa           <  6:37996/149‑1 (MQ=255)
                                                                                                                             agcaagCCGACCCCATGCACGTGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGc         >  5:11744/6‑149 (MQ=255)
                                                                                                                                 tGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGCCTgg     <  2:45546/149‑1 (MQ=255)
                                                                                                                                  gCCGACCCCATGCACGTGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCg                                               >  6:177078/1‑106 (MQ=255)
                                                                                                                                  gCCGACCCCATGCACGTGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCg                                               <  5:177078/106‑1 (MQ=255)
                                                                                                                                  gCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGCCTgg     >  1:325743/1‑148 (MQ=255)
                                                                                                                                     gACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGCCTGGAAg  >  7:34199/1‑148 (MQ=255)
                                                                                                                                         |                                                                                                                                               
GCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGTGGAGGTGGCCGGTTTTCTCGCGTGTGGCCAAAGCCTGCATATCGGCAGCCTGGAAG  >  NC_002947/1416797‑1417077

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