标签：fa NBS 模型 pep qua hmm ARC 构建

下载植物gff、cds、dna、pep
pfam中下载hmm模型
搜索基因家族并以1e-20筛选

hmmsearch --cut_tc --domtblout NBS-ABC.out NBS-ARC.hmm Arabidopsis_thaliana.TAIR10.pep.all.fa
grep -v "#" NBS-ABC.out|awk '($7 + 0) < 1E-20'|cut -f1 -d  " "|sort -u > NBS-ARC_qua_id.txt
less Arabidopsis_thaliana.TAIR10.pep.all.fa | /home/jun_qzt/miniconda3/envs/rna-seq/bin/seqkit grep -f NBS-ARC_qua_id.txt > NBS-ARC_qua.fa

多序列比对，构建模型

/home/jun_qzt/miniconda3/envs/rna-seq/bin/clustalw
弹出clustalw的操作界面，以下展示具体输入过程：
选择1（输入待比对序列）→ 输入待比对序列的文件名：NBS-ARC_qua.fa → 选择2（开始进行序列比对）→选择9（选择输出比对结构的格式为aligned）→ 按enter键 → 选择1（选择比对模式为全局比对） → 指定输出的比对结果的文件名称:NBS-ARC_qua.aln → 回车后开始比对 → 输入一个树文件名(new GUIDE TREE file):NBS-ARC_qua.dnd （最后才能得到NBS-ARC.aln，否则NBS-ARC.aln为空）

hmmbuild NBS-ARC_qua.hmm  NBS-ARC_qua.aln
hmmsearch --cut_tc --domtblout NBS-ARC.second.out NBS-ARC_qua.hmm Arabidopsis_thaliana.TAIR10.pep.all.fa

再次筛选符合植物的序列

grep -v "#" NBS-ARC.second.out|awk '($7 + 0) < 1E-03' | cut -f1 -d " "|sort -u >final.NBS.list
less Arabidopsis_thaliana.TAIR10.pep.all.fa.gz | /data1/spider/ytbiosoft/seqkit grep -f final.NBS.list > final_NBS-ARC_qua.fa

标签：fa,NBS,模型,pep,qua,hmm,ARC,构建
来源： https://www.cnblogs.com/qinziting/p/15575094.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。