Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,212,988 T→A 7.8% L19M (TTG→ATG)  mexC → multidrug RND transporter membrane fusion protein MexC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,212,9880TA 7.8% 137.7 / 3.1 51L19M (TTG→ATG) mexCmultidrug RND transporter membrane fusion protein MexC
Reads supporting (aligned to +/- strand):  ref base T (24/23);  new base A (2/2);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.69e-01

GTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGCCGAGGTGCG  >  NC_002947/3212848‑3213133
                                                                                                                                            |                                                                                                                                                 
gTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCAAGGCGATCGCCATGGccaat                                                                                                                                           <  1:9859/149‑4 (MQ=255)
 tCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGgt                                                                                                                                          <  6:146320/149‑1 (MQ=255)
   ttttATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGACCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGgtgt                                                                                                                                        >  8:84058/1‑149 (MQ=255)
           ttttAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCg                                                                                                                                 >  3:30563/1‑148 (MQ=255)
                    gaCAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGCGCCCGCagcagagc                                                                                                                        <  4:98803/148‑1 (MQ=255)
                     aCAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagagcag                                                                                                                      >  8:148948/1‑149 (MQ=255)
                       aaTATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGa                                                                                                                    >  6:30397/1‑149 (MQ=255)
                         tatTGACTGATGTGGGTTTGATAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCtt                                                                                                                                                  <  6:367433/117‑1 (MQ=255)
                         tatTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCAGTTTTCAATCAACGTTCGGGTGTGTGATGGGCAAGACAGGATCGATAAGGGTGGGCCCATTGGACATGGCGATCGCCtt                                                                                                                                                  >  5:367433/1‑117 (MQ=255)
                             gACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTg                                                                                                               >  1:234787/1‑148 (MQ=255)
                             gACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTg                                                                                                               <  7:148948/148‑1 (MQ=255)
                             gACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTg                                                                                                               <  7:29700/148‑1 (MQ=255)
                                tCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGc                                                                                                           >  7:206234/1‑149 (MQ=255)
                                tCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGc                                                                                                           >  8:119501/1‑149 (MQ=255)
                                  aTGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGcc                                                                                                          >  4:37103/1‑148 (MQ=255)
                                  aTGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGc                                                                                                           <  4:30563/147‑1 (MQ=255)
                                     tgGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAAt                                                                                                       <  3:372088/148‑1 (MQ=255)
                                             aaaaGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTAc                                                                                              >  2:206363/1‑149 (MQ=255)
                                               aaGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTAccc                                                                                            >  8:12571/1‑149 (MQ=255)
                                               aaGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGAGAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTAccc                                                                                            >  2:93185/1‑149 (MQ=255)
                                                aGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCg                                                                                           >  2:195479/1‑149 (MQ=255)
                                                  tAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGt                                                                                         <  6:230225/149‑1 (MQ=255)
                                                  tAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGt                                                                                         >  1:146433/1‑149 (MQ=255)
                                                    gCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTg                                                                                        <  7:84058/148‑1 (MQ=255)
                                                       tCCGCCCCGTTTTCAATCAAGCTTCGGCGGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGGGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGaaa                                                                                    <  2:79109/149‑1 (MQ=255)
                                                         cGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATc                                                                                  <  2:234787/149‑1 (MQ=255)
                                                         cGCCCAGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATc                                                                                  <  3:37103/149‑1 (MQ=255)
                                                           ccccGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGt                                                                                >  5:107949/1‑149 (MQ=255)
                                                           ccccGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGt                                                                                <  2:146433/149‑1 (MQ=255)
                                                                     aatcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGgc                                                                       <  5:30397/148‑1 (MQ=255)
                                                                      atcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGa                                                                              <  2:230390/140‑1 (MQ=255)
                                                                      atcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGa                                                                              >  1:230390/1‑140 (MQ=255)
                                                                                         aTGCGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTggc                                                  >  6:299599/1‑149 (MQ=255)
                                                                                            ggCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTggcggc                                               >  5:260951/1‑149 (MQ=255)
                                                                                                  aCAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACt                                         <  6:260951/149‑1 (MQ=255)
                                                                                                    aGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACACGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGTCTTGGCGGCTGAACTgc                                       >  3:28433/1‑149 (MQ=255)
                                                                                                       ttCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgccgg                                    >  5:318386/1‑149 (MQ=255)
                                                                                                           aTAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTTGCCTTGGCGGCTGAACTgccgggccg                                <  1:206363/149‑1 (MQ=255)
                                                                                                              aGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGt                             <  8:206234/149‑1 (MQ=255)
                                                                                                                ggTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTc                            >  4:157045/1‑148 (MQ=255)
                                                                                                                      acacccggCCAAGGCGATCGCCATGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCAGGGTCGCGCCGa                     >  3:41166/7‑149 (MQ=255)
                                                                                                                       gcATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagagc                                                                                                                        >  5:301133/1‑49 (MQ=255)
                                                                                                                       gcATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagagc                                                                                                                        <  6:301133/49‑1 (MQ=255)
                                                                                                                       gcATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGAt                    <  1:195479/149‑1 (MQ=255)
                                                                                                                         aTTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATgc                  <  7:119501/149‑1 (MQ=255)
                                                                                                                         aTTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATgc                  <  6:318386/149‑1 (MQ=255)
                                                                                                                          ccggCCAAGGCGATCGCCATGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgcc                                      <  4:215030/126‑1 (MQ=255)
                                                                                                                          ccggCCAAGGCGATCGCCATGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgcc                                      >  3:215030/3‑128 (MQ=255)
                                                                                                                                   gCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGCCGa        <  3:157045/149‑1 (MQ=255)
                                                                                                                                     gATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGCCGAgg      >  1:208153/1‑149 (MQ=255)
                                                                                                                                         gCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGc           >  5:190739/1‑140 (MQ=255)
                                                                                                                                         gCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGc           <  6:190739/140‑1 (MQ=255)
                                                                                                                                         gCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGCCGAGGTGCg  >  4:154509/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
GTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGCCGAGGTGCG  >  NC_002947/3212848‑3213133

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