Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,880,419 C→G 9.2% A737G (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,4190CG 9.2% 115.7 / 2.7 40A737G (GCC→GGC) mexFmultidrug RND transporter MexF
Reads supporting (aligned to +/- strand):  ref base C (17/19);  new base G (2/2);  total (19/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 = 9.15e-01

TCTGCAAATCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCAACGTCCAGGCCGAGCAGC  >  NC_002947/3880279‑3880565
                                                                                                                                            |                                                                                                                                                  
tCTGCAAATCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGa                                                                                                                                           >  1:1099886/1‑150 (MQ=255)
   gCAAATCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGGCGGCATCGActt                                                                                                                                        <  1:1098605/150‑3 (MQ=255)
    cAAATCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCgg                                                                                                                                       >  1:142183/1‑150 (MQ=255)
     aaaTCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGATCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCAt                                                                                                                                              >  1:799304/1‑142 (MQ=255)
     aaaTCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCAt                                                                                                                                              <  2:799304/142‑1 (MQ=255)
     aaaTCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCggg                                                                                                                                      <  1:961716/150‑1 (MQ=255)
      aaTCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTGTCGGAACGGTCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACACACCGAAGCCGCCATCGACCGGGa                                                                                                                                     >  1:448733/1‑150 (MQ=255)
      aaTCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGa                                                                                                                                     >  1:876167/1‑150 (MQ=255)
          gAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAg                                                                                                                                 <  2:1343850/150‑1 (MQ=255)
           aaGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAgg                                                                                                                                <  2:876167/150‑1 (MQ=255)
            aGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGc                                                                                                                               <  1:1058308/150‑1 (MQ=255)
                         cTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGGCGGCATCGActt                                                                                                                                        >  2:18356/1‑126 (MQ=255)
                         cTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGGCGGCATCGActt                                                                                                                                        <  1:18356/128‑3 (MQ=255)
                             gTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTgg                                                                                                              <  2:142183/150‑1 (MQ=255)
                                 cGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGAt                                                                                                          <  2:366333/150‑1 (MQ=255)
                                     gcgcTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCAcc                                                                                                      >  2:384717/1‑150 (MQ=255)
                                        cTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGAc                                                                                                   >  2:904079/1‑150 (MQ=255)
                                        cTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGAc                                                                                                   >  1:1067808/1‑150 (MQ=255)
                                           tACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATgcc                                                                                                                                                   <  2:405343/99‑1 (MQ=255)
                                           tACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATgcc                                                                                                                                                   >  1:405343/1‑99 (MQ=255)
                                              aaGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTc                                                                                             <  2:1099886/150‑1 (MQ=255)
                                               aGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCg                                                                                            >  2:41428/1‑150 (MQ=255)
                                                      ccAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCt                                                                                     >  1:767298/1‑150 (MQ=255)
                                                       cAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTg                                                                                    >  2:102160/1‑150 (MQ=255)
                                                                  cGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCGACAAGTCGATGCTGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCt                                                                         <  2:339453/150‑1 (MQ=255)
                                                                            cacaACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTg                                                               <  1:384717/150‑1 (MQ=255)
                                                                              caACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTa                                                             <  2:1067808/150‑1 (MQ=255)
                                                                                     ccAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAAc                                                      >  1:607640/1‑150 (MQ=255)
                                                                                       agagCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGa                                                    <  1:904079/150‑1 (MQ=255)
                                                                                                  gCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCTCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCt                                         >  2:45207/1‑150 (MQ=255)
                                                                                                  gCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGgtg                                                                                                                 >  1:456792/1‑78 (MQ=255)
                                                                                                  gCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGgtg                                                                                                                 <  2:456792/78‑1 (MQ=255)
                                                                                                    cTGTTCACCAGCTACCAGACCAACGTGCCACACATGATTGGCACGACCGACCGTGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGCTTACCTGCGCTCGCTGTACACCAACGACTTCAACCGCttt                                       >  2:15021/1‑150 (MQ=255)
                                                                                                    aggcctgccAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCttt                                       <  1:807433/143‑1 (MQ=255)
                                                                                                           ccAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCa                                >  2:51795/1‑150 (MQ=255)
                                                                                                            cAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCAc                               >  1:877854/1‑150 (MQ=255)
                                                                                                              gCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGATCCACGGTGTGGCGATCACCGACATCTTCTACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGGCTTCAACCGCTTTGGCCGCACtt                             <  2:448733/150‑1 (MQ=255)
                                                                                                                         aaCGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACAcc                                                                                       <  1:39935/81‑1 (MQ=255)
                                                                                                                         aaCGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACAcc                                                                                       >  2:39935/1‑81 (MQ=255)
                                                                                                                               ccACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCAACGTCCAg            <  2:607640/150‑1 (MQ=255)
                                                                                                                               ccACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCAACGTCCAg            <  1:102160/150‑1 (MQ=255)
                                                                                                                                         aTGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCAACGTCCAGGCCGagcagc  <  2:877854/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
TCTGCAAATCGAAGACCGTGGCAACCTGGGTTACGAAGCGCTGTACAAGGAAACCCAGAACATCATCGCCAAGAGCCACAACGTCCCAGAGCTGGCAGGCCTGTTCACCAGCTACCAGGTCAACGTGCCACAAGTCGATGCCGCCATCGACCGGGAAAAGGCCAAGACCCACGGTGTGGCGATCACCGACATCTTCGACACCCTGCAGGTTTACCTGGGCTCGCTGTACACCAACGACTTCAACCGCTTTGGCCGCACTTACCAGGTCAACGTCCAGGCCGAGCAGC  >  NC_002947/3880279‑3880565

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