1000geno_maf
2024-04-03 16:20:00
千人基因组计划相关SNP下载
在网址https://ftp-trace.ncbi.nih.gov/1000genomes/ftp/release/20100804/ALL.2of4intersection.20100804.genotypes.vcf.gz获得千人基因组计划数据
这是一个大小为64G的文件,为了加快下载进程,使用ascp下载命令进行下载
首先进入aspera官方网站,找到「IBM Aspera Connect」,进入下载页面,找到对应的版本与平台,这里以Linux平台最新的3.10.0版本为例进行介绍。
# 下载
wget https://ak-delivery04-mul.dhe.ibm.com/sar/CMA/OSA/092u0/0/ibm-aspera-connect-3.10.0.180973-linux-g2.12-64.tar.gz
# 解压
tar xvf ibm-aspera-connect-3.10.0.180973-linux-g2.12-64.tar.gz
# 解压后得到一个脚本文件,运行该脚本,即可完成自动安装
sh ibm-aspera-connect-3.10.0.180973-linux-g2.12-64.sh
# 所有安装文件都在~/.aspera/connect目录下,添加环境变量
echo 'export PATH=~/.aspera/connect/bin/:$PATH' >> ~/.bashrc
# 使环境变量生效
source ~/.bashrc
# 查看ascp可执行文件所在的路径,应该是:~/.aspera/connect/bin/ascp
which ascp
# 看程序是否能正常运行
ascp -h
#人类基因组计划相关SNP下载
ascp -i ~/.aspera/connect/etc/asperaweb_id_dsa.openssh -l 100M -T -P33001 fasp-g1k@fasp.1000genomes.ebi.ac.uk:vol1/ftp/release/20100804/ALL.2of4intersection.20100804.genotypes.vcf.gz .
在千人基因组文件所在目录下打开终端后,输入plink
#下一行代码将输出bed,bim,fam三个文件,其中bed是一个二进制文件
plink --vcf ALL.2of4intersection.20100804.genotypes.vcf.gz --make-bed --out ALL.2of4intersection.20100804.genotypes
#下面这条命令中的sample.txt由两列构成,第一列为提取的样本Family ID,第二列为Within-family ID(IID)
plink --bfile ALL.2of4intersection.20100804.genotypes --keep sample.txt --recode --make-bed --out ch_sample
#计算等位基因频率
plink --bfile ch_sample --freq --out ch_MAF_check
1000geno_maf相关推荐
最新文章
- 如何指定在Git上执行shell命令时使用的私有SSH密钥?
- 机房漏水后,漏水报警器是如何检测应急处理的?
- 【AJAX】AJAX实现搜索信息自己主动推荐并补全
- 超详细JDK下载与安装步骤
- PP苹果助手 v2.3.0 官方版
- python应用学习(五)——requests爬取网页图片
- CCKS2020基于本体的金融知识图谱自动化构建技术评测第五名方法总结
- 0708-/etc/login.defs
- visual studio 实时调试器问题
- js 浏览器窗口活跃监听
- 利用Python+xarray+cartopy+matplotlib 实现遥感地形图制图绘制 —— xarray 学习文档01
- hashcat破解WiFi显示No hashes loaded的解决方法
- SLCP验厂辅导,工厂进行社会劳工整合项目(SLCP认证)验证的具体步骤
- 2021春招面经系列--计算机网络和操作系统
- java中构造方法(或者叫构造方法)
- VS2019加QT5.14中Please assign a Qt installation in ‘Qt Project Settings‘.问题的解决
- MC30P6280晟矽微原厂是一款性能好低功耗8位OTP+IO型MCU
- Python上位机,监视串口,逐字符匹配字头
- 关于卸载公司强制安装的杀毒软件的详细步骤,亲测有效。
- 【伺服电机/舵机与Arduino使用教程】