Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,691,815 C→G 42.7% intergenic (+40/‑131) PP_1487 → / → PP_1488 hypothetical protein/methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,691,8150CG42.7% 0.2 / 14.1 15intergenic (+40/‑131)PP_1487/PP_1488hypothetical protein/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base G (4/2);  total (9/6)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

TTCATGGAGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGA  >  NC_002947/1691674‑1691958
                                                                                                                                             |                                                                                                                                               
ttCATGGAGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGACGGAGcgccg                                                                                                                                         >  1:191646/1‑150 (MQ=255)
           tGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTgcggcg                                                                                                                              >  1:181740/1‑150 (MQ=255)
                                              tGGACCCGGGCTTGCGCCAGGAACTGGCAAGGCATGTGTGGTTGGCTGCAGCGTAGTGGCGGTGCTTCCCTTGTCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTTCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAAt                                                                                           <  1:188803/150‑1 (MQ=38)
                                                   ccGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATACACCTCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACAt                                                                                      >  1:117744/1‑150 (MQ=255)
                                                                                            tGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTc                                             >  2:104994/1‑150 (MQ=255)
                                                                                                                   cACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCATAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAgg                      >  2:287815/1‑150 (MQ=255)
                                                                                                                      gggcTATTGAGCGCCTGTGAGACCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCt                   <  2:94346/147‑1 (MQ=255)
                                                                                                                         gaattttgcgcCTGTGAGACGTAGGGCCGCACGGGCGTTGCTCGATCTCACAGGGGCAAAATTCCTGACGACAATCACATCAATACGGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCCCGTTGTGAAGTCCTTAGGGCTGaa                >  2:213956/8‑150 (MQ=21)
                                                                                                                               aGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTgaa          <  1:104994/150‑1 (MQ=255)
                                                                                                                                    cTGTGAGACGGAGCGCCGCACGGGCGGCGCTCGATCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACt     >  1:247100/1‑150 (MQ=38)
                                                                                                                                    cTGTGAGACGGAGCGCCGCACGGGCGGCGCTCGATCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACt     >  2:426127/1‑150 (MQ=38)
                                                                                                                                     tgtgAGACGGAGCGCCGCACGGGCGGCGCTCGATCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTg    <  1:426127/150‑1 (MQ=38)
                                                                                                                                     tgtgAGACGGAGCGCCGCACGGGCGGCGCTCGATCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTg    <  2:247100/150‑1 (MQ=38)
                                                                                                                                       tgAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGa  <  2:181740/150‑1 (MQ=255)
                                                                                                                                       tgAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGa  >  2:266367/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                               
TTCATGGAGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGA  >  NC_002947/1691674‑1691958

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: