Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,212,988 T→A 9.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 9.8% 104.0 / 3.5 41L19M (TTG→ATG) mexCmultidrug RND transporter membrane fusion protein MexC
Reads supporting (aligned to +/- strand):  ref base T (18/19);  new base A (2/2);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-01

AAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGCCGAG  >  NC_002947/3212844‑3213128
                                                                                                                                                |                                                                                                                                            
aaGAGTCTTTTATGATTTTAATGAGACAATATTGACTCAAGTGGGTTTGAAAAGTAGCATCCGCCACGTTTTCAATCAAGGTTCGTGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGAGCATTGGCAATGCCTATAGCCTTTGc                                                                                                                                          >  3:164217/1‑149 (MQ=255)
 agagTCTTTTATGATTTTAATGAGACAATATTGGCTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGcc                                                                                                                                         <  4:192691/149‑1 (MQ=255)
 agagTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGcc                                                                                                                                         >  1:31790/1‑149 (MQ=255)
 agagTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCAAGGCGATCGCCATGGc                                                                                                                                          <  1:12491/148‑1 (MQ=255)
  gagTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCg                                                                                                                                        >  2:84908/1‑149 (MQ=255)
     tCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGgt                                                                                                                                     <  2:31790/149‑1 (MQ=255)
          tATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGa                                                                                                                                 <  6:236101/148‑1 (MQ=255)
            tgtttttAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTATCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGc                                                                                                                               <  4:164217/145‑1 (MQ=255)
             gATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGccc                                                                                                                             >  1:61565/1‑149 (MQ=255)
               ttttAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGa                                                                                                                                 <  3:167529/143‑1 (MQ=255)
               ttttAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGa                                                                                                                                 >  4:167529/1‑143 (MQ=255)
                   aaTGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCgcag                                                                                                                         >  3:155800/1‑147 (MQ=255)
                    aggagaCAATATTGACGCATGTGGGTTTGAATAGTAGCATCCGCCCCGTTTTCACTCTAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGcagcag                                                                                                                      <  4:29319/147‑1 (MQ=255)
                        gaCAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagagc                                                                                                                   <  4:278320/148‑1 (MQ=255)
                         aCAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagagcag                                                                                                                 >  5:207202/1‑149 (MQ=255)
                          cAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAgg                                                                                                                <  1:84908/149‑1 (MQ=255)
                                           ggTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCg                                                                                               <  2:61565/149‑1 (MQ=255)
                                               tGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGc                                                                                            >  5:192892/1‑148 (MQ=255)
                                                        gCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTg                                                                                   >  2:104373/1‑148 (MQ=255)
                                                        gCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTg                                                                                   <  2:257337/148‑1 (MQ=255)
                                                        gCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTg                                                                                   <  6:207202/148‑1 (MQ=255)
                                                              gCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCTGGTGGAAATCg                                                                            >  1:108216/1‑149 (MQ=255)
                                                                         aatcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGgc                                                                  <  1:104373/148‑1 (MQ=255)
                                                                          atcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCGATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGgcgc                                                                <  2:108216/149‑1 (MQ=255)
                                                                          atcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGgcgc                                                                >  5:178275/1‑149 (MQ=255)
                                                                                         tgtgATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTg                                                                                   >  2:18466/1‑115 (MQ=255)
                                                                                         tgtgATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTg                                                                                   <  1:18466/115‑1 (MQ=255)
                                                                                             aTGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTggc                                             <  4:155800/149‑1 (MQ=255)
                                                                                               gggCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATc                                                                             >  4:117268/1‑115 (MQ=255)
                                                                                               gggCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATc                                                                             <  3:117268/115‑1 (MQ=255)
                                                                                                 gaaaGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGTCGATGGGCTTGGCCGGGTGTGAGCACGCAGCAGAGCATGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAg                                                           <  8:179237/129‑1 (MQ=255)
                                                                                                 gCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGaa                                                                                >  2:76182/1‑110 (MQ=255)
                                                                                                 gCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGaa                                                                                <  1:76182/110‑1 (MQ=255)
                                                                                                 gCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAg                                                           >  7:179237/1‑131 (MQ=255)
                                                                                                           ttCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgccgg                               >  5:204535/1‑149 (MQ=255)
                                                                                                                            acccggCCAAGGCGATCGCCATGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATg              >  6:316149/5‑149 (MQ=255)
                                                                                                                             aTTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATgc             >  8:288381/1‑149 (MQ=255)
                                                                                                                                    aaggCGATCGCCATGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTg                                        <  7:274564/113‑1 (MQ=255)
                                                                                                                                    aaggCGATCGCCATGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTg                                        >  8:274564/3‑115 (MQ=255)
                                                                                                                                           tCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGCCGAg  >  4:269982/1‑146 (MQ=255)
                                                                                                                                           tCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGCCGAg  <  3:269982/146‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                            
AAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGCCGAG  >  NC_002947/3212844‑3213128

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