Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,912,334 A→C 26.6% E824A (GAG→GCG)  mexB_2 → Multidrug resistance protein MexB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,912,3340AC26.6% 72.6 / 7.5 54E824A (GAG→GCG) mexB_2Multidrug resistance protein MexB
Reads supporting (aligned to +/- strand):  ref base A (11/27);  new base C (11/5);  total (22/32)
Fisher's exact test for biased strand distribution p-value = 1.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAACATTCACGTGCGCAACGACAGCGGCAAGATGGTGCCACTGGGCGCGTTCGTGCAGGCCAGGTGGGTCAGCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTCGATCTCAGGCGAGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGTCGAGCGGCTGGTGGCGCAGCTACCGGCCGGTACCGGTCTGGAATGGACCGGCCTGTCACTGCAGGAGCGGCTGTCCGGTAGCCAGGCCCC  >  ACB122/3912205‑3912466
                                                                                                                                 |                                                                                                                                    
tGAACATTCACGTGCGCAACGACAGCGGCAAGATGGTGCCACTGGGCGCGTGCGGGCAGGCCAGGTGGGTCAGCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTCGATCTCAGGCGAGCCTg                                                                                                                                 <  2:126529/135‑1 (MQ=255)
 gAACATTCACGTGCGCAACGACAGCGGCAAGATGGTGCCACTGGGCGCGTTCGTGCAGGCCAGGTGGGTCAGCGGCCCGGTGCAACTGCCCCGCTACAACGGTTACCCGGCGGTGTCGATCGCAGGCGCGCCTgc                                                                                                                                >  5:330850/1‑135 (MQ=255)
  aaCAATCACGTGCGCAACGACAGCGGCAAGATGGTGCCACTGGGCGCGTTCGTGCAGGCCTGGTGGGTCAGCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTCGATCTCAGGCGAGCCTgcc                                                                                                                               <  3:30105/135‑1 (MQ=255)
     aTTCACGTGCGCACCGCCAGCGGCAAGATGGTGCCACTGGGCGCGTTCGTGCAGGCCAGGTGGGGCAGCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTCGATCTCAGGCGAGCCTgccgcc                                                                                                                            <  5:254334/135‑1 (MQ=255)
     aTTCACGTGCGCAACGACAGCGGCAAGATGGTGCCACTGGGCGCGTTCGTGCAGGCCAGGTGGGTCAGCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTCGATCTCAGGCGAGCCTgccgcc                                                                                                                            <  5:123295/135‑1 (MQ=255)
                  aCGACAGCGGCAAGATGGTGCCACTGGGCGCGTTCGTGCAGGCCAGGTGGGTCAGCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTCGATCTCAGGCGAGCCTGCCGCCGGCTACAGCTCgg                                                                                                               <  4:333039/135‑1 (MQ=255)
                       aGCGGCAAGATGGTGCCACTGGGCGCGTTCGTGCAGGCCAGGTGGGTCAGCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTCGATCTCAGGCGAGCCTGCCGCCGGCTACAGCTCGGGCGAg                                                                                                          <  7:143286/135‑1 (MQ=255)
                            cAAGATGGTGCCACTGGGCGCGTTCGTGCAGGCCAGGTGGGTCAGCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTCGATCTCAGGCGAGCCTGCCGCCCGCTACAGCGCGGGCGAGGCCAt                                                                                                     >  4:226451/1‑135 (MQ=255)
                              aGATGGTGCCACTGGGCGCGTTCGTGCAGGCCAGGTGGGTCAGCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTCGATCTCACGCGAGCCTGCCGCCGGCGACCCCGCCCGCGAGGCCCTgg                                                                                                   >  4:68303/1‑135 (MQ=255)
                                     gCCACTGGGCGCGGTCGTGCAGGCCAGGTGGGTCAGCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTCGATCTCAGGCGCGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGt                                                                                            <  8:247592/135‑1 (MQ=255)
                                      ccACGGGGCGCGTTCGTGCAGGCCAGGTGGGTCAGCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTCGATCTCAGGCGAGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAggc                                                                                           <  4:180854/135‑3 (MQ=255)
                                       cACTGGGCGCGTTCGTGCAGGCCAGGTGGGTCAGCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTGGCTATCAGGCGAGCCAGCCGCCCGCTACAGCTCGGGCGCGGCCATGGCCGAAGt                                                                                            >  6:330923/1‑133 (MQ=255)
                                                 gTTCGTGCAGGCCAGGTGGGTCAGCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTCGCTCTCAGGCGAGCCAGCCGCCGGCGACAGCGCGGGCGAGGCCATGGCCGAAGTAGCGCGGCtggt                                                                                >  1:50220/1‑135 (MQ=255)
                                                  ttCGTGCAGGCCAGGGGGGTCAGCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTCGATCTCAGGCGAGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGTCGAGCGGCtggtg                                                                               <  8:410537/135‑1 (MQ=255)
                                                        cAGGCCCGGGGGGGCAGCGGCCCGGGGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTCGAGCTCAGGCGAGCCGGCCGCCGGCTACAGCTCGGGCGCGGCCAGGGCCGAAGGCGAGCGGCTGGTGGCGCAg                                                                         <  2:240599/135‑1 (MQ=255)
                                                        cAGGCCAGGGGGGTCAGCGGCCCGGTGCAACTGACCCGCTACACCGGTTACCCGGCGGTGTCGATCTCAGGCGAGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGGCGAGCGGCTGGTGGCGCAg                                                                         <  6:222047/135‑1 (MQ=255)
                                                         cggCCGGGGGGGGCCGCGGCCCGGTCCACCTCCCCCGCAACAACGGTTCCCCGGCGGGGTCGATCGCCGGCCAGCCGGCCGCCGGCTACAGCCCGGGCGAGGCCAGGCCCGAAGTCGAGCGGCTGGTGGCGCAGc                                                                        <  6:56655/134‑1 (MQ=255)
                                                         aGGCCAGGTGGGTCAGCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTCGCTCTCAGGCGCCCCGGCCCCCCGCGACCCCACCGGCGCGGCCCAGGCCGCACGCCGGCGGCCGGGGgcgccg                                                                         >  5:436679/1‑132 (MQ=255)
                                                         aGGCCAGGTGGGTCAGCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTCGCTCGCCGGCGCGCCAGCCGCCGGCTACCGCGCGGGCGCCGCCATGGCCGCACTCGAGCGGCTGGTGgcgccg                                                                         >  8:109162/1‑132 (MQ=255)
                                                           gCCAGGTGGGGCAGCGGCCCGGTGCAACTGACCCGCTACAACGGGTACCCGGCGGTGTCGATCTCAGGCGAGCCGGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGTCGAGCGGCTGGTGGCGCAGCTa                                                                      <  7:280274/135‑1 (MQ=255)
                                                              aGGTGGGTCAGCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGGCGATCTCAGGCGCGTCTGCCGCCCGCGACAGCGCGGGCGCGGCCCTGGGCGAAGTCGGGCGGCTGGGGGcgccgcggacg                                                                   >  1:419844/1‑130 (MQ=255)
                                                                  gggTCCGCGGCCCGGGGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTCGATCTCCGGCGAGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGTCGAGCGGCTGGTGGCGCAGCTAccggccg                                                               <  4:155461/135‑1 (MQ=255)
                                                                   gggcAGCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGGGTCGATCTCAGGCGAGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGTCGAGCGGCTGGTGGCGCAGCTAccggccgg                                                              <  4:157706/132‑1 (MQ=255)
                                                                   ggTCAGCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTCGAGCTCCGGCGAGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGTCGAGCGGCTGGTGGCGCAGCTAccggccgg                                                              <  4:216503/135‑1 (MQ=255)
                                                                       aGCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGGGTCGATCTCAGGCGCGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGTCGAGCGGCTGGTGGCGCAGCTACCGGCCGGTAcc                                                          <  4:312802/135‑1 (MQ=255)
                                                                        gCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTGGCTCTCAGGGGCGCCCGCCCCCGGCTACAGCTCGGGGGCGGCCATTGCCGGAGTCGCGCGGCGGGTGGCGCCCCTACCGGCCCGTACCg                                                         >  4:439199/1‑135 (MQ=255)
                                                                        gCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTCGATCTCCGGCGAGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGTCGAGCGGCTGGTGGCGCAGCTACCGGCCGGTTCCg                                                         >  6:223059/1‑135 (MQ=255)
                                                                        gCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTCGATCTCAGGCGCGCCTGCCCCCGGCTACAGCACGGGCGCGGCCATGGCCGCAGTCGAGCGGCAGGTGGCGCCGCTACCGGCCGGTACCg                                                         >  7:86853/1‑135 (MQ=255)
                                                                        gCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGGCGCTCTCCGGCGCGCCTGCCGCCGGCTACCGCGCGGGCGCGGCCCGGGCCGCAGTCGGGCGGCGGGTGGCGCCGCTACCGGCCGGTGCCg                                                         >  6:27477/1‑135 (MQ=255)
                                                                          ggCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTGGCTCTCAGGCGCGCCTGCCGCCGGCGACAGCTCGGGCGAGGCCCTGGGCGAAGTCGAGCGGCTGGTGGCGCAGCTACCGGCCGGTTCCggg                                                       >  3:35286/1‑134 (MQ=255)
                                                                                  gCAACTGACGCTCTACAAGGGTTACCGGTCGGGGGCGATCGCAGGCGAGCCGGCGGCCGGCTACAGCTCGGGGGAGGCCTTGGCGTAAGTCGAGCGGCTGGTGGCGCAGCTACCGGCCGGTACCGGTCTGGAATg                                               <  8:263227/135‑1 (MQ=255)
                                                                                  gCAACTCACCCGCTACAACGGTTACCCGGCGGTGTCGATCGCAGGCGAGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGTCGAGCGGCTGGTGGCGCAGCTACCGGCCGGTACCGGTCTGGAATg                                               <  6:176926/135‑1 (MQ=255)
                                                                                      cTGACCCGCTACAACGGTTACCCGGCGGTGTCGATCTCAGGCGAGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGTCGAGCGGCTGGTGGCGCAGCTACCGGCCGGTACCGGTCTGGAATGGAcc                                           <  8:343522/135‑1 (MQ=255)
                                                                                          cccGCTACAACGGTTACCCGGCGGTGGCGATCTCAGGCGAGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGTCGAGCGGCTGGTGGCGCAGCTACCGGCCGGTACCGGTCTGGAATGGACCGGcc                                       <  6:217505/135‑1 (MQ=255)
                                                                                              cTACAACGGTTACCCGGCGGTGTCGCTCTCAGGCGCGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGCAGTCGAGCGGCTGGTTGCGCAGCTACCGGCCGGTACCCGTCTGGAATGGACCGGCCTGtg                                   >  4:51343/1‑134 (MQ=255)
                                                                                              cTACAACGGTTACCCGGCGGTGTCGATCTCAGGCGAGCCTGCCGCCGGCTACAGCTCGGGCGCCGCCATGGCCGAAGTCGACCGGCTGGTGGCGCCCCTACCCGCCGGTACCCGGCTGGAATGGACCGGCCGGTc                                   >  4:253107/1‑135 (MQ=255)
                                                                                              cTACAACGGTTACCCGGCGGTGGCGCTCTCAGGCGAGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGTCGCGCGGCTGGTGGCGCAGCTACCGGCCGGTCCCGGTCTGGAATGGACCGGCCTGTc                                   >  8:40185/1‑135 (MQ=255)
                                                                                                   aCGGTTACCGGGCGGTGTCGATCTCAGGCGAGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGCCGGAAGGCGAGCGGCTGGTGGCGCAGCTACCGGCCGGTACCGGTCTGGAATGGACCGGCCTGTCACTGc                              <  3:402149/135‑1 (MQ=255)
                                                                                                        tccccGGGGGTGTGGATGGCAGGGGAGCGGGCCGCCGGCTACAGGCCGGGCGAGGCCATGGCGGAAGTCGAGCGGCTGGTTGCGCAGCTACCGGCCGGTACCGGTCTGGAATGGACCGGCCTGTCACTGCAGGAg                         <  4:129049/133‑1 (MQ=255)
                                                                                                        tACCCGGCGGTGTCGATCTCAGGCGAGCCTGCCGCCGGCTACAGCGCGGGCGAGGCCATGGCCGAAGTCGAGCGGCTGGTGGCGCAGCTACCGGCCGGTACCGGTCTGGAATGGACCGGCCTGTCACTGCAGGAg                         <  7:332350/135‑1 (MQ=255)
                                                                                                        tACCCGGCGGGGTCGATCGCAGGCGAGCCTGCCGCCGGCTACCGCTCGGGCGAGGCCATGGCCGAAGTCGAGCGGCTGGTGGCGCAGCTACCGGCCGGTACCGGTCTGGAATGGACCGGCCTGTCACTGCAGGAg                         <  8:185426/135‑1 (MQ=255)
                                                                                                        gtcccGGGGGTGGCGAGCTCAGGGGAGCCGGCCGCGGGCTACCGCTCGGGCGAGGCCTTGGCGGAAGTCGCGCGGCTGGTGGCGCAGCTACCGGCCGGTACCGGTCTGGAATGGACCGGCCTGTCACTGCAGGAg                         <  4:379459/133‑1 (MQ=255)
                                                                                                          cccGGCGGGGGCGATCCCCGGCGCGCCTGCCGCCGGCTACTGCTGGGGCGAGGCCGTGGCCGAAGGGGAGCGGCTGGTGGCGCAGCTACCGGCCGGTACCGGTCTGGAATGGACCGGCCTGTCACTGCAGGAGCg                       <  1:464712/135‑1 (MQ=255)
                                                                                                           ccGGCGGTGTCGATCTCAGGCGAGCCGGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGTCGAGCGGCTGGTGGCGCAGCTACCGGCCGGTACCGGTCTGGAATGGACCGGCCTGTCACTGCAGGAGCg                       <  7:375765/134‑1 (MQ=255)
                                                                                                            cggcggTGTCGATCTCAGGCGCGCCTGCCGCCGGCTACAGCTCGGGCGCGGCCATGGCCGAAGTCGCGCGGCTGGTGGCGCAGCTACCGGCCGGTACCGGTCTGGAATGGACCGGCCTGTCACTGCAGGAGCGGc                     >  3:228050/1‑135 (MQ=255)
                                                                                                            cggcggTGGCGATCTCAGGCGCGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGTCGAGCGGCTGGTGCCGCAGCTACCGGCCGGTACCGGTCTGGAATGGCCCGGCCGGTGACTGCAGGCGCGGc                     >  6:439504/1‑135 (MQ=255)
                                                                                                               cgttgtCGATATCAGGCTCGCCTGCCGCCGGATACAGCTCGGGCTAGGCCATGCCCTACGGAGAGCGGCTGGTTGCGCAGCTACAGGCCGGTACCGGTCTGTAATGGACCGGCCTGTCACTGCAGGAGCGGCTGt                  <  4:435236/132‑1 (MQ=255)
                                                                                                                ggTGTCGAGCGCCGGCGAGCCTGCCGCCGGCTACAGCGCGGGCGAGGCCAGGGCCGAAGGCGAGCGGCTGGTGGCGCAGCTACCGGCCGGTACCGGTCTGGAATGGACCGGCCTGTCACTGCAGGAGCGGCTGTc                 <  6:403123/135‑1 (MQ=255)
                                                                                                                      gATCTCAGGCGCGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCTTGGCCGAAGTCGAGCGGCTGGTGGCGCAGCTACCGGCCGGTACCGGTCTGGAATGGACCGGCCTGTCACTGCAGGAGCGGCTGTCCGGTAg           <  5:151464/135‑1 (MQ=255)
                                                                                                                          tcAGGCGAGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGTCGAGCGGCTGGTGGCGCAGCTACCGGCCGGTACCGGTCTGGAATGGACCGGCCTGTCACTGCAGGAGCGGCTGTCCGGTAGCCAg       >  5:447455/1‑135 (MQ=255)
                                                                                                                            aGGCGAGCCGGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGGCGAGCGGCTGGTGGCGCAGCTACCGGCCGGTACCGGTCTGGAATGGACCGGCCTGTCACTGCAGGAGCGGCTGTCCGGTAGCCAGGc     <  7:161154/135‑1 (MQ=255)
                                                                                                                               cGAGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGTCGCGCGGCTGGTGGCGCAGCTACCGGCCGGTTCCGGTCTGGAATGGACCCGCCTGTCACAGCACGAGCGGGTGTCCGGTAGCCCGGcccc  >  5:305925/1‑135 (MQ=255)
                                                                                                                               cGAGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGTCGAGCTGCTGGTGGCGCAGCTACCGGCAGGTACCGGGCTGGAATGGACCGGCCTGTCACTGCAGGAGCGGATCTCCGGTAGCCAGTcccc  >  2:325967/1‑135 (MQ=255)
                                                                                                                               cGAGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGTCGAGCGGCTGGTGGCGCAGCTACCGGCCGGTACCGGTCTGGAATGGACCGGCCTGTCACAGCAGGAGCGGCTGTCCGGTAGGCAGGcccc  >  5:387033/1‑135 (MQ=255)
                                                                                                                                 |                                                                                                                                    
TGAACATTCACGTGCGCAACGACAGCGGCAAGATGGTGCCACTGGGCGCGTTCGTGCAGGCCAGGTGGGTCAGCGGCCCGGTGCAACTGACCCGCTACAACGGTTACCCGGCGGTGTCGATCTCAGGCGAGCCTGCCGCCGGCTACAGCTCGGGCGAGGCCATGGCCGAAGTCGAGCGGCTGGTGGCGCAGCTACCGGCCGGTACCGGTCTGGAATGGACCGGCCTGTCACTGCAGGAGCGGCTGTCCGGTAGCCAGGCCCC  >  ACB122/3912205‑3912466

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