Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,880,422 C→G 7.3% A738G (GCC→GGC)  mexF → multidrug RND transporter MexF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,880,4220CG 7.3% 173.3 / 4.6 55A738G (GCC→GGC) mexFmultidrug RND transporter MexF
Reads supporting (aligned to +/- strand):  ref base C (30/21);  new base G (2/2);  total (32/23)
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.31e-01

GTTCCGTCTGCAAATCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCAACGTCCAG  >  NC_002947/3880273‑3880555
                                                                                                                                                     |                                                                                                                                     
gTTCCGTCTGCAAATCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCAAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGGGCCACAAGTCGATgccgc                                                                                                                                       >  1:789659/1‑150 (MQ=255)
         gCAAATCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCg                                                                                                                              >  2:1491070/1‑150 (MQ=255)
          cAAATCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCgg                                                                                                                             <  2:1570036/150‑1 (MQ=255)
            aaTCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGa                                                                                                                           >  2:876874/1‑150 (MQ=255)
                gAATACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAg                                                                                                                       <  1:569004/150‑1 (MQ=255)
                       gTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGa                                                                                                                >  2:1044027/1‑150 (MQ=255)
                       gTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGa                                                                                                                >  1:1494779/1‑150 (MQ=255)
                        tGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGGCGGCATCGActtgt                                                                                                                            >  1:1410567/1‑133 (MQ=255)
                        tGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGGCGGCATCGActtgt                                                                                                                            <  2:1410567/137‑5 (MQ=255)
                            aaCCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCAc                                                                                                           >  1:1360483/1‑150 (MQ=255)
                                 gggTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGGGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGgtgt                                                                                                      <  2:104798/150‑1 (MQ=255)
                                   cttACTAAGCGCCCTACACGCAAACCACGACCATCACCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAACTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTgg                                                                                                    <  2:789659/149‑1 (MQ=255)
                                     tACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCg                                                                                                  >  1:1620033/1‑150 (MQ=255)
                                       cGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGAt                                                                                                >  1:1157806/1‑150 (MQ=255)
                                        gAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATc                                                                                               >  2:487461/1‑150 (MQ=255)
                                             gcTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGa                                                                                          <  1:73751/150‑1 (MQ=255)
                                                         aaCCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCa                                                                                                            <  1:1485757/120‑1 (MQ=255)
                                                         aaCCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCa                                                                                                            >  2:1485757/1‑120 (MQ=255)
                                                            ccAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCt                                                                           >  2:160329/1‑150 (MQ=255)
                                                              aGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGAAATCTTCGACACCCTGc                                                                         <  1:56047/150‑1 (MQ=255)
                                                                         gCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTg                                                              <  1:1491070/150‑1 (MQ=255)
                                                                              gagCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTc                                                         >  1:1141238/1‑150 (MQ=255)
                                                                                gCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGc                                                       <  1:1171233/150‑1 (MQ=255)
                                                                                  cacaACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGcc                                                                                                                    >  2:848860/1‑87 (MQ=255)
                                                                                  cacaACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCGAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGCGCTCGCTg                                                     >  2:1376314/1‑150 (MQ=255)
                                                                                  cacaACGCCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGcc                                                                                                                    <  1:848860/87‑1 (MQ=255)
                                                                                           ccAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCg                                                                                                                              <  2:696322/68‑1 (MQ=255)
                                                                                           ccAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCg                                                                                                                              >  1:696322/1‑68 (MQ=255)
                                                                                           ccAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAAc                                            >  1:505561/1‑150 (MQ=255)
                                                                                            cAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACg                                           >  1:1349249/1‑150 (MQ=255)
                                                                                              gagCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGAc                                         <  2:1349249/150‑1 (MQ=255)
                                                                                                  tGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGc                                                                                                                     >  1:849215/1‑70 (MQ=255)
                                                                                                  tGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGc                                                                                                                     <  2:849215/70‑1 (MQ=255)
                                                                                                   ggCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCaa                                    <  1:487461/150‑1 (MQ=255)
                                                                                                    gCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGAc                                                                                         >  1:1602539/1‑96 (MQ=255)
                                                                                                    gCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGAc                                                                                         <  2:1602539/96‑1 (MQ=255)
                                                                                                    gCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAAc                                   >  2:1624024/1‑150 (MQ=255)
                                                                                                          cTGTTCACCAGCTACCAGGTCAACGTGCCACGAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCttt                             >  1:776255/1‑150 (MQ=255)
                                                                                                           tGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTg                            >  1:202057/1‑150 (MQ=255)
                                                                                                            gTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTgg                           <  2:1360483/150‑1 (MQ=255)
                                                                                                               cACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCg                        >  1:1089073/1‑150 (MQ=255)
                                                                                                                aCCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACAGCTTTggagct                       >  1:14821/1‑146 (MQ=255)
                                                                                                                      tACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTAc                 >  1:126210/1‑150 (MQ=255)
                                                                                                                        ccAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCa               >  2:727527/1‑150 (MQ=255)
                                                                                                                         cAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCATGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTAc                                                                 <  1:1124924/99‑1 (MQ=255)
                                                                                                                         cAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCATGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTAc                                                                 >  2:1124924/1‑99 (MQ=255)
                                                                                                                          aGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAgg             >  1:413398/1‑150 (MQ=255)
                                                                                                                             tCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCttt                             >  2:1576773/1‑131 (MQ=255)
                                                                                                                             tCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCttt                             <  1:1576773/131‑1 (MQ=255)
                                                                                                                                   tGCCACAAGTCGATGCCGCCATCAACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCAACGTcc    <  1:160329/150‑1 (MQ=255)
                                                                                                                                     ccACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCAACGTCCAg  >  1:597801/1‑150 (MQ=255)
                                                                                                                                     ccACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCAACGTCCAg  <  1:1624024/150‑1 (MQ=255)
                                                                                                                                         cggTCGATGGCGGCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTggg                                                            <  2:1331304/86‑1 (MQ=255)
                                                                                                                                         cggTCGATGGCGGCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTggg                                                            >  1:1331304/3‑88 (MQ=255)
                                                                                                                                           gTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTc                                      <  1:729393/108‑1 (MQ=255)
                                                                                                                                           gTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTc                                      >  2:729393/1‑108 (MQ=255)
                                                                                                                                                   cgccATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAAc                                            <  1:981420/94‑1 (MQ=255)
                                                                                                                                                   cgccATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAAc                                            >  2:981420/1‑94 (MQ=255)
                                                                                                                                                     |                                                                                                                                     
GTTCCGTCTGCAAATCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCAACGTCCAG  >  NC_002947/3880273‑3880555

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: