Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,221,967 C→A 7.1% L964M (CTG→ATG)  mexD → multidrug RND transporter MexD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,221,9670CA 7.1% 194.7 / 2.0 56L964M (CTG→ATG) mexDmultidrug RND transporter MexD
Reads supporting (aligned to +/- strand):  ref base C (36/15);  new base A (2/2);  total (38/18)
Fisher's exact test for biased strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-01

GACGTGTACTTCAAGGTCGGCCTGATCACCATCATTGGCCTGGCGGCGAAGAACGCCATCCTGATCGTCGAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCGTGATCGGCGGCATGCTCAGCGCC  >  NC_002947/3221832‑3222101
                                                                                                                                       |                                                                                                                                      
gacgTGTACTTCAAGGTCGGCCTGATCACCAACAATGGCCTGGCGGCGAAGAACGCCATCCTGATCGTCGAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGcc                                                                                                                                        >  2:78798/1‑136 (MQ=255)
               gTCGGCCTGATCACCATCATTGGCCTGGAGGCGAAGAACGCCATCCTGATCGTCGAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCTTCCGa                                                                                                                         >  1:438383/1‑136 (MQ=255)
                 cGGCCTGTTCTCCTTGATTGGCCTGGCGGCGAAGAACGCCATCCTGATCGTCGAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGACGCCATCGAGGCCGCGCGCGTGCGTTtccgtcagttc                                                                                                                       <  2:475238/136‑6 (MQ=37)
                   gCCTGATCACCATCATTGGCCTGGCGGCGAAGAACGCCATCCTGATCGTCGAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGt                                                                                                                     <  2:72458/136‑1 (MQ=255)
                   gCCTGATCACCATCATTGGCCTGGCGGCGAAGAACGCCATCCTGATCGTCGAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGt                                                                                                                     <  1:26701/136‑1 (MQ=255)
                      tGATCACCATCATTGGCCTGGCGGCGAAGAACGCCATCCTGATCGTCGAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGtgat                                                                                                                  >  2:511719/1‑136 (MQ=255)
                      tGATCACCATCATTGGCCTGGCGGCGAAGAACGCCATCCTGATCGTCGAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGtgat                                                                                                                  >  2:381125/1‑136 (MQ=255)
                      tGATCACCATCATTGGCCTGGCGGCGAAGAACGCCATCCTGATCGTCGAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCAGCGTTTCCGTCCGATCGtgat                                                                                                                  >  2:233751/1‑136 (MQ=255)
                      tGATCACCATCATTGGCCTGGCGGCGAAGAACGCCATCCTGATCGTCGAGTTCGCCAAGGAACTGTGGGAGAAGGCCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTACCGTCCGATCGtgaa                                                                                                                  >  2:404568/1‑135 (MQ=255)
                           aCCATCATTGGCCTGGCGGCGAAGAACGCCATCCTGATCGTCGAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCt                                                                                                             >  1:272466/1‑136 (MQ=255)
                                     gCCTGGCGGCGAAGAACGCCATCCTGATCGTCGAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCa                                                                                                          >  1:165225/1‑129 (MQ=255)
                                     gCCTGGCGGCGAAGAACGCCATCCTGATCGTCGAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCa                                                                                                          <  2:165225/129‑1 (MQ=255)
                                        tGGCGGCGAAGAACGCCATCCTGATCGTCGAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCAt                                                                                                <  2:545092/136‑1 (MQ=255)
                                        gggcggcGAAGAACGCCATCCTGATCGTCGAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCAt                                                                                                <  1:517380/135‑1 (MQ=255)
                                            ggcgAAGAACGCCATCCTGATCGTCGAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTc                                                                                            >  2:440014/1‑136 (MQ=255)
                                             gcgAAGAACGCCATCCTGATCGTCGAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCg                                                                                           >  1:142291/1‑136 (MQ=255)
                                             gcgAAGAACGCCATCCTGATCGTCGAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCg                                                                                           >  1:373177/1‑136 (MQ=255)
                                                           ccTGATCGTCGAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTg                                                                             <  2:511987/136‑1 (MQ=255)
                                                                  gtcgAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCg                                                                      >  2:41816/1‑136 (MQ=255)
                                                                  gtcgAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCg                                                                      >  1:611951/1‑136 (MQ=255)
                                                                  gtcgAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCg                                                                      >  1:607470/1‑136 (MQ=255)
                                                                  gtcgAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCg                                                                      >  1:292122/1‑136 (MQ=255)
                                                                              aaGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATg                                                                                                        <  1:162143/90‑1 (MQ=255)
                                                                              aaGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATg                                                                                                        >  2:162143/1‑90 (MQ=255)
                                                                                         ggAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCATGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGc                                               >  1:514009/1‑136 (MQ=255)
                                                                                         ggAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCATGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGc                                               >  1:517243/1‑136 (MQ=255)
                                                                                          gagaAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGcc                                              >  1:558288/1‑136 (MQ=255)
                                                                                               gggCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCATGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAgcgc                                         <  2:517243/136‑1 (MQ=255)
                                                                                               gggCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCATGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAgcgc                                         <  2:514009/136‑1 (MQ=255)
                                                                                                    aCAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCAt                                    >  1:335958/1‑136 (MQ=255)
                                                                                                    aCAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCAt                                    >  1:128588/1‑136 (MQ=255)
                                                                                                    aCAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCAt                                    >  1:152020/1‑136 (MQ=255)
                                                                                                             gCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCgg                           >  1:646008/1‑136 (MQ=255)
                                                                                                             gCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCgg                           >  1:20388/1‑136 (MQ=255)
                                                                                                             gCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCgg                           >  1:471770/1‑136 (MQ=255)
                                                                                                             gCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCgg                           >  1:290582/1‑136 (MQ=255)
                                                                                                             gCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCgg                           >  1:376486/1‑136 (MQ=255)
                                                                                                             gCGACGCCGCCATAGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTCGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCgg                           >  1:588616/1‑136 (MQ=255)
                                                                                                               gACGCCGCCATCGAGGCCGCGGGCCTGCGGTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCGGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCg                         >  2:140127/1‑136 (MQ=255)
                                                                                                               gACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCg                         >  1:106344/1‑136 (MQ=255)
                                                                                                               gACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCg                         >  2:444248/1‑136 (MQ=255)
                                                                                                                  gccgccATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCGTGa                      <  1:122658/136‑1 (MQ=255)
                                                                                                                  gccgccATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCGTGa                      <  1:89724/136‑1 (MQ=255)
                                                                                                                    cgccATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCGTGATc                    >  1:203864/1‑136 (MQ=255)
                                                                                                                    cgccATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCGTGATc                    >  2:583135/1‑136 (MQ=255)
                                                                                                                    cgccATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCGTGATc                    >  2:672658/1‑136 (MQ=255)
                                                                                                                    cgccATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTACGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCGTGAta                    >  2:401162/1‑135 (MQ=255)
                                                                                                                       cATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCGTGATcggc                 <  2:26136/136‑1 (MQ=255)
                                                                                                                        aTCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCGTGATcggcg                <  1:41816/136‑1 (MQ=255)
                                                                                                                         tCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCGTGa                      >  1:32806/1‑129 (MQ=255)
                                                                                                                         tCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCGTGa                      <  2:32806/129‑1 (MQ=255)
                                                                                                                         tCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCGTGATcggcgg               >  1:186655/1‑136 (MQ=255)
                                                                                                                            aGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCGTGATCGGCGGCAt            <  1:344620/136‑1 (MQ=255)
                                                                                                                                 gcgcgcCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCGTGATCGGCGGCATGCTCa       <  2:106344/136‑1 (MQ=255)
                                                                                                                                   gcgcCTGCATTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCGTGATCGGCGGCATGCTCAgc     >  1:570145/1‑136 (MQ=255)
                                                                                                                                     gcCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCGTGATCGGCGGCATGCTCAgcgc   >  1:566390/1‑136 (MQ=255)
                                                                                                                                     gcCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTGCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCGTGATCGGCGGCATGCTCAgcgc   <  2:163537/136‑1 (MQ=255)
                                                                                                                                      ccTGCGTTTCCGTCCGATCGTGATGACGTGCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCTTCGGTACCGGCGTGATCGGCGGCATGCTCAGCGcc  <  2:167613/136‑1 (MQ=255)
                                                                                                                                      ccTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCGTGATCGGCGGCATGCTCAGCGcc  <  2:20388/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
GACGTGTACTTCAAGGTCGGCCTGATCACCATCATTGGCCTGGCGGCGAAGAACGCCATCCTGATCGTCGAGTTCGCCAAGGAACTGTGGGAGAAGGGCTACAGCCTGTGCGACGCCGCCATCGAGGCCGCGCGCCTGCGTTTCCGTCCGATCGTGATGACCTCCATGGCGTTCATCCTCGGCGTGGTGCCGCTGGCCATCGCCTCCGGCGCGGGTGCGGCAAGCCAGCGCGCCATCGGTACCGGCGTGATCGGCGGCATGCTCAGCGCC  >  NC_002947/3221832‑3222101

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