HLAforest:使用RNA-seq数据进行HLA分型

HLAforest:使用RNA-seq数据进行HLA分型

欢迎关注"生信修炼手册"！

HLAforest是一款针对RNA_seq数据进行HLA分型的软件，github地址如下

​​clone 修改config.sh

主要是修改​​HLAFOREST_HOME​​​和​​NUM_THREADS​​这两个选项，示例如下

#!/bin/bash# User modifiable variablesHLAFOREST_HOME=/soft/hlaforest/NUM_THREADS=10 # the number of threads bowtie should use

​​HLAFOREST_HOME​​​指定软件的安装目录，​​NUM_THREADS​​指定bowtie比对的线程数。

2. 修改 CallSimulation.sh

调整​​CONFIG_PATH​​的设置

CONFIG_PATH=/soft/hlaforest/scripts/config.sh

3. 修改 CallHaplotypesPE.sh

调整​​CONFIG_PATH​​的设置

CONFIG_PATH=/soft/hlaforest/scripts/config.sh

4. 添加PATH环境变量

在环境变量的配置文件中,将​​scripts​​目录添加到PATH变量中，方便程序调用，写法如下

export PATH=/soft/hlaforest/scripts:$PATH

上述参数都调整好之后，软件就可以运行了。需要注意的是该软件依赖bioperl和bowtie,由于我的系统是已经装过了的，所以这里没给出这两个软件的安装过程。

软件自带测试数据集，用法如下

CallHaplotypesPE.sh test2/ test2/gm12878_short_1.fastq test2/gm12878_short_2.fastq

只需要指定双端测序的原始文件就可以了，默认分型结果保存在​​haplotypes.txt​​中，该文件内容示意如下

ROOT:A ROOT:A:11 ROOT:A:11:01 ROOT:A:11:01:04ROOT:A ROOT:A:01 ROOT:A:01:01 ROOT:A:01:01:01 ROOT:A:01:01:01:01ROOT:B ROOT:B:55 ROOT:B:55:01 ROOT:B:55:01:01ROOT:B ROOT:B:08 ROOT:B:08:01 ROOT:B:08:01:01

在上述文件中，HLA-A基因的分型结果为​​A:11:01:04​​​和​​A:01:01:01:01​​。

更多参数和用法请参考软件的帮助文档。

版权声明：本文内容由网络用户投稿，版权归原作者所有，本站不拥有其著作权，亦不承担相应法律责任。如果您发现本站中有涉嫌抄袭或描述失实的内容，请联系我们jiasou666@gmail.com 处理，核实后本网站将在24小时内删除侵权内容。