Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,416,880 T→C 9.6% L41P (CTG→CCG)  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,8800TC 9.6% 140.4 / 5.2 52L41P (CTG→CCG) PP_1239metallo‑beta‑lactamase family protein
Reads supporting (aligned to +/- strand):  ref base T (27/20);  new base C (3/2);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

CTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGT  >  NC_002947/1416736‑1417027
                                                                                                                                                  |                                                                                                                                                    
cTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACCCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTC‑‑CCTGGAACGCGCCGGGGCTCTTCACTCCGTGAAACCGAGCGGCGCCTGGCGCTGCTCg                                                                                                                                                 >  2:578176/1‑150 (MQ=255)
 tATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTc                                                                                                                                                  <  2:255632/148‑1 (MQ=255)
 tATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTc                                                                                                                                                  >  1:255632/1‑148 (MQ=255)
  aTTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCggg                                                                                                                                               <  1:637418/150‑1 (MQ=255)
      aGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGttt                                                                                                                                           >  2:962989/1‑150 (MQ=255)
        gACCACCTGGGCTAAGTGCGCTTCGCGGGACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGACCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCt                                                                                                                                         >  1:309510/1‑150 (MQ=255)
                    tAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCagctcagctc                                                                                                                             >  2:784623/1‑150 (MQ=255)
                    tAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCagctcagctc                                                                                                                             >  2:557509/1‑150 (MQ=255)
                        tGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTcagtg                                                                                                                         <  2:428659/150‑1 (MQ=255)
                             ttCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCAGCGCCAGGCGCCGCTCggt                                                                                                                                               <  2:418792/123‑2 (MQ=255)
                             ttCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCAGCGCCAGGCGCCGCTCggt                                                                                                                                               >  1:418792/1‑122 (MQ=255)
                                gcgGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTg                                                                                                                 >  2:1019853/1‑150 (MQ=255)
                                    tACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCa                                                                                                             <  1:499782/150‑1 (MQ=255)
                                         ggAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCt                                                                                                                  <  1:1023464/140‑1 (MQ=255)
                                         ggAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCt                                                                                                                  >  2:1023464/1‑140 (MQ=255)
                                         ggAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCAc                                                                                                        >  1:221595/1‑150 (MQ=255)
                                               ggCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTTATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGCTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACAt                                                                                                  >  2:879845/1‑150 (MQ=255)
                                                        gggAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCAc                                                                                         <  1:962989/150‑1 (MQ=255)
                                                         ggAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGAcccca                                                                                        >  2:930771/1‑147 (MQ=255)
                                                                   gCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGt                                                                              >  2:804758/1‑150 (MQ=255)
                                                                       aTCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCCCATGCACGTGGTCGGc                                                                          <  2:221595/150‑1 (MQ=255)
                                                                          gCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTGTTTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTg                                                                       >  2:1112965/1‑150 (MQ=255)
                                                                            cAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCt                                                                     <  1:557509/150‑1 (MQ=255)
                                                                                ggtgACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCa                                                                 <  2:981849/150‑1 (MQ=255)
                                                                                    acacGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCggg                                                                                                                                               >  1:329230/1‑68 (MQ=255)
                                                                                    acacGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCggg                                                                                                                                               <  2:329230/68‑1 (MQ=255)
                                                                                      acGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGCAACCGAGCGGCGCCTGGCGCTGCTCGGGGTATCTGCAGCTCCGATCAGTGCAGTGCTGGTCATCGACGAACATGTCGACTTCTTGCTTTTGGTCGGCTTGCTTTCACGgcgc                                                           >  1:715368/1‑150 (MQ=255)
                                                                                      acGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGcag                                                                                                                                     >  1:539268/1‑76 (MQ=255)
                                                                                      acGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGcag                                                                                                                                     <  2:539268/76‑1 (MQ=255)
                                                                                      acGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGgcgc                                                           >  1:175723/1‑150 (MQ=255)
                                                                                        gTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTa                                                         >  1:1053344/1‑150 (MQ=255)
                                                                                         ttCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACAGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTAc                                                        <  1:755334/150‑1 (MQ=255)
                                                                                                tGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCAGCGCCAGGCGCCGCTCggtt                                                                                                                                              >  1:978008/1‑55 (MQ=21)
                                                                                                tGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCAGCGCCAGGCGCCGCTCggtt                                                                                                                                              <  2:978008/57‑3 (MQ=18)
                                                                                                  gTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCg                                               <  1:784623/150‑1 (MQ=255)
                                                                                                      gATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTc                                            >  1:616887/1‑150 (MQ=255)
                                                                                                         ttGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCa                                                                                                         <  1:75773/87‑1 (MQ=255)
                                                                                                         ttGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCa                                                                                                         >  2:75773/1‑87 (MQ=255)
                                                                                                              gCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCa                                    >  2:670105/1‑150 (MQ=255)
                                                                                                                    cccTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAg                              <  1:983687/150‑1 (MQ=255)
                                                                                                                     ccTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCGGCATCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAgg                             >  2:770899/1‑150 (MQ=255)
                                                                                                                         cGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTGGATGCATTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACt                         <  2:309510/150‑1 (MQ=255)
                                                                                                                            gAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTg                      >  1:646320/1‑150 (MQ=255)
                                                                                                                              ccccGAGCAGCGCCAGGCGCCGCTCGGTGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTgcg                    >  2:137176/3‑150 (MQ=255)
                                                                                                                                 cGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGAAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGACTTGCTGTCACGGCGCGACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCgg                 >  2:519006/1‑150 (MQ=255)
                                                                                                                                  gAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGc                >  2:915708/1‑150 (MQ=255)
                                                                                                                                     cggcgCCTGGCGCTGCTCGGTGTTTCTGCAGCTCAGCTCAGTACAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTTGCTTGCAGCCAAGGGACTTTGCGCGGCATg             <  1:578176/150‑1 (MQ=255)
                                                                                                                                      ggcgCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATgc            <  2:616887/150‑1 (MQ=255)
                                                                                                                                      ggcgCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATgc            <  1:1019853/150‑1 (MQ=255)
                                                                                                                                          ccTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCaa        >  2:566050/1‑150 (MQ=255)
                                                                                                                                           cTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGt                                                                                                               >  2:219607/1‑47 (MQ=255)
                                                                                                                                           cTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGt                                                                                                               <  1:219607/47‑1 (MQ=255)
                                                                                                                                                gctgctCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGt  <  2:175723/150‑1 (MQ=255)
                                                                                                                                                  tgctCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCAc                                                                >  1:417069/1‑87 (MQ=255)
                                                                                                                                                  tgctCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCAc                                                                <  2:417069/87‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATGCGCAAGCCGGT  >  NC_002947/1416736‑1417027

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: