2025-08-26提取最长转录本2

1.先提取最长转录本gff（conda install -c bioconda agat -y）

perl ~/miniconda3/bin/agat_sp_keep_longest_isoform.pl --gff A11.gff -o longest_isoforms.gff

2.基于最长转录本gff提取转录本序列

gffread longest_isoforms.gff -g A11_chr.fa -w longest_isoforms.fasta

最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。