Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,917,206 C→A 6.8% P330T (CCG→ACG)  mexB → multidrug resistance protein MexB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,917,2060CA 6.8% 190.6 / 3.7 59P330T (CCG→ACG) mexBmultidrug resistance protein MexB
Reads supporting (aligned to +/- strand):  ref base C (28/27);  new base A (2/2);  total (30/29)
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.15e-01

GCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATGG  >  NC_002947/3917063‑3917347
                                                                                                                                               |                                                                                                                                             
gCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTc                                                                                                                                          <  2:70570/149‑1 (MQ=255)
 cAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCg                                                                                                                                         <  8:175804/149‑1 (MQ=255)
   aaCCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGt                                                                                                                                        >  4:276189/1‑148 (MQ=255)
     ccGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCaa                                                                                                                                     <  8:283078/149‑1 (MQ=255)
      cGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAg                                                                                                                                    >  5:263983/1‑149 (MQ=255)
      cGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAg                                                                                                                                    >  3:88775/1‑149 (MQ=255)
             ggCCTTCAGCGTGCATTTGTCACCCGGCGCCAACGCCATGGAAACCTCTACCCTGGTACGGGCGAAGATGCAGGACCTTGCGCGCTACTTACCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGtcta                                                                                                                             <  3:205166/149‑3 (MQ=255)
                    aGCGTGCAGTTGTCACCCGGCGCCAACGCCATGTAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTTCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGc                                                                                                                       <  5:287480/148‑1 (MQ=255)
                    aGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCa                                                                                                                      <  7:140156/149‑1 (MQ=255)
                           aGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGTCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGtcatc                                                                                                               <  2:8857/149‑1 (MQ=255)
                               gTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAaca                                                                                                           >  5:318448/1‑149 (MQ=255)
                                  aCCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACAccc                                                                                                        >  1:251166/1‑149 (MQ=255)
                                  aCCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACAccc                                                                                                        >  4:33021/1‑149 (MQ=255)
                                   cccGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCt                                                                                                       >  3:88687/1‑149 (MQ=255)
                                    cgggtgacAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTg                                                                                                      >  2:23516/8‑149 (MQ=255)
                                       gcgcCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTc                                                                                                   <  3:276189/149‑1 (MQ=255)
                                       gcgcCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTc                                                                                                   >  5:33198/1‑149 (MQ=255)
                                          ccAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAg                                                                                                <  3:126531/149‑1 (MQ=255)
                                           cAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAgg                                                                                               <  6:263983/149‑1 (MQ=255)
                                                 cATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGtcatc                                                                                                               <  2:234285/127‑1 (MQ=255)
                                                 cATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGtcatc                                                                                                               >  1:234285/1‑127 (MQ=255)
                                                 cATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCAtgc                                                                                         >  7:134404/1‑149 (MQ=255)
                                                    ggAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCAtgctgc                                                                                      >  5:279433/1‑149 (MQ=255)
                                                        aCCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGgt                                                                                  >  6:45296/1‑149 (MQ=255)
                                                          cGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATACCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGgtgt                                                                                >  7:240515/1‑149 (MQ=255)
                                                           gCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTTGACGAACGGCGAGGTGtcgta                                                                                                                            <  8:259759/104‑3 (MQ=37)
                                                           gCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGtt                                                                               <  6:318448/149‑1 (MQ=255)
                                                               cccTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCAGCTGGTGTTCGCa                                                                           <  2:17364/149‑1 (MQ=255)
                                                                ccTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGACGCCATGCTGCTGGTGTTCGCAg                                                                          <  6:261664/149‑1 (MQ=255)
                                                                    gTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCa                                                                       <  8:134404/148‑1 (MQ=255)
                                                                    gTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCAt                                                                      <  2:75758/149‑1 (MQ=255)
                                                                        gggCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCAtgttc                                                                  <  6:33198/149‑1 (MQ=255)
                                                                                  gCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCAtgttc                                                                  >  6:190536/1‑139 (MQ=255)
                                                                                  gCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCAtgttc                                                                  <  5:190536/139‑1 (MQ=255)
                                                                                  gCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATACCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTTATCAACACCCTCTTCGAGGCCATGCTGCTGGTGTTAGCTGTCATGttcctgttcctg                                                         >  7:268372/1‑148 (MQ=255)
                                                                                   cAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCa                                                       <  4:88687/149‑1 (MQ=255)
                                                                                    aGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGttcctgttcctg                                                         >  1:203263/1‑146 (MQ=255)
                                                                                      gACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTGGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGaa                                                    >  3:298618/1‑149 (MQ=255)
                                                                                         cTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCt                                                 >  6:63793/1‑149 (MQ=255)
                                                                                          tGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTg                                                >  2:284763/1‑149 (MQ=255)
                                                                                            gcgcgcTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTgc                                               <  5:63793/148‑1 (MQ=255)
                                                                                                    cTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGc                                      <  3:33021/149‑1 (MQ=255)
                                                                                                       cgcGGAAGGCGTCAAGTACGACATCGCCTACGACACGTCGCCGTTCGTCAAGTTTTCGATCGATCTGG‑‑‑TCAACTCCCTGTTCGAGGCCATGCTGCTTTTGTTCGGAGTGATGTTCGTGTTCCTGCAGAACCTGCGCTTCACGCtta                                   <  2:50132/144‑3 (MQ=255)
                                                                                                         cGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATc                                 >  8:168103/1‑149 (MQ=255)
                                                                                                                       tACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGctgctg                                                                                    <  1:211377/84‑1 (MQ=255)
                                                                                                                       tACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGctgctg                                                                                    >  2:211377/1‑84 (MQ=255)
                                                                                                                       tACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCa                                                       <  4:150283/113‑1 (MQ=255)
                                                                                                                       tACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCa                                                       >  3:150283/1‑113 (MQ=255)
                                                                                                                              tCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTg            >  5:11034/1‑149 (MQ=255)
                                                                                                                               ccccTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTgg           >  1:125855/1‑149 (MQ=255)
                                                                                                                                 ccTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGACCCTGCGCTACACGCTGATCCCGACCCTGGTCGTGACAGTggga         >  8:237154/1‑147 (MQ=255)
                                                                                                                                     cGACACCTTGACGAACGGCGAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGa     >  6:122006/1‑149 (MQ=37)
                                                                                                                                      gACACCTTGACGAACGGCGAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGAt                                  <  2:3642/119‑1 (MQ=25)
                                                                                                                                      gACACCTTGACGAACGGCGAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGAt                                  >  1:3642/1‑119 (MQ=25)
                                                                                                                                        cacCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATgg  <  6:279433/149‑1 (MQ=255)
                                                                                                                                           cTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAAc                                                   <  3:290242/97‑1 (MQ=255)
                                                                                                                                           cTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAAc                                                   >  4:290242/1‑97 (MQ=255)
                                                                                                                                             cGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTgg                       >  2:212670/1‑123 (MQ=255)
                                                                                                                                             cGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTgg                       <  1:212670/123‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                             
GCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATGG  >  NC_002947/3917063‑3917347

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: ‑