2025-11-24從fasta和gff文件中提取gDNA

1.從gff2bed

awk'$3=="gene" {print $1 "\t" $4-1 "\t" $5 "\t" $9 "\t0\t" $7}' hap1.gff3 > hap1.genes.bed

2.基于bed文件，利用bedtools提取基因gDNA序列

bedtools getfasta -fi hap1.fa -bed hap1.genes.1000.1.bed -fo hap1.genes.gdna.fasta -name -s

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成，瀏覽時請結(jié)合常識與多方信息審慎甄別。
平臺聲明：文章內(nèi)容（如有圖片或視頻亦包括在內(nèi)）由作者上傳并發(fā)布，文章內(nèi)容僅代表作者本人觀點(diǎn)，簡書系信息發(fā)布平臺，僅提供信息存儲服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

友情鏈接更多精彩內(nèi)容