超级会员免费看
牛转录组分析的实用生物信息学流程
1. 数据提取与基因组索引生成
1.1 特征提取
可以使用以下命令从 Angus 的 gfff/gtf 文件中提取特征：
zgrep -P ‘^29\t’ Bos_indicus_hybrid.UOA_Brahman_1.104.gtf.gz > UOA_Brahman_1 _chr29/Brahman_chr29.gtf
或者使用以下替代代码：
unzip Bos_indicus_hybrid.UOA_Brahman_1.104.gtf.gz
grep -P ‘^29\t’ Bos_indicus_hybrid.UOA_Brahman_1.104.gtf > UOA_Brahman_1 _chr29/Brahman_chr29.gtf
1.2 生成基因组索引
在进行 RNA 读取比对之前，需要对参考序列进行索引。使用以下命令：
hisat2-build Bos_indicus_hybrid.UOA_Brahman_1.cds.all.fa.gz UOA_Brahman_1 _chr29/ Brahman_chr29.ref
2. 读取比对
接下来，使用 HISAT2 比对器进行读取比对，并为 Brahman 和 Angus 牛的每个生物重复生成相应的比对文
订阅专栏 解锁全文
扫一扫
225
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
