Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,022,806 C→G 11.8% intergenic (‑501/+167) PP_3547 ← / ← emrB NAD(P)‑binding oxidoreductase/multidrug transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,022,8060CG11.8% 106.4 / 6.7 42intergenic (‑501/+167)PP_3547/emrBNAD(P)‑binding oxidoreductase/multidrug transporter permease
Reads supporting (aligned to +/- strand):  ref base C (18/19);  new base G (3/2);  total (21/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.98e-01

TTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCG  >  NC_002947/4022659‑4022951
                                                                                                                                                   |                                                                                                                                                 
ttGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATc                                                                                                                                                   >  4:53082/1‑148 (MQ=255)
   tCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATg                                                                                                                                               >  1:65469/1‑149 (MQ=255)
     aaaCGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGc                                                                                                                                              <  8:24386/148‑1 (MQ=255)
        cgcgACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGt                                                                                                                                          <  3:56411/149‑1 (MQ=255)
          cgACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCa                                                                                                                                        <  5:94325/149‑1 (MQ=255)
             cAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATGGATGCCCGTCaagc                                                                                                                                     <  7:240168/149‑4 (MQ=255)
              aaaaaGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATGGATGCCCGTCaagc                                                                                                                                     >  8:240168/1‑145 (MQ=255)
              aaaaaGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTc                                                                                                                                    <  2:65469/149‑1 (MQ=255)
              aaaaaGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTc                                                                                                                                    >  5:207442/1‑149 (MQ=255)
                   gTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTc                                                                                                                               <  3:53082/149‑1 (MQ=255)
                   gTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCCCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTc                                                                                                                               >  6:224518/1‑149 (MQ=255)
                        cAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATGGATGCCCGTc                                                                                                                                         <  4:42605/134‑1 (MQ=255)
                        cAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATGGATGCCCGTc                                                                                                                                         >  3:42605/1‑134 (MQ=255)
                                    gCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTc                                                                                                                <  6:207442/147‑1 (MQ=37)
                                    gCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCgg                                                                                                              <  6:47418/149‑1 (MQ=37)
                                         tATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTACCGGGAGGACCGCGCTGCCGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATGCCTGCCCGTCCCCTCCCTTCGCAACCCCAGCGGTCGgccgcc                                                                                                         >  7:21963/1‑149 (MQ=11)
                                                           cTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCgcgg                                                                                       >  8:74325/1‑149 (MQ=38)
                                                            tACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCt                                                                                                                      >  8:174945/1‑117 (MQ=35)
                                                            tACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCt                                                                                                                      <  7:174945/117‑1 (MQ=35)
                                                                 gccGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTg                                                                                 <  8:230561/148‑1 (MQ=39)
                                                                   cGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTgg                                                                                <  8:29513/148‑1 (MQ=40)
                                                                           ccTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATcgcg                                                                       >  1:256595/1‑149 (MQ=255)
                                                                            cTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATcgcgc                                                                      >  3:238011/1‑149 (MQ=255)
                                                                            cTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATcgcgc                                                                      <  7:220246/149‑1 (MQ=255)
                                                                            cTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCAGCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATcgcgc                                                                      >  5:91760/1‑149 (MQ=255)
                                                                                 tCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCt                                                                   >  6:156080/1‑147 (MQ=255)
                                                                                       cTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCa                                                           <  2:256595/149‑1 (MQ=255)
                                                                                              ccGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCa                                                     >  2:181792/1‑148 (MQ=255)
                                                                                                gggAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCaaaa                                                  <  5:253587/149‑1 (MQ=255)
                                                                                                 ggaggaCCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAg                                                 <  5:224518/149‑1 (MQ=255)
                                                                                                        cgcgCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagc                                          <  7:74325/149‑1 (MQ=255)
                                                                                                             gctgggCCCCAGGATGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagag                                     >  1:34209/5‑149 (MQ=255)
                                                                                                               cAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagc                                    >  2:227556/1‑148 (MQ=255)
                                                                                                                  gCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAAt                                                       >  6:161429/1‑126 (MQ=255)
                                                                                                                  gCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAAt                                                       <  5:161429/126‑1 (MQ=255)
                                                                                                                      cATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagc                              >  3:243572/1‑147 (MQ=255)
                                                                                                                                 ccAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgagagcga                      <  1:223757/144‑1 (MQ=255)
                                                                                                                                 ccAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgagagcga                      >  2:223757/1‑144 (MQ=255)
                                                                                                                                   aGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgagagcgagagcgag               <  7:36845/149‑1 (MQ=255)
                                                                                                                                          gACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgagagcgagagcgagagcga          >  2:215575/1‑147 (MQ=255)
                                                                                                                                           aCGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAgcgagagcgagagcgagagcgagagcgagagcgaga        <  5:156080/148‑1 (MQ=255)
                                                                                                                                             gggCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCa     >  3:171429/1‑149 (MQ=255)
                                                                                                                                                cATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCg  >  2:155700/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCG  >  NC_002947/4022659‑4022951

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