Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,917,215 A→G 6.9% K333E (AAG→GAG)  mexB → multidrug resistance protein MexB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,917,2150AG 6.9% 170.5 / 2.4 58K333E (AAG→GAG) mexBmultidrug resistance protein MexB
Reads supporting (aligned to +/- strand):  ref base A (28/26);  new base G (2/2);  total (30/28)
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.60e-01

CCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATGGGCACCTTCGCCGTG  >  NC_002947/3917068‑3917361
                                                                                                                                                   |                                                                                                                                                  
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                                                                                                                      >  4:33021/1‑149 (MQ=255)
                             aCCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACAccc                                                                                                                      >  1:251166/1‑149 (MQ=255)
                              cccGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCt                                                                                                                     >  3:88687/1‑149 (MQ=255)
                               cgggtgacAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTg                                                                                                                    >  2:23516/8‑149 (MQ=255)
                                  gcgcCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTc                                                                                                                 >  5:33198/1‑149 (MQ=255)
                                  gcgcCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTc                                                                                                                 <  3:276189/149‑1 (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)
                                                                                                                                                gTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCAtgctgc                                                                                                    <  3:234549/52‑1 (MQ=255)
                                                                                                                                                gTCAAGGTGTCGATCGAGCAGGTCAACAAAACACTGTTCGAGGCCAtgctgc                                                                                                    >  4:234549/1‑52 (MQ=255)
                                                                                                                                                 tCAAGGTGTCGATGGATCAGGTCATCAACACGCTTTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTTCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCGAGTGTCATTGATGGGCACCTTCGCCGTg  <  8:240515/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATGGGCACCTTCGCCGTG  >  NC_002947/3917068‑3917361

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