Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 858,206 A→G 30.8% A207A (GCT→GCC phr ← deoxyribodipyrimidine photolyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947858,2060AG30.8% 11.4 / 5.9 13A207A (GCT→GCCphrdeoxyribodipyrimidine photolyase
Reads supporting (aligned to +/- strand):  ref base A (2/7);  new base G (2/2);  total (4/9)
Fisher's exact test for biased strand distribution p-value = 5.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGGGCGAAATCACGCCGGCGGCCAGGTAGGGGGAAAGCTGGCTGGTACCAGGCTTGGCAGGCAAGTCGCGCAAGTGGTGGTAGTCGTCGATGGTTTCGTCCAGGAAGCGTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTG  >  NC_002947/858074‑858341
                                                                                                                                    |                                                                                                                                       
gCGGGGCGAAATCACGCCGGCGGCCAGGTAGGGGGAAAGCTGGCTGGTACCAGGCTTGGCAGGCAAGTCGCGCAAGTGGTGGTAGTCGTCGATGGTTTCGTCCAGGAAGCGTGACAGGCGGGACTGTGCCTCGGCt                                                                                                                                      >  2:22813/1‑136 (MQ=255)
gCGGGGCGAAATCACGCCGGCGGCCAGGTAGGGGGAAAGCTGGCTGGTACCAGGCTTGGCAGGCAAGTCGCGCAAGTGGTGGTAGTCGTCGATGGTTTCGTCCAGGAAGCGTGACAGGCGGGACTGTGCCTCGGCt                                                                                                                                      >  2:25457/1‑136 (MQ=255)
       gAAATCACGCCGGCGGCCAGGTAGGGGGAAAGCTGGCTGGTACCAGGCTTGGCAGGCAAGTCGCGCAAGTGGTGGTAGTCGTCGATGGTTTCGTCCAGGAAGCGTGACAGGCGGGACTGTGCCTCAGCTTCACCgg                                                                                                                               <  2:515126/136‑1 (MQ=255)
                            tAGGGGGAAAGCTGGCTGGTACCAGGCTTGGCAGGCAAGTCGCGCAAGTGGTGGTAGTCGTCGATGGTTTCGTCCAGGAAGCGTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATg                                                                                                          >  2:519621/1‑136 (MQ=255)
                                        tggctggTACCAGGCTTGGCAGGCAAGTCGCGCAAGTGGTGGTAGTCGTCGATGGTTTCGTCCAGGAAGCGTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGttt                                                                                              >  1:87106/1‑136 (MQ=255)
                                                      cTTGGCAGGCAAGTCGCGCAAGTGGTGGTAGTCGTCGATGGTTTCGTCCAGGAAGCGTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCctct                                                                                <  1:266798/136‑1 (MQ=255)
                                                           caggcaAGTCGCGCAAGTGGTGGTAGTCGTCGATGGTTTCGTCCAGGAAGCGTGACAGGCGGGACTGTGCCTCGGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATg                                                                           <  1:22813/136‑1 (MQ=255)
                                                           caggcaAGTCGCGCAAGTGGTGGTAGTCGTCGATGGTTTCGTCCAGGAAGCGTGACAGGCGGGACTGTGCCTCGGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATg                                                                           <  1:25457/136‑1 (MQ=255)
                                                           caggcaAGTCGCGCAAGTGGTGGTAGTCGTCGATGGTTTCGTCCAGGAAGCGTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATg                                                                           <  2:505497/136‑1 (MQ=255)
                                                                                                                      cGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCa                <  2:14461/136‑1 (MQ=255)
                                                                                                                          aCTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAgggc            <  1:68162/136‑1 (MQ=255)
                                                                                                                                 cTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAgg     <  1:443236/136‑1 (MQ=255)
                                                                                                                                 cTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCAAGGGCGGGCAgg     <  2:87106/136‑1 (MQ=255)
                                                                                                                                    aGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTg  <  2:39043/136‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                       
GCGGGGCGAAATCACGCCGGCGGCCAGGTAGGGGGAAAGCTGGCTGGTACCAGGCTTGGCAGGCAAGTCGCGCAAGTGGTGGTAGTCGTCGATGGTTTCGTCCAGGAAGCGTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTG  >  NC_002947/858074‑858341

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