运行和验证
操作步骤
- 使用PuTTY工具,以root用户登录服务器。
- 下载测试算例包“POST_GATK_DATA.zip”。
- 使用SFTP工具将测试算例包上传至服务器“/path/to/TEST”路径下。
- 执行以下命令解压算例包。
cd /path/to/TEST unzip POST_GATK_DATA.zip
- 执行以下命令,用XHMM合并GATK覆盖深度输出的多个样本(在同一位点)。
xhmm --mergeGATKdepths -o ./DATA.RD.txt \ --GATKdepths group1.DATA.sample_interval_summary \ --GATKdepths group2.DATA.sample_interval_summary \ --GATKdepths group3.DATA.sample_interval_summary
运行结果如下:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Command-line parameter values: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% mergeGATKdepths GATKdepths="group1.DATA.sample_interval_summary" GATKdepths="group2.DATA.sample_interval_summary" GATKdepths="group3.DATA.sample_interval_summary" outputMatrix="./DATA.RD.txt" %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Reading GATK depth file group1.DATA.sample_interval_summary Reading GATK depth file group2.DATA.sample_interval_summary Reading GATK depth file group3.DATA.sample_interval_summary Writing GATK read-depth matrix of 30 samples by 300 targets to ./DATA.RD.txt Total CPU time for processing this job: 6 seconds [user time: 6 system time: 0 children user time: 0 children system time: 0] Total time for processing this job: 0 seconds
生成文件:
DATA.RD.txt