生信技能34 - 通过UCSC单条染色体Fasta序列构建hg19参考序列及其索引

生信与基因组学

于 2024-01-08 12:11:17 发布

阅读量638

点赞数 10

分类专栏：生信分析项目实战技能合集文章标签： bash 数据分析

本文链接：https://blog.csdn.net/LittleComputerRobot/article/details/135453724

版权

生信分析项目实战技能合集专栏收录该内容

75 篇文章 ¥29.90 ¥99.00

订阅专栏

超级会员免费看

这篇博客介绍了如何通过UCSC下载染色体Fasta序列，然后使用bash脚本合并成hg19参考序列，并利用bwa构建其索引，涉及的数据分析和bash技术。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

1. UCSC文件下载

# 染色体长度文件
wget http://hgdownload.soe.ucsc.edu/goldenPath/hg19/bigZips/hg19.chrom.sizes
# md5校验文件
wget http://hgdownload.soe.ucsc.edu/goldenPath/hg19/bigZips/md5sum.txt
wget http://hgdownload.soe.ucsc.edu/goldenPath/hg19/bigZips/chromFa.tar.gz

# 检查chromFa.tar.gz完整性
cat md5sum.txt |grep chromFa.tar.gz > chromFa.tar.gz.md5.check
md5sum