Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,297,646 | C→A | 23.7% | intergenic (‑11/‑181) | PP_2900 ← / → pvdQ | hypothetical protein/acyl‑homoserine lactone acylase PvdQ |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,297,646 | 0 | C | A | 23.7% | 21.4 / 4.9 | 17 | intergenic (‑11/‑181) | PP_2900/pvdQ | hypothetical protein/acyl‑homoserine lactone acylase PvdQ |
Reads supporting (aligned to +/- strand): ref base C (4/9); new base A (2/2); total (6/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.84e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01 |
GGTGCGCAGGCCGAAGGTGCGCAACAGCTCGCGGATCTGTTCGTTGCCGATTCGTTCTTCAGCAAACACGTCCTGTCCCTTGGACGTTTTGCCTGGGACGCGCCCGGCCTGGACCGAGGTTTGCGGGTTGTACATGCGAATGCTCCTTTTCGTGGACGAACCCGAAATGTGGGTGGCTCACTTACTAGGACGAATGAGCACGCAAAAACTGCAACCCATTCGTGAAAAAAAATGCCTGTTTCTCATCAAGACGTTCCATCGCCTT > NC_002947/3297501‑3297765 | ggTGCGCAGGCCGAAGGTGCGCAACAGCTCGCGGATCTGTTCGTTGCCGATTCGTTCTTCAGCAAACACGTCCTGTCCCTTGGACGTTTTGCCTGGGACGCGCCCGGCCTGGACCGAGGTTTGCGGGTTGTACATGCGAATGCTCCtttt < 2:97263/150‑1 (MQ=255) gcAACAGCTCGCGGATCTGTTCGTTGCCGATTCGTTCTTCAGCAAACACGTCCTGTCCCTTGGACGTTTTGCCTGGGACGCGCCCGGCCTGGACCGAGGTTTGCGGGTTGTACATGCGAATGCTCCTTTTCGTGGACGAACCCGAAAtgt < 2:118507/150‑1 (MQ=255) cAACAGCTCGCGGATCTGTTCGTTGCCGATTCGTTCTTCAGCAAACACGTCCTGTCCCTTGGACGTTTTGCCTGGGACGCGCCCGGCCTGGACCGAGGTTTGCGGGTTGTACATGCGAATGCTCCTTTTCGTGGACGAACCCGAAAtgtg < 2:115667/150‑1 (MQ=255) ttCGTTGCCGATTCGTTCTTCAGCAAACACGTCCTGTCCCTTGGACGTTTTGCCTGGGACGCGCCCGGCCTGGACCGAGGTTTGCGGGTTGTACATGCGAATGCTCCTTTTCGTGGACGAACCCGAAATGTGGGTGGCTCACTTACTAgg > 2:124214/1‑150 (MQ=255) tCGTTCTTCAGCAAACACGTCCTGTCCCTTGGACGTTTTGCCTGGGACGCGCCCGGCCTGGACCGAGGTTTGCGGGTTGTACATGCGAATGCTCCTTTTCGTGGACGAACCCGAAATGTGGGTGGCTCACTTACTAGGACGAATGAGCAc < 1:124214/150‑1 (MQ=255) gCAAACACGTCCTGTCCCTTGGACGTTTTGCCTGGGACGCGCCCGGCCTGGACCGAGGTTTGCGGGTTGTACATGCGAATGCTCATTTTCGTGGACGAACCCGAAATGTGGGTGGCTCACt < 1:366392/121‑1 (MQ=255) gCAAACACGTCCTGTCCCTTGGACGTTTTGCCTGGGACGCGCCCGGCCTGGACCGAGGTTTGCGGGTTGTACATGCGAATGCTCATTTTCGTGGACGAACCCGAAATGTGGGTGGCTCACt < 1:366388/121‑1 (MQ=255) gCAAACACGTCCTGTCCCTTGGACGTTTTGCCTGGGACGCGCCCGGCCTGGACCGAGGTTTGCGGGTTGTACATGCGAATGCTCATTTTCGTGGACGAACCCGAAATGTGGGTGGCTCACt > 2:366388/1‑121 (MQ=255) gCAAACACGTCCTGTCCCTTGGACGTTTTGCCTGGGACGCGCCCGGCCTGGACCGAGGTTTGCGGGTTGTACATGCGAATGCTCATTTTCGTGGACGAACCCGAAATGTGGGTGGCTCACt > 2:366392/1‑121 (MQ=255) gTCCTGTCCCTTGGACGTTTTGCCTGGGACGCGCCCGGCCTGGACCGAGGTTTGCGGGTTGTACATGCGAATGCTCCTTTTCGTGGACGAACCCGAAATgtgg > 1:303555/1‑103 (MQ=255) gTCCTGTCCCTTGGACGTTTTGCCTGGGACGCGCCCGGCCTGGACCGAGGTTTGCGGGTTGTACATGCGAATGCTCCTTTTCGTGGACGAACCCGAAATgtgg < 2:303555/103‑1 (MQ=255) tCCTGTCCCTTGGACGTTTTGCCTGGGACGCGCCCGGCCTGGAACGAGGTTTGCGGGTTGTTCACGCGAATGCTCCTTTTCTTGGACGATCCCGAAATGTGGGTGGCTCACTTACTAGGACGAATGAGCACGCAAAAACTGCAACCCAtt < 2:487461/150‑1 (MQ=255) gTTTTGCCTGGGACGCGCCCGGCCTGGACCGAGGTTTGCGGGTTGTACATGCGAATGCTCCTTTTCGTGGACGAACCCGAAATGTGGGTGGCTCACTTACTAGGACGAATGAGCACGCAAAAACTGCAACCCATTCGTGAAAAAAAATGc < 1:559406/150‑1 (MQ=255) gACCGAGGTTTGCGGGTTGTACATGCGAATGCTCCTTTTCGTGGACGAACCCGGAATGTGGGTGGCTCACTTACTAGGACGAATGAGCa > 2:299359/1‑89 (MQ=255) gACCGAGGTTTGCGGGTTGTACATGCGAATGCTCCTTTTCGTGGACGAACCCGGAATGTGGGTGGCTCACTTACTAGGACGAATGAGCa < 1:299359/89‑1 (MQ=255) aggggTTTGCGGGTTTTACAGGCGAATGCGCCGTTTCGTGGAGGAACCCGAAATGTGGTTGGCTCACTTACTAGGACCAATGAGCACGCAAAAACTGCAACCAAAGCATAAAAAAAAATGCCTGTTTCTCATCAAGACGTTCCATCGCCt < 2:166206/147‑1 (MQ=255) gAGGTTTGCGGGTTGTACATGCGAATGCTCCTTTTCGTGGACGAACCCGAAATTTGGGTGGCCCACTTACTAGGACGAATGAGCACGCAAAAACTGCAACCCATTCGTGAAAAAAAATGCCTGTTTCTCATCAAGACGTTCCATCGCCtt > 2:200216/1‑150 (MQ=255) | GGTGCGCAGGCCGAAGGTGCGCAACAGCTCGCGGATCTGTTCGTTGCCGATTCGTTCTTCAGCAAACACGTCCTGTCCCTTGGACGTTTTGCCTGGGACGCGCCCGGCCTGGACCGAGGTTTGCGGGTTGTACATGCGAATGCTCCTTTTCGTGGACGAACCCGAAATGTGGGTGGCTCACTTACTAGGACGAATGAGCACGCAAAAACTGCAACCCATTCGTGAAAAAAAATGCCTGTTTCTCATCAAGACGTTCCATCGCCTT > NC_002947/3297501‑3297765 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |