| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 3,554,933 | A→C | 15.4% | K169T (AAA→ACA) | PP_3139 → | group 1 family glycosyl transferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 3,554,933 | 0 | A | C | 15.4% | 52.0 / 5.7 | 26 | K169T (AAA→ACA) | PP_3139 | group 1 family glycosyl transferase |
| Reads supporting (aligned to +/- strand): ref base A (12/10); new base C (2/2); total (14/12) | |||||||||||
| 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.46e-01 | |||||||||||
GGCGACACATGTGCACGCGCATTTCGGCACCAACTCTACTGAAGTGGTCATGCTGGCCAACGCCCTGGGCGGGCCGGCATACAGCTTTACCGTGCACGGCCCTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTGGGTGGCGCCCGAGCAGTGGGGCAAGGTGAAGGTGGTGCATTGCGGCCTGGAGCGTGGGTTTCATGAGG > NC_002947/3554790‑3555064 | ggCGACACATGTGCACGCGCATTTCGGCACCAACTCTACTGAAGTGGTCATGCTGGCCAACGCCCTGGGCGGGCCGGCATACAGCTTTACCGTGCACGGCCCTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACg < 1:370299/150‑1 (MQ=255) gCGACACATGTGCACGCGCATTTCGGCACCAACTCTACTGAAGTGGTCATGCTGGCCAACGCCCTGGGCGGGCCGGCATACAGCTTTACCGTGCACGGCCCTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACgg < 2:384580/150‑1 (MQ=255) aTGTGCACGCGCATTTCGGCACCAACTNTACTGAAGTGGTCATGCTGGCCAACGCCCTGGGCGGGCCGGCATACAGCTTTACCGTGCACGGCCCTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGcgccg > 1:499366/1‑150 (MQ=255) aTGTGCACGCGCATTTCGGCACCAACTCTACTGAAGTGGTCATGCTGGCCAACGCCCTGGGCGGGCCGGCATACAGCTTTACCGTGCACGGCCCTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGcgccg > 1:413485/1‑150 (MQ=255) gCACGCGCATTTCGGCACCAACTCTACTGAAGTGGTCATGCTGGCCAACGCCCTGGGCGGGCCGGCATACAGCTTTACCGTGCACGGCCCTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCtt < 2:413485/150‑1 (MQ=255) aTTTCGGCACCAACTCTACTGAAGTGGTCATGCTGGCCAACGCCCTGGGCGGGCCGGCATACAGCTTTACCGTGCACGGCCCTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCgcgg > 1:120621/1‑150 (MQ=255) tACTGAAGTGGTCATGCTGGCCAACGCCCTGGGCGGGCCGGCATACAGCTTTACCGTGCACGGCCCTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACgg < 2:499366/150‑1 (MQ=255) tGAAGTGGTCATGCTGGCCAACGCCCTGGGCGGGCCGGCATACAGCTTTACCGTGCACGGCCCTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGgcg > 1:419804/1‑150 (MQ=255) agccagCCCTGGGCGGGCCGGCATACAGCTTTACCGTGCACGGCCCTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGtt < 2:483597/145‑1 (MQ=255) aCGCCCTGGGCGGGCCGGCATACAGCTTTACCGTGCACGGCACTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTggg > 1:417247/1‑150 (MQ=255) gCCCTGGGCGGGCCGGCATACAGCTTTACCGTGCACGGCCCTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCATTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTGGGTg > 1:502534/1‑150 (MQ=255) ggcgggcCGGCATACAGCTTTACCGTGCACGGCCCTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGACAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTGGGTGGCGccc > 2:419119/1‑150 (MQ=255) cGGCATACAGCTTTACCGTGCACGGCACTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTGGGTGGCGCCCGAGCAGt < 2:417247/150‑1 (MQ=255) cATACAGCTTTACCGTGCACGGCCCTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGACAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTGGGTGGCGCCCGAGCAGTggg < 1:419119/150‑1 (MQ=255) gaagaaTTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGc < 2:466590/90‑1 (MQ=255) gaagaaTTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGc > 1:466590/1‑90 (MQ=255) aaTTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTGGGTGGCGCCCGAGCAGTGGGGCaag > 1:426178/1‑125 (MQ=255) aaTTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTGGGTGGCGCCCGAGCAGTGGGGCaag < 2:426178/125‑1 (MQ=255) ttCGACAAACCACAGTTCCTGCACTTGGGTGAGAAGGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTGGGTGGCGCCCGAGCAGTGGGGCAAGGTGAAGGTGGTGCATTGCGGCCTGGAg > 2:95669/1‑150 (MQ=255) aaaCCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTGGGTGGCGCCCGAGCAGTGGGGCAAGGTGAAGGTGGTGCATTGCGGCCTGGAGCGTggg > 2:466050/1‑150 (MQ=255) ccACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCATTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTGGGTGGCGCCCGAGCAGTGGGGCAAGGTGAAGGTGGTGCATTGCGGCCTGGAGCGTGGGttt < 2:502534/150‑1 (MQ=255) cacaGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTGGGTGGCGCCCGAGCAGTGGGGCAAGGTGAAGGTGGTGCATTGCGGCCTGGAGCGTGGGTTTc > 2:876/1‑150 (MQ=255) acaGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTGGGTGGCGCCCGAGCAGTGGGGCAAGGTGAAGGTGGTGCATTGCGGCCTGGAGCGTGGGTTTCa > 1:22240/1‑150 (MQ=255) tCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTGGGTGGCGCCCGAGCAGTGGGGCAAGGTGAAGGTGGTGCATTGCGGCCTGGAGCGTGGGTTTCATGAgg < 2:22240/150‑1 (MQ=255) ttGGGTGAGACAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTGGGTGGCGCCCGAGCAGTGGGGCAAGGTGAAGGTGGTGCa > 2:164458/1‑113 (MQ=255) ttGGGTGAGACAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTGGGTGGCGCCCGAGCAGTGGGGCAAGGTGAAGGTGGTGCa < 1:164458/113‑1 (MQ=255) | GGCGACACATGTGCACGCGCATTTCGGCACCAACTCTACTGAAGTGGTCATGCTGGCCAACGCCCTGGGCGGGCCGGCATACAGCTTTACCGTGCACGGCCCTGAAGAATTCGACAAACCACAGTTCCTGCACTTGGGTGAGAAAGTACGGCGCGCCGCCTTTGTCGCGGCGGTGAGCTCGTACGGGCGCAGCCAGTTGTACCGTTGGGTGGCGCCCGAGCAGTGGGGCAAGGTGAAGGTGGTGCATTGCGGCCTGGAGCGTGGGTTTCATGAGG > NC_002947/3554790‑3555064 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |