Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,883,770 C→A 18.1% R167S (CGC→AGC)  pucL → 2‑oxo‑4‑hydroxy‑4‑carboxy‑5‑ureidoimidazoline decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,883,7700CA18.1% 60.6 / 3.2 22R167S (CGC→AGC) pucL2‑oxo‑4‑hydroxy‑4‑carboxy‑5‑ureidoimidazoline decarboxylase
Reads supporting (aligned to +/- strand):  ref base C (6/12);  new base A (2/2);  total (8/14)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

CAGTTCCCGTTCATCATGGCGGTGAAGGGCAGCAACCGGCACCAGATTCTCGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATGCGCACCCTTATGATCGAGCCCTTGACCAAAGAAGCCTTCGCCCAGTTCGGAGACGTGA  >  NC_002947/4883635‑4883897
                                                                                                                                       |                                                                                                                               
cAGTTCCCGTTCATCATGGCGGTGAAGGGCAGCAACCGGCACCAGATTCTCGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTcc                                                                                                                                 <  1:13947/136‑1 (MQ=255)
     cccGTTCATCATGGCGGTGAAGGGCAGCAACCGGCACCAGATTCTCGCCACCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTtg                                                                                                                            <  1:95831/136‑1 (MQ=255)
         ttCATCATGGCGGTGAAGGGCAGCAACCGGCACCAGATTCTCGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTtgctgc                                                                                                                        >  2:304694/1‑136 (MQ=255)
          tcatcaTGGCGGTGAAGGGCAGCAACCGGCACCAGATTCTCGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCa                                                                                                                       >  1:151188/1‑136 (MQ=255)
          tcatcaTGGCGGTGAAGGGCAGCAACCGGCACCAGATTCTCGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCa                                                                                                                       >  1:457548/1‑136 (MQ=255)
            atcatGGCGGTGAAGGGCAGCAACCGGCACCAGATTCTCGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGc                                                                                                                     <  1:263051/136‑1 (MQ=255)
                 ggCGGTGAAGGGCAGCAACCGGCACCAGATTCTCGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTaa                                                                                                                <  1:33645/136‑1 (MQ=255)
                                       cACCAGATTCTCGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAg                                                                                          <  1:106905/136‑1 (MQ=255)
                                          cAGATTCTCGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGaaa                                                                                       <  2:50482/136‑1 (MQ=255)
                                               tctcGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGaa                                                                                  <  1:304694/136‑1 (MQ=255)
                                                   gccgccTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAAt                                                                              <  2:357438/136‑1 (MQ=255)
                                                                    cATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCAt                                                             <  1:215105/136‑1 (MQ=255)
                                                                     aTCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATg                                                            <  1:427970/136‑1 (MQ=255)
                                                                     aTCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATg                                                            <  1:238488/136‑1 (MQ=255)
                                                                     aTCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATg                                                            <  2:4804/136‑1 (MQ=255)
                                                                                tgccgatgccgaGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCAGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATg                                                            <  1:449990/125‑1 (MQ=255)
                                                                                tgccgatgccgaGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCAGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATg                                                            <  1:452102/125‑1 (MQ=255)
                                                                                tgccgatgccgaGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCAGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATg                                                            >  2:449990/1‑125 (MQ=255)
                                                                                tgccgatgccgaGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCAGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATg                                                            >  2:452102/1‑125 (MQ=255)
                                                                                   cgatgccgaGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATGCGCACCCTTATGAt                                              >  1:215036/1‑136 (MQ=255)
                                                                                                                             cGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGATCAATAGAAGAGATAACCGCATGCGCACCCTTATGATCGAGCCCTTGACCAAAGAAGCCTTCGCCCAGTTCGGAGACGt    <  2:476368/136‑1 (MQ=255)
                                                                                                                               cccTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATGCGCACCCTTATGATCGAGCCCTTGACCAAAGAAGCCTTCGCCCAGTTCGGAGACGTGa  >  1:481714/1‑136 (MQ=255)
                                                                                                                               cccTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATGCGCACCCTTATGATCGAGCCCTTGACCAAAGAAGCCTTCGCCCAGTTCGGAGACGTGa  >  2:321999/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                               
CAGTTCCCGTTCATCATGGCGGTGAAGGGCAGCAACCGGCACCAGATTCTCGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATGCGCACCCTTATGATCGAGCCCTTGACCAAAGAAGCCTTCGCCCAGTTCGGAGACGTGA  >  NC_002947/4883635‑4883897

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: