Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,481,134 C→A 36.4% P173P (CCC→CCA ybdR → Zn‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,481,1340CA36.4% 1.1 / 7.7 11P173P (CCC→CCAybdRZn‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (5/2);  new base A (2/2);  total (7/4)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-03

CAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTC  >  NC_002947/4481002‑4481257
                                                                                                                                    |                                                                                                                           
cAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCAcc                                                                                                                          >  1:169671/1‑136 (MQ=255)
    cATCTGTATGGAGGGATCCCTGGTGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCAATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCtt                                                                                                                      <  1:177170/136‑1 (MQ=255)
                       cTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCAACCGCTTGGc                                                                                                                   >  1:46719/1‑120 (MQ=255)
                       cTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCAACCGCTTGGc                                                                                                                   >  1:48077/1‑120 (MQ=255)
                       cTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCAACCGCTTGGc                                                                                                                   <  2:46719/120‑1 (MQ=255)
                       cTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCAACCGCTTGGc                                                                                                                   <  2:48077/120‑1 (MQ=255)
                                  gCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCg                                                                                        >  2:248529/1‑136 (MQ=255)
                                                       aaaGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCt                                                                   >  1:133245/1‑136 (MQ=255)
                                                           gCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATgg                                                               <  2:169671/136‑1 (MQ=255)
                                                                                                                        cGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTc  >  1:232270/1‑136 (MQ=255)
                                                                                                                        cGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTc  >  1:91816/1‑136 (MQ=255)
                                                                                                                                    |                                                                                                                           
CAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTC  >  NC_002947/4481002‑4481257

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: