Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,623,245 | C→T | 21.4% | intergenic (+38/‑208) | PP_5517 → / → PP_2295 | hypothetical protein/antirestriction protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,623,245 | 0 | C | T | 21.4% | 43.9 / 8.3 | 28 | intergenic (+38/‑208) | PP_5517/PP_2295 | hypothetical protein/antirestriction protein |
Reads supporting (aligned to +/- strand): ref base C (11/11); new base T (4/2); total (15/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.55e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01 |
TTCGGGGATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGT > NC_002947/2623099‑2623382 | ttCGGGGATGCGTGTGGAGGTGTTCCACCAGAATGGGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGCGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATCCGCAATCAGACATCAAGGGGCCCAATCggg > 2:98834/1‑150 (MQ=255) tCGGGGATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCGATCGGGc < 2:189196/150‑1 (MQ=255) gggATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCGATTGGGcccc < 2:174297/150‑2 (MQ=255) tgtgGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCGATTGGGcccct < 2:94633/143‑3 (MQ=255) tgtgGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCGATTGGGcccct > 1:94633/1‑141 (MQ=255) gAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTc > 2:9708/1‑143 (MQ=255) gAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTc < 1:9708/143‑1 (MQ=255) tACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTc > 2:404522/1‑150 (MQ=255) gTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAggg < 1:98834/150‑1 (MQ=255) ttccttTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAg > 1:590662/1‑150 (MQ=255) tcCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAg > 2:307236/1‑150 (MQ=255) ggTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCACTCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATc < 2:551556/149‑1 (MQ=255) gCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGaa < 1:404522/150‑1 (MQ=255) cTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAgc < 2:590662/150‑1 (MQ=255) ctacaccAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCaaaa > 1:705841/1‑150 (MQ=255) ctacaccAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCaaaa > 1:467178/1‑150 (MQ=255) cAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGa > 1:190417/1‑150 (MQ=255) gcgcTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGACCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTcc > 1:800740/1‑150 (MQ=255) gACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACaaa < 2:467178/150‑1 (MQ=255) cGCAATCACACCTCAAGGGGCGCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAGGACACGTCATTTGCGCAAGTGGAAATAAATCa < 2:438733/150‑1 (MQ=255) gCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCAc > 1:801055/1‑150 (MQ=255) cAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACt < 2:705841/150‑1 (MQ=255) ccTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATggg < 1:307236/150‑1 (MQ=255) ccTCAAGGGGCCCAATCGGGCCCTTTTTCGCGGCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCTCACATGTAAATAAATCACTTGCaggtg < 2:10599/150‑5 (MQ=255) cTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGCCAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAACGGAAATAAATCACTTGCATgggg > 1:716672/1‑150 (MQ=255) aaaagggCCCGATTGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGgtg > 1:484708/5‑150 (MQ=255) aaaagggCCCGATTGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAAATCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGgtg > 1:478243/5‑150 (MQ=255) aaagggCCCGATTGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGgtgt > 2:342664/4‑150 (MQ=255) | TTCGGGGATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGT > NC_002947/2623099‑2623382 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |