Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,416,868 G→A 10.4% R37Q (CGG→CAG)  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,8680GA10.4% 127.6 / 5.9 48R37Q (CGG→CAG) PP_1239metallo‑beta‑lactamase family protein
Reads supporting (aligned to +/- strand):  ref base G (26/17);  new base A (3/2);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

CCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCT‑CAGCCAAGGGACTTTGCGCGGCATGC  >  NC_002947/1416719‑1417017
                                                                                                                                                       |                                                                                                                                                      
ccgccgCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAACGGAGCGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAgcg                                                                                                                                                        >  1:979857/1‑150 (MQ=255)
                 cTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACCCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTC‑‑CCTGGAACGCGCCGGGGCTCTTCACTCCGTGAAACCGAGCGGCGCCTGGCGCTGCTCg                                                                                                                                       >  2:578176/1‑150 (MQ=255)
                  tATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTc                                                                                                                                        >  1:255632/1‑148 (MQ=255)
                  tATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTc                                                                                                                                        <  2:255632/148‑1 (MQ=255)
                   aTTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCggg                                                                                                                                     <  1:637418/150‑1 (MQ=255)
                       aGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGttt                                                                                                                                 >  2:962989/1‑150 (MQ=255)
                         gACCACCTGGGCTAAGTGCGCTTCGCGGGACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGACCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCt                                                                                                                               >  1:309510/1‑150 (MQ=255)
                                     tAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCagctcagctc                                                                                                                   >  2:784623/1‑150 (MQ=255)
                                     tAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCagctcagctc                                                                                                                   >  2:557509/1‑150 (MQ=255)
                                         tGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTcagtg                                                                                                               <  2:428659/150‑1 (MQ=255)
                                              ttCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCAGCGCCAGGCGCCGCTCggt                                                                                                                                     >  1:418792/1‑122 (MQ=255)
                                              ttCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCAGCGCCAGGCGCCGCTCggt                                                                                                                                     <  2:418792/123‑2 (MQ=255)
                                                 gcgGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTg                                                                                                       >  2:1019853/1‑150 (MQ=255)
                                                     tACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCa                                                                                                   <  1:499782/150‑1 (MQ=255)
                                                          ggAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCt                                                                                                        <  1:1023464/140‑1 (MQ=255)
                                                          ggAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCt                                                                                                        >  2:1023464/1‑140 (MQ=255)
                                                          ggAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCAc                                                                                              >  1:221595/1‑150 (MQ=255)
                                                                ggCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTTATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGCTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACAt                                                                                        >  2:879845/1‑150 (MQ=255)
                                                                      cAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTgg                                                                                                                                                <  2:437272/88‑1 (MQ=255)
                                                                      cAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTgg                                                                                                                                                >  1:437272/1‑88 (MQ=255)
                                                                         gggAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCAc                                                                               <  1:962989/150‑1 (MQ=255)
                                                                          ggAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGAcccca                                                                              >  2:930771/1‑147 (MQ=255)
                                                                                    gCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGt                                                                    >  2:804758/1‑150 (MQ=255)
                                                                                        aTCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCCCATGCACGTGGTCGGc                                                                <  2:221595/150‑1 (MQ=255)
                                                                                           gCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTGTTTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTg                                                             >  2:1112965/1‑150 (MQ=255)
                                                                                             cAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCt                                                           <  1:557509/150‑1 (MQ=255)
                                                                                                 ggtgACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCa                                                       <  2:981849/150‑1 (MQ=255)
                                                                                                     acacGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCggg                                                                                                                                     >  1:329230/1‑68 (MQ=255)
                                                                                                     acacGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCggg                                                                                                                                     <  2:329230/68‑1 (MQ=255)
                                                                                                       acGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGCAACCGAGCGGCGCCTGGCGCTGCTCGGGGTATCTGCAGCTCCGATCAGTGCAGTGCTGGTCATCGACGAACATGTCGACTTCTTGCTTTTGGTCGGCTTGCTTTCACGgcgc                                                 >  1:715368/1‑150 (MQ=255)
                                                                                                       acGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGcag                                                                                                                           >  1:539268/1‑76 (MQ=255)
                                                                                                       acGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGcag                                                                                                                           <  2:539268/76‑1 (MQ=255)
                                                                                                       acGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGgcgc                                                 >  1:175723/1‑150 (MQ=255)
                                                                                                         gTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTa                                               >  1:1053344/1‑150 (MQ=255)
                                                                                                          ttCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACAGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTAc                                              <  1:755334/150‑1 (MQ=255)
                                                                                                                 tGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCAGCGCCAGGCGCCGCTCggtt                                                                                                                                    <  2:978008/57‑3 (MQ=18)
                                                                                                                 tGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCAGCGCCAGGCGCCGCTCggtt                                                                                                                                    >  1:978008/1‑55 (MQ=21)
                                                                                                                   g‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCg                                     <  1:784623/150‑1 (MQ=255)
                                                                                                                       g‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTc                                  >  1:616887/1‑150 (MQ=255)
                                                                                                                          ttGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCa                                                                                               <  1:75773/87‑1 (MQ=255)
                                                                                                                          ttGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCa                                                                                               >  2:75773/1‑87 (MQ=255)
                                                                                                                               gCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCT‑Ca                          >  2:670105/1‑150 (MQ=255)
                                                                                                                                     cccTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCT‑CAGCCAAg                    <  1:983687/150‑1 (MQ=255)
                                                                                                                                      ccTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCGGCATCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCT‑CAGCCAAgg                   >  2:770899/1‑150 (MQ=255)
                                                                                                                                          cGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTGGATGCATTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCT‑CAGCCAAGGGACt               <  2:309510/150‑1 (MQ=255)
                                                                                                                                             gAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCT‑CAGCCAAGGGACTTTg            >  1:646320/1‑150 (MQ=255)
                                                                                                                                               ccccGAGCAGCGCCAGGCGCCGCTCGGTGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCT‑CAGCCAAGGGACTTTgcg          >  2:137176/3‑150 (MQ=255)
                                                                                                                                                  cGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGAAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGACTTGCTGTCACGGCGCGACAATGTACCGGTCTACCT‑CAGCCAAGGGACTTTGCGCgg       >  2:519006/1‑150 (MQ=255)
                                                                                                                                                   gAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCT‑CAGCCAAGGGACTTTGCGCGGc      >  2:915708/1‑150 (MQ=255)
                                                                                                                                                      cggcgCCTGGCGCTGCTCGGTGTTTCTGCAGCTCAGCTCAGTACAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGT‑TGCTTGCAGCCAAGGGACTTTGCGCGGCATg   <  1:578176/150‑1 (MQ=255)
                                                                                                                                                       ggcgCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCT‑CAGCCAAGGGACTTTGCGCGGCATgc  <  2:616887/150‑1 (MQ=255)
                                                                                                                                                       ggcgCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCT‑CAGCCAAGGGACTTTGCGCGGCATgc  <  1:1019853/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
CCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGG‑TCG‑ATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCT‑CAGCCAAGGGACTTTGCGCGGCATGC  >  NC_002947/1416719‑1417017

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