Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,223,266 G→T 22.2% intergenic (‑530/+30) PP_2823 ← / ← PP_2824 methyl‑accepting chemotaxis transducer/TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,223,2660GT22.2% 32.9 / 5.3 18intergenic (‑530/+30)PP_2823/PP_2824methyl‑accepting chemotaxis transducer/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (10/4);  new base T (2/2);  total (12/6)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

GACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGAGCCAG  >  NC_002947/3223118‑3223407
                                                                                                                                                    |                                                                                                                                             
gACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACTGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAgg                                                                                                                                              <  2:498852/150‑1 (MQ=255)
    cAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAg                                                                                                                                          <  2:602999/150‑1 (MQ=255)
      aTGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGgc                                                                                                                                        >  1:269258/1‑150 (MQ=255)
       tGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGgcg                                                                                                                                       >  1:122183/1‑150 (MQ=255)
          cgcCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCg                                                                                                                                    <  2:269258/150‑1 (MQ=255)
                    tGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCgg                                                                                                                          >  2:533059/1‑150 (MQ=255)
                                                            ggTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCa                                                                                  <  1:533059/150‑1 (MQ=255)
                                                                   tCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAAt                                                                           >  1:501646/1‑150 (MQ=255)
                                                                                                  gacgGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTGAGGACTTGGCGGTTATGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCgcgccag                                            >  1:64693/1‑147 (MQ=255)
                                                                                                                     acGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACt                         >  2:198603/1‑150 (MQ=255)
                                                                                                                                 gtggGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGGTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGgcagca             >  1:96944/1‑150 (MQ=255)
                                                                                                                                         ggcGCCTTCCCTGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGAc                          >  1:484951/3‑129 (MQ=255)
                                                                                                                                         ggcGCCTTCCCTGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGAc                          <  2:484951/127‑1 (MQ=255)
                                                                                                                                          gcGCCTTCCCTGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAg                                                           >  2:444569/2‑95 (MQ=255)
                                                                                                                                          gcGCCTTCCCTGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAg                                                           <  1:444569/94‑1 (MQ=255)
                                                                                                                                            gCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGGAGTGCAGCATCATGagccag  >  2:203389/1‑150 (MQ=255)
                                                                                                                                            gCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGCTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGagccag  >  2:482506/1‑150 (MQ=255)
                                                                                                                                               ttCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTaa                                     >  1:504168/1‑112 (MQ=255)
                                                                                                                                               ttCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTaa                                     <  2:504168/112‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGAGCCAG  >  NC_002947/3223118‑3223407

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