Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,915,797 C→G 8.0% P256A (CCG→GCG)  PP_3455 → multidrug RND transporter membrane fusion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,915,7970CG 8.0% 155.2 / 4.6 50P256A (CCG→GCG) PP_3455multidrug RND transporter membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base C (21/25);  new base G (2/2);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

AGCAGTTGGACCCGATCCATGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGCGGGTACGCCTGGAGCAAGCGGT  >  NC_002947/3915651‑3915936
                                                                                                                                                  |                                                                                                                                           
agcagTTGGACCCGATCCATGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCg                                                                                                                                           <  7:19714/149‑1 (MQ=255)
     ttGGACCCGATCCATGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGcc                                                                                                                                      >  7:334662/1‑149 (MQ=255)
       ggACCCGATCCATGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGGGTAGGCGCTGCCg                                                                                                                                     <  4:138706/148‑1 (MQ=255)
        gACCCGATCCATGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCgg                                                                                                                                    <  2:206121/148‑1 (MQ=255)
         aCCCGATCCATGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCggg                                                                                                                                   <  8:334662/148‑1 (MQ=255)
         aCCCGATCCATGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGc                                                                                                                                  <  8:24849/149‑1 (MQ=255)
           ccGATCCATGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGGGTAGGCGCTGCCgtcag                                                                                                                                 >  3:99146/1‑144 (MQ=255)
           ccGATCCATGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGGGTAGGCGCTGCCgt                                                                                                                                    <  4:99146/145‑2 (MQ=255)
            cGATCCATGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAg                                                                                                                               <  3:239179/149‑1 (MQ=255)
            cGATCCATGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAg                                                                                                                               >  6:243838/1‑149 (MQ=255)
                  aTGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGctgctg                                                                                                                         >  6:150281/1‑149 (MQ=255)
                  aTGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGctgctg                                                                                                                         <  4:319985/149‑1 (MQ=255)
                     ccGATGTCACCCAGTCTACCCGCGAGCTCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGtt                                                                                                                       <  1:93119/148‑1 (MQ=255)
                                 aGTCTACCCGCGAGCTCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCa                                                                                                            >  6:275649/1‑147 (MQ=255)
                                 aGTCTACCCGCGAGCTCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCa                                                                                                            <  5:275649/147‑1 (MQ=255)
                                 aGTCTACCCGCGAGCTCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGc                                                                                                          <  5:150281/149‑1 (MQ=255)
                                       cccGCGAGCTCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGAc                                                                                                    <  5:243838/149‑1 (MQ=255)
                                        ccGCGAGCTCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGAcc                                                                                                   >  7:148094/1‑149 (MQ=255)
                                                tCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCacc                                                                                           >  6:210103/1‑149 (MQ=255)
                                                tCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCa                                                                                             >  8:64697/1‑147 (MQ=255)
                                                tCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCa                                                                                             <  7:64697/147‑1 (MQ=255)
                                                     gCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCaccaac                                                                                        >  7:149325/1‑147 (MQ=255)
                                                            gccgcgccCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGTCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCAcc                                                                               <  8:149325/149‑1 (MQ=255)
                                                              cgcgcaCTGCGCGCCGGAGAATTGAAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCTGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCt                                                                             >  2:285877/1‑149 (MQ=255)
                                                                       cgcgcCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGTGa                                                                    <  8:148094/149‑1 (MQ=255)
                                                                                  aTTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAAcc                                                         >  7:45095/1‑149 (MQ=255)
                                                                                             tGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGAc                                              <  5:210103/149‑1 (MQ=255)
                                                                                               ggCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACct                                            >  7:191614/1‑149 (MQ=255)
                                                                                                 cGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCAcc                                                                               >  4:305826/1‑112 (MQ=255)
                                                                                                 cGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCAcc                                                                               <  3:305826/112‑1 (MQ=255)
                                                                                                  gACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACctgct                                         >  6:326409/1‑149 (MQ=255)
                                                                                                        cAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAgg                                   >  3:289303/1‑149 (MQ=255)
                                                                                                        cAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGGTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCCCCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAgg                                   >  6:361725/1‑149 (MQ=255)
                                                                                                              cgcgCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCTGGCAGcag                             <  8:45095/149‑3 (MQ=255)
                                                                                                                   cACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGa                                               <  1:205726/126‑1 (MQ=255)
                                                                                                                   cACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGa                                               >  2:205726/1‑126 (MQ=255)
                                                                                                                   cACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGa                                               <  3:27763/126‑1 (MQ=255)
                                                                                                                   cACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGa                                               >  4:27763/1‑126 (MQ=255)
                                                                                                                     cGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGCgg                      <  7:302629/149‑1 (MQ=255)
                                                                                                                        tGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGc                        <  3:121094/144‑1 (MQ=255)
                                                                                                                        tGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGc                        >  4:121094/1‑144 (MQ=255)
                                                                                                                         gATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGCGGGTAc                  >  1:64802/1‑149 (MQ=255)
                                                                                                                             cAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGCGGGTACGCCt              <  2:64802/149‑1 (MQ=255)
                                                                                                                              aGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGCGGGTACGCCTg             <  2:193582/149‑1 (MQ=255)
                                                                                                                                gATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGCGGGTACGCCTGGa           <  4:289303/149‑1 (MQ=255)
                                                                                                                                    ccGGCAGCGGGTAGGCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGCGGGTACGCCTGGAGc         >  3:346893/2‑147 (MQ=255)
                                                                                                                                     cGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGCGGGTACGCCTGGAGCAAg      >  4:69951/1‑149 (MQ=255)
                                                                                                                                      ggCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCa                                                                        <  3:160995/82‑1 (MQ=255)
                                                                                                                                      ggCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCa                                                                        >  4:160995/1‑82 (MQ=255)
                                                                                                                                         aGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGCGGGTACGCCTGGAGCAAGCGGt  <  5:326409/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
AGCAGTTGGACCCGATCCATGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCCCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGCGGGTACGCCTGGAGCAAGCGGT  >  NC_002947/3915651‑3915936

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