Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,318,934 G→C 9.8% G399R (GGC→CGC)  PP_3789 → efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,318,9340GC 9.8% 109.6 / 4.3 41G399R (GGC→CGC) PP_3789efflux transporter
Reads supporting (aligned to +/- strand):  ref base G (21/16);  new base C (2/2);  total (23/18)
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.28e-01

ATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCG  >  NC_002947/4318788‑4319079
                                                                                                                                                   |                                                                                                                                                 
aTCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGc                                                                                                                                                 <  6:192019/149‑1 (MQ=255)
  cTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGCGCa                                                                                                                                                >  3:178192/1‑148 (MQ=255)
             gggggCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCgg                                                                                                                                     <  3:177239/148‑1 (MQ=255)
               gggCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGCGCAAAGTTGCCCggggt                                                                                                                                  <  4:178192/149‑2 (MQ=255)
                 gcggcAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCggg                                                                                                                                    >  6:204300/1‑145 (MQ=255)
                          aGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAAt                                                                                                                        <  6:179527/148‑1 (MQ=255)
                             ccGGCGTGGTGTGGATTCATACTATCGCTCGCGGGACCCCAGGGTTTGTGCTTCTCACTGATCGCATGCGTGCTGGCGCGTAGCAGCCATTGTGACAAGCGCACTCGACTGACACCGGGCAACttt                                                                                                                                           >  7:104714/1‑126 (MQ=255)
                                           aTTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTg                                                                                                      <  5:204300/149‑1 (MQ=255)
                                               aTACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCTGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGt                                                                                                                <  1:65285/136‑1 (MQ=255)
                                               aTACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGt                                                                                                                >  2:65285/1‑135 (MQ=255)
                                                    aTCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAAt                                                                                             >  7:29661/1‑149 (MQ=255)
                                                            cgGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCAcc                                                                                     >  8:228284/1‑149 (MQ=255)
                                                                     aGGGTTTGTGCTTCTCACTGATCGCCTGGGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGt                                                                                           <  1:298401/134‑1 (MQ=255)
                                                                     aGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGt                                                                                           >  2:298401/1‑134 (MQ=255)
                                                                       ggTTTGTGGTTCTCACTGATCGCCTGCGTGCGGGCGGGTAGCAGCGTTTGTGCCATGCGCCCTCGACTGACCGCGGGCTACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTtctc                                                                           <  2:290436/148‑1 (MQ=255)
                                                                            gtgCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGtctc                                                                                        <  7:207886/130‑1 (MQ=255)
                                                                            gtgCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGtctc                                                                                        >  8:207886/1‑130 (MQ=255)
                                                                              gCTTCTCACTGATCGCCTGCGTGCTGGCGGATAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCAtt                                                                   >  1:52919/1‑149 (MQ=255)
                                                                                ttCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCt                                                                 >  3:279968/1‑149 (MQ=255)
                                                                                 tctcACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGCGCAAAGTTGCCCggggt                                                                                                                                  >  2:259975/1‑82 (MQ=38)
                                                                                 tctcACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGCGCAAAGTTGCCCggggt                                                                                                                                  <  1:259975/83‑2 (MQ=38)
                                                                                     aCTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCAccc                                                                                    >  8:720/1‑125 (MQ=255)
                                                                                     aCTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCAccc                                                                                    <  7:720/125‑1 (MQ=255)
                                                                                          tCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAg                                                       >  8:30153/1‑149 (MQ=255)
                                                                                            gCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGTGTTATTTAATAATACAGAAGGGGCTTTGGTCGTCAATGTCTCACCGTCAGTTGTGTGATCATTATGGTGCCTTAGc                                                      >  6:178184/1‑148 (MQ=255)
                                                                                              cTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTaa                                                   >  5:246830/1‑149 (MQ=255)
                                                                                              cTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTaa                                                   >  2:23030/1‑149 (MQ=255)
                                                                                                     cTGGCGGATAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGt                                             <  2:52919/148‑1 (MQ=255)
                                                                                                      tGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGttt                                           >  4:5888/1‑149 (MQ=255)
                                                                                                       ggCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACGCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGttt                                           >  6:265519/1‑148 (MQ=255)
                                                                                                        gCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCt                                                                                   <  5:39318/107‑1 (MQ=255)
                                                                                                        gCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCt                                                                                   >  6:39318/1‑107 (MQ=255)
                                                                                                             tAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAAt                                    <  7:228284/149‑1 (MQ=255)
                                                                                                                cagcCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATAcc                                 <  8:29661/149‑1 (MQ=255)
                                                                                                                   ccTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATg                              <  6:246830/149‑1 (MQ=255)
                                                                                                                           ccATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTgg                       >  2:53646/1‑148 (MQ=255)
                                                                                                                             aTGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGaaa                    >  3:217785/1‑149 (MQ=255)
                                                                                                                              tGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAAt                   <  2:281928/149‑1 (MQ=255)
                                                                                                                                       gACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCtt          >  8:274166/1‑149 (MQ=255)
                                                                                                                                       gACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCtt          <  1:23030/149‑1 (MQ=255)
                                                                                                                                               ccGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCg  >  8:174155/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCG  >  NC_002947/4318788‑4319079

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: