aln2snps
msa_position | p | q | ref_position | ref_seq | type | gene | start | end | protein_id | feature_position | feature_position_mod3 | triplet_start | triplet_end | ref_triplet | mut_triplet | ref_AA | mut_AA | ref_freq | mut_freq | mut_seq |
3653
|
36626
|
1.98774342659889E-68
|
4.81792614569765E-68
|
23402
|
A
|
CDS
|
S
|
21563
|
25384
|
YP_009724390.1
|
1841
|
2
|
23401
|
23403
|
GAT
|
GGT
|
D
|
G
|
11231
|
30413
|
G
|
D614 to G614 change ~ Dewpoint2meter (humidity)
614x3=1842, so, position 1841 is the 2nd nucleotide position. This is consistent with mod3=2.
Code:
https://github.com/QinLab/aln2snps/blob/master/annotate_gwas_with_GFF.ipynb
Output:
https://github.com/QinLab/aln2snps/blob/master/0613/gwas-dewpoint-meta0613-positions-pq-gff-draft20200713.csv
There are other spike mutations. Are they mentioned in previous publications?
No comments:
Post a Comment