Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,067,659 C→A 26.6% A15E (GCG→GAG)  PP_4460 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,067,6590CA26.6% 28.9 / 4.3 15A15E (GCG→GAG) PP_4460LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (5/6);  new base A (2/2);  total (7/8)
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.27e-01

AAGCTTTTGATTAGTAGAAGGAGCTCCCGTAAGCTAAAAAGATCATGTCTTTATGCTATCGGCAGGTGGAGTACCATGGAGAACAGAATCACCCTTCGTCACCTAGAAGCGTTTCGGGCGATCATGGTGCGTAAATCCGTCACCGGTGCTGCAGAAATGCTTGAAGTGACCCAACCGGTAGTCACGAGGTTAATTGCCGATCTTGAGGAACGAATAGCACTGACACTGTTTACAAGAACTAAGGGGCGCCTGGT  >  NC_002947/5067541‑5067794
                                                                                                                      |                                                                                                                                       
aaGCTTTTGATTAGTAGAAGGAGCTCCCGTAAGCTAAAAAGATCATGTCTTTATGCTATCGGCAGGTGGAGTACCATGGAGAACAGAATCACCCTTCGTCACCTAGAAGCGTTTCGGGCGATCATGGTGCGTAAAt                                                                                                                        >  1:62936/1‑136 (MQ=255)
           tagtagAAGGAGCTCCCGTAAGCTAAAAAGATCATGTCTTTATGCTATCGGCAGGTGGAGTACCATGGAGAACAGAATCACCCTTCGTCACCTAGAAGCGTTTCGGGCGATCATGGTGCGTAAATCCGTCACCggg                                                                                                             >  2:640165/1‑135 (MQ=255)
                gAAGGAGCTCCCGTAAGCTAAAAAGATCATGTCTTTATGCTATCGGCAGGTGGAGTACCATGGAGAACAGAATCACCCTTCGTCACCTAGAAGCGTTTCGGGCGATCATGGTGCGTAAATCCGTCACCGGtgctgc                                                                                                        <  1:65746/136‑1 (MQ=255)
                gAAGGAGCTCCCGTAAGCTAAAAAGATCATGTCTTTATGCTATCGGCAGGTGGAGTACCATGGAGAACAGAATCACCCTTCGTCACCTAGAAGCGTTTCGGGCGATCATGGTGCGTAAATCCGTCACCGGtgctgc                                                                                                        <  2:373675/136‑1 (MQ=255)
                                               tCTTTATGCTATCGGCAGGTGGAGTACCATGGAGAACAGAATCACCCTTCGTCACCTAGAAGCGTTTCGGGCGATCATGGTGCGTAAATCCGTCACCGGTGCTGCAGAAATGCTTGAAGTGACCCAACCGGTAGTc                                                                         <  1:32539/136‑1 (MQ=255)
                                                             gCAGGTGGAGTACCATGGAGAACAGAATCACCCTTCGTCACCTAGAAGCGTTTCGGGAGATCATGGTGCGTAAATCCGTCACCGGTGCTGCAGAAATGCTTGAAGTGACCCAACCGGTAGTCACGAGGTTAATTGc                                                           >  1:16853/1‑136 (MQ=255)
                                                             gCAGGTGGAGTACCATGGAGAACAGAATCACCCTTCGTCACCTAGAAGCGTTTCGGGAGATCATGGTGCGTAAATCCGTCACCGGTGCTGCAGAAATGCTTGAAGTGACCCAACCGGTAGTCACGAGGTTAATTGc                                                           >  1:61545/1‑136 (MQ=255)
                                                                               agaACAGAATCACCCTTCGTCACCTAGAAGCGTTTCGGGCGATCATGGAGCGTAAATCCGTCACCGGTGCTGCAGAAATGCTTGAAGTGACCCAACCGGTAGTCACGAGGTTAATTGCCGATCTTGAGGAACGAAt                                         >  1:475460/1‑136 (MQ=255)
                                                                                      aaTCACCCTTCGTCACCTAGAAGCGTTTCGGGCGATCATGGTGCGTAAATCCGTCACCGGTGCTGCAGAAATGCTTGAAGTGACCCAACCGGTAGTCACGAGGTTAATTGCCGATCTTGAGGAACGAATAGCACTg                                  <  1:504677/136‑1 (MQ=255)
                                                                                              ttCGTCACCTAGAAGCGTTTCGGGCGATCATGGTGCGTAAATCCGTCACCGGTGCTGCAGAAATGCTTGAAGTGACCCAACCGGTAGTCACGAGGTTAATTGCCGATCTTGAGGAACGAATAGCACTGACACTGtt                          <  1:258954/136‑1 (MQ=255)
                                                                                                  tCACCTAGAAGCGTTTCGGGCGATCATGGTGCGTAAATCCGTCACCGGTGCTGCAGAAATGCTTGAAGTGACCCAACCGGTAGTCACGAGGTTAATTGCCGATCTTGAGGAACGAATAGCACTGACACTGTTTACa                      <  1:665641/136‑1 (MQ=255)
                                                                                                         gAAGCGTTTCGGGCGATCATGGTGCGTAAATCCGTCACCGGTGCTGCAGAAATGCTTGAAGTGACCCAACCGGTAGTCACGAGGTTATTTGCCGATCTTGAGGAACGAATAGCACTGACACTGTTTACAAGAACTa               >  2:426377/1‑136 (MQ=255)
                                                                                                            gCGTTTCGGGCGATCATGGTGCGTAAATCCGTCACCGGTGCTGCAGAAATGCTTGAAGTGACCCAACCGGTAGTCACGAGGTTAATTGCCGATCTTGAGGAACGAATAGCACTGACACTGTTTACAAGAACTAAgg            >  2:408780/1‑136 (MQ=255)
                                                                                                               tttCGGGAGATCATGGTGCGTAAATCCGTCACCGGTGCTGCAGAAATGCTTGAAGTGACCCAACCGGTAGTCACGAGGTTAATTGCCGATCTTGAGGAACGAATAGCACTGACACTGTTTACAAGAACTAAGGGgc         <  2:16853/136‑1 (MQ=255)
                                                                                                               tttCGGGAGATCATGGTGCGTAAATCCGTCACCGGTGCTGCAGAAATGCTTGAAGTGACCCAACCGGTAGTCACGAGGTTAATTGCCGATCTTGAGGAACGAATAGCACTGACACTGTTTACAAGAACTAAGGGgc         <  2:61545/136‑1 (MQ=255)
                                                                                                                      cGATCATGGTGCGTAAATCCGTCACCGGTGCTGCAGAAATGCTTGAAGTGACCCAACCGGTAGTCACGAGGTTAATTGCCGATCTTGAGGAACGAATAGCACTGACACTGTTTACAAGAACTAAGGGGCGCCTGGt  <  1:379462/136‑1 (MQ=255)
                                                                                                                      |                                                                                                                                       
AAGCTTTTGATTAGTAGAAGGAGCTCCCGTAAGCTAAAAAGATCATGTCTTTATGCTATCGGCAGGTGGAGTACCATGGAGAACAGAATCACCCTTCGTCACCTAGAAGCGTTTCGGGCGATCATGGTGCGTAAATCCGTCACCGGTGCTGCAGAAATGCTTGAAGTGACCCAACCGGTAGTCACGAGGTTAATTGCCGATCTTGAGGAACGAATAGCACTGACACTGTTTACAAGAACTAAGGGGCGCCTGGT  >  NC_002947/5067541‑5067794

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: