sh 从CDS #fasta中提取fasta序列

Posted

tags:

篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了sh 从CDS #fasta中提取fasta序列相关的知识，希望对你有一定的参考价值。

exon_file="N2_GFF_WS245/CDS.gff.gz"
reference="reference/c_elegans.PRJNA13758.WS245.genomic.fa"
paste <(gunzip -kfc ${exon_file})  <(bedtools getfasta -fi ${reference} -bed ${exon_file} -tab -fo -) | grep 'WormBase'

以上是关于sh 从CDS #fasta中提取fasta序列的主要内容，如果未能解决你的问题，请参考以下文章

如何从大fasta文件中找出自己想要的序列

Linux根据基因ID提取fasta序列

请教如何从FASTA文件中批量查找序列

python文本处理---fasta文件提取指定ID的序列

fasta文件拆分与合并

perl 多fasta文件匹配，并提取匹配文件第一条序列