Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,538,571 G→C 20.9% intergenic (+61/+60) PP_4871 → / ← PP_4872 lysine decarboxylase family protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,538,5710GC20.9% 46.6 / 9.5 24intergenic (+61/+60)PP_4871/PP_4872lysine decarboxylase family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (12/7);  new base C (3/2);  total (15/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

GGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAG  >  NC_002947/5538424‑5538707
                                                                                                                                                   |                                                                                                                                        
gggCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAgg                                                                                                                                          <  5:2390/148‑1 (MQ=255)
   cagcagCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAAt                                                                                                                                      <  4:242084/149‑1 (MQ=255)
    agcagcCGGACGAGCTGCTGGAAGGGATGGTCAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAAt                                                                                                                                      >  1:31656/1‑148 (MQ=255)
    agcagcCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAAt                                                                                                                                      <  2:31656/148‑1 (MQ=255)
        gcCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATggggg                                                                                                                                 <  1:23046/149‑1 (MQ=255)
                  ctgctgGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCg                                                                                                                             >  8:235739/1‑143 (MQ=255)
                  ctgctgGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCg                                                                                                                             <  7:235739/143‑1 (MQ=255)
                    gctgGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCCCATTGCCTTTGTGGGAGCGGCCtt                                                                                                                        <  8:180219/146‑1 (MQ=38)
                    gctgGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCCCATTGCCTTTGTGGGAGCGGCCtt                                                                                                                        >  7:180219/1‑146 (MQ=38)
                                            ggCGCCGGCCGGGCCTATTGGGGTCCACACGCCGCCTCAGTGACCCTGCCATTTTTGGGACCGCTGCGCGGTCCTTTGGCGACAAACGGCGGCGCCCACAAAGGCAATGGGGTTGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGccc                                                                                             <  8:255466/149‑1 (MQ=16)
                                                       cgcCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGGCACAAGGCCCCTGCCACAAAGGCAATGGGGGGGCGGCCGTGTGTCGCGATGGGCGGGAAATCCGGCCCAATATTGGAt                                                                                  >  3:175584/1‑149 (MQ=255)
                                                                                       aGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAagca                                                                                                 >  4:97451/1‑102 (MQ=255)
                                                                                                                       ttCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCctct                  >  6:174839/1‑149 (MQ=255)
                                                                                                                          gcgACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCttt                <  5:174839/148‑1 (MQ=255)
                                                                                                                            gACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTgcg             >  1:86932/1‑149 (MQ=255)
                                                                                                                             acacAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTgcgc            >  1:48102/1‑149 (MQ=255)
                                                                                                                               acaAGGCCGCTCCCCCATTGCCTTTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCa                                                                                <  5:93729/79‑1 (MQ=14)
                                                                                                                               acaAGGCCGCTCCCCCATTGCCTTTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCa                                                                                >  6:93729/1‑79 (MQ=14)
                                                                                                                               acaAGGCCGCTCCCCCATTGCCTTTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTc          >  6:182653/1‑149 (MQ=38)
                                                                                                                                 aaGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCt        >  3:326384/1‑149 (MQ=255)
                                                                                                                                    gCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGAGCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTgg     >  5:137208/1‑149 (MQ=255)
                                                                                                                                     ccGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTgg     >  1:9780/1‑148 (MQ=255)
                                                                                                                                       gCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGAAg  >  2:248974/1‑149 (MQ=255)
                                                                                                                                         tCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGc                                   <  4:150169/114‑1 (MQ=255)
                                                                                                                                         tCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGc                                   >  3:150169/1‑114 (MQ=255)
                                                                                                                                                   |                                                                                                                                        
GGGCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAG  >  NC_002947/5538424‑5538707

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