Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,880,630 A→C 6.9% D822A (GAC→GCC)  mexF → multidrug RND transporter MexF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,880,6300AC 6.9% 172.7 / 3.5 58D822A (GAC→GCC) mexFmultidrug RND transporter MexF
Reads supporting (aligned to +/- strand):  ref base A (24/30);  new base C (2/2);  total (26/32)
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.31e-01

GCCGCACTTACCAGGTCAACGTCCAGGCCGAGCAGCAGTTCCGCCTTGATGCCGAGCAGATCGGCCAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAGCTGCTGAAGGAAGAACTGCCCAACGGCATGACCTTCGAGTGGACCGACCT  >  NC_002947/3880486‑3880777
                                                                                                                                                |                                                                                                                                                   
gCCTCACTTACCAGGGAAACGTCCAGGCCGAGCAGCAGTTCCGCCTTGATGCCGAGCAGATCGGCCAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGTGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCgg                                                                                                                                                 <  5:13307/149‑1 (MQ=255)
gCCGCACTTACCAGGTCAACGTCCAGGCCGAGCAGCAGTTCCGCCTTGATGCCGAGCAGATCGGCCAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCgg                                                                                                                                                 <  6:220894/149‑1 (MQ=255)
   gCACTTACCAGGTCAACGTCCAGGCCGAGCAGCAGTTCCGCCTTGATGCCGAGCAGATCGGCCAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTc                                                                                                                                              <  6:114483/149‑1 (MQ=255)
    cACTTACCAGGTCAACGTCCAGGCCGAGCAGCAGTTCCGCCTTGATGCCGAGCAGATCGGCCAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCa                                                                                                                                             <  6:172605/149‑1 (MQ=255)
       ttACCAGGTCAACGTCCAGGCCGAGCAGCAGTTCCGCCTTGATGCCGAGCAGATCGGCCAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATg                                                                                                                                           >  3:16407/1‑148 (MQ=255)
          ccAGGTCAACGTCCAGGCCGAGCAGCAGTTCCGCCTTGATGCCGAGCAGATCGGCCAACTGAAAGCGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACt                                                                                                                                       <  1:313636/149‑1 (MQ=255)
               tCAACGTCCAGGCCGAGCAGCAGTTCCGCCTTGATGCCGAGCAGATCGGCCAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAAc                                                                                                                                  >  3:222880/1‑149 (MQ=255)
                 aaCGTCCAGGCCGAGCAGCAGTTCCGCCTTGATGCCGAGCAGATCGGCCAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACgg                                                                                                                                >  4:187094/1‑149 (MQ=255)
                       cAGGCCGTGGAGCAGTTCCGCCTTGATGCAGAGCTGATCGGCCATCTGAAAGTGCGCAACAAACTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTtgat                                                                                                                          <  7:77813/149‑1 (MQ=255)
                       cAGGCCGAGCAGCAGTTCCGCCTTGATGCCGAGCAGATCGGCCAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTtcat                                                                                                                          <  1:162067/149‑1 (MQ=255)
                       cAGGCCGAGCAGCAGTTCCGACTTGATGCCGAGCAGATCGGACAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGCGGCATGCACCACAACTGCTacat                                                                                                                          >  8:252417/1‑149 (MQ=255)
                            cGAGCAGCAGTTCCGCCTTGATGCCGAGCAGATCGGCCAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTg                                                                                                                     >  2:318544/1‑149 (MQ=255)
                                 agcagTTCCGCCTTGATGCCGAGCAGATCGGCCAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCgag                                                                                                                <  7:6846/149‑1 (MQ=255)
                                      ttCCGCCTTGATGCCGAGCAGATCGGCCAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCaa                                                                                                           >  1:241051/1‑149 (MQ=255)
                                              tGATGCCGAGCAGATCGGCCAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACaa                                                                                                                                   >  6:177180/1‑117 (MQ=255)
                                              tGATGCCGAGCAGATCGGCCAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACaa                                                                                                                                   <  5:177180/117‑1 (MQ=255)
                                                     gAGCAGATCGGCCAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGcccc                                                                                            <  4:214848/149‑1 (MQ=255)
                                                           aTCGGCCAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCggg                                                                                         <  8:39793/146‑1 (MQ=255)
                                                           aTCGGCCAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCggg                                                                                         >  7:39793/1‑146 (MQ=255)
                                                                 cAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCa                                                                                                                                         <  4:108051/92‑1 (MQ=255)
                                                                 cAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCa                                                                                                                                         >  3:108051/1‑92 (MQ=255)
                                                                 cAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTgc                                                                                                  <  6:165835/131‑1 (MQ=255)
                                                                 cAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTgc                                                                                                  >  5:165835/1‑131 (MQ=255)
                                                                   aCTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTg                                                                              <  3:248319/149‑1 (MQ=255)
                                                                      gAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGt                                                                            >  5:313487/1‑148 (MQ=255)
                                                                      gAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGt                                                                            >  4:22539/1‑148 (MQ=255)
                                                                                aacaacCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGc                                                                                        >  6:310770/1‑126 (MQ=255)
                                                                                aacaacCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGc                                                                                        <  5:310770/126‑1 (MQ=255)
                                                                                aacaacCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGc                                                                 >  8:146376/1‑149 (MQ=255)
                                                                                       ttGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATc                                                          >  5:235401/1‑149 (MQ=255)
                                                                                            gagaTGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGaa                                                     <  3:187094/149‑1 (MQ=255)
                                                                                            gagaTGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGaa                                                     >  8:202968/1‑149 (MQ=255)
                                                                                             agaTGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAg                                                    <  6:235401/149‑1 (MQ=255)
                                                                                                tgatCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAgctg                                                 <  2:241051/149‑1 (MQ=255)
                                                                                                ccgccccGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATGGAGAAgctg                                                 >  3:235524/5‑149 (MQ=255)
                                                                                                 gatCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATACACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGCAGCTGCGATCGAGAtgctgc                                                >  8:292624/1‑149 (MQ=255)
                                                                                                     ccGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAGCTGCTgaa                                            <  3:22539/149‑1 (MQ=255)
                                                                                                           gCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAGCTGCTgaaggaa                                        <  8:58613/147‑1 (MQ=255)
                                                                                                           gCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAGCTGCTgaaggaa                                        <  4:222880/147‑1 (MQ=255)
                                                                                                              aCCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAGCTGCTGAAGGAAGAACt                                   <  7:146376/149‑1 (MQ=255)
                                                                                                                   ccTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAGCTGCTGAAGGAAGAACTGccc                               >  4:281743/1‑148 (MQ=255)
                                                                                                                     tCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCga                                                        >  5:49850/1‑121 (MQ=255)
                                                                                                                     tCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCga                                                        <  6:49850/121‑1 (MQ=255)
                                                                                                                       aaGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAGCTGCTGAAGGAAGAACTGCCCAACg                           <  3:281743/148‑1 (MQ=255)
                                                                                                                           tCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAGCTGCTGAAGGAAGAACTGCCCAACGGCATg                      >  7:222338/1‑149 (MQ=255)
                                                                                                                              gCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAGCTGCTGAAGGAAGAACTGCCCAACGGCATGAcc                   <  2:183004/149‑1 (MQ=255)
                                                                                                                               cGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAGCTGCTGAAGGAAGAACTGCCCAACGGCATGACCt                  >  7:114150/1‑149 (MQ=255)
                                                                                                                                gACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAGCTGCTGAAGGAAGAACTGCCCAACGGCATGACCtt                 <  7:202968/149‑1 (MQ=255)
                                                                                                                                   agctgcGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAGCTGCTGAAGGAAGAACTGCCCAACGGCATGACCTTCGa              <  7:252417/144‑1 (MQ=255)
                                                                                                                                   acCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAGCTGCTGAAGGAAGAACTGCCCAACGGCATGACCTTc                >  5:68564/1‑147 (MQ=255)
                                                                                                                                     cTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATc                                                          >  4:7263/1‑103 (MQ=255)
                                                                                                                                     cTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATc                                                          <  3:7263/103‑1 (MQ=255)
                                                                                                                                      cccGGTCCGGCCCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAGCTGCTGAAGGAAGAACTGCCCAACGGCATGACCTTCGAGTg           >  1:215006/2‑149 (MQ=255)
                                                                                                                                      cccGGTCCGGCCCGGGTCATGCACTACAACGGCGTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAGCTGCTGAAGGAAGAACTGCCCAACGGCATGACCTTCGAGTg           >  1:219963/2‑149 (MQ=255)
                                                                                                                                        cGGGCCGGACCGGGTCATGCACTACAACGTCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAGCTGCTGAAGGAAGAACTGCCCAACGGCATGACCTTCGAGTGga         <  6:342549/149‑1 (MQ=255)
                                                                                                                                             cGGCCCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAGCTGCTGAAGGAAGAACTGCCCAACGGCATGACCTTCGAGTGgaccgac    <  2:219963/149‑1 (MQ=255)
                                                                                                                                             cGGCCCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAGCTGCTGAAGGAAGAACTGCCCAACGGCATGACCTTCGAGTGgaccgac    <  2:215006/149‑1 (MQ=255)
                                                                                                                                               gACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAGCTGCTGAAGGAAGAACTGCCCAACGGCATGACCTTCGAGTGGACCgacct  <  8:114150/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCCGCACTTACCAGGTCAACGTCCAGGCCGAGCAGCAGTTCCGCCTTGATGCCGAGCAGATCGGCCAACTGAAAGTGCGCAACAACCTTGGCGAGATGATCCCGCTGGCGACCTTCCTCAAGGTCAGCGACACCTCCGGGCCGGACCGGGTCATGCACTACAACGGCTTCATCACTGCCGAGATCAACGGCGCTGCCGCCCCGGGCTACAGCTCTGGTCAGGCGGAAGCTGCGATCGAGAAGCTGCTGAAGGAAGAACTGCCCAACGGCATGACCTTCGAGTGGACCGACCT  >  NC_002947/3880486‑3880777

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