Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,691,814 T→C 46.4% intergenic (+39/‑132) 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,8140TC46.4% ‑1.3 / 11.1 15intergenic (+39/‑132)PP_1487/PP_1488hypothetical protein/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base T (5/3);  new base C (4/3);  total (9/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

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: