vcf文件新建索引(idx)

这篇博客记录了如何利用GATK的IndexFeatureFile参数来为Mills_and_1000G_gold_standard.indels.hg38.vcf文件创建索引。执行过程中详细展示了使用的软件版本、Java运行环境以及操作系统的相关信息,并在完成后成功写入了索引文件。
摘要由CSDN通过智能技术生成

原文件如下:

 直接使用gatk 的 IndexFeatureFile 参数:

#gatk IndexFeatureFile -F  vcf_file 
$ gatk IndexFeatureFile -F Mills_and_1000G_gold_standard.indels.hg38.vcf
Using GATK jar /share_data/wujm/software/gatk-4.1.4.0/gatk-package-4.1.4.0-local.jar
Running:
    java -Dsamjdk.use_async_io_read_samtools=false -Dsamjdk.use_async_io_write_samtools=true -Dsamjdk.use_async_io_write_tribble=false -Dsamjdk.compression_level=2 -jar /share_data/wujm/software/gatk-4.1.4.0/gatk-package-4.1.4.0-local.jar IndexFeatureFile -F Mills_and_1000G_gold_standard.indels.hg38.vcf
13:23:38.650 INFO  NativeLibraryLoader - Loading libgkl_compression.so from jar:file:/share_data/wujm/software/gatk-4.1.4.0/gatk-package-4.1.4.0-local.jar!/com/intel/gkl/native/libgkl_compression.so
Dec 06, 2021 1:23:38 PM shaded.cloud_nio.com.google.auth.oauth2.ComputeEngineCredentials runningOnComputeEngine
INFO: Failed to detect whether we are running on Google Compute Engine.
13:23:38.875 INFO  IndexFeatureFile - ------------------------------------------------------------
13:23:38.875 INFO  IndexFeatureFile - The Genome Analysis Toolkit (GATK) v4.1.4.0
13:23:38.875 INFO  IndexFeatureFile - For support and documentation go to https://software.broadinstitute.org/gatk/
13:23:38.875 INFO  IndexFeatureFile - Executing as wujianming@pang on Linux v3.10.0-957.el7.x86_64 amd64
13:23:38.876 INFO  IndexFeatureFile - Java runtime: OpenJDK 64-Bit Server VM v1.8.0_152-release-1056-b12
13:23:38.876 INFO  IndexFeatureFile - Start Date/Time: 2021年12月6日 下午01时23分38秒
13:23:38.876 INFO  IndexFeatureFile - ------------------------------------------------------------
13:23:38.876 INFO  IndexFeatureFile - ------------------------------------------------------------
13:23:38.876 INFO  IndexFeatureFile - HTSJDK Version: 2.20.3
13:23:38.876 INFO  IndexFeatureFile - Picard Version: 2.21.1
13:23:38.876 INFO  IndexFeatureFile - HTSJDK Defaults.COMPRESSION_LEVEL : 2
13:23:38.876 INFO  IndexFeatureFile - HTSJDK Defaults.USE_ASYNC_IO_READ_FOR_SAMTOOLS : false
13:23:38.876 INFO  IndexFeatureFile - HTSJDK Defaults.USE_ASYNC_IO_WRITE_FOR_SAMTOOLS : true
13:23:38.876 INFO  IndexFeatureFile - HTSJDK Defaults.USE_ASYNC_IO_WRITE_FOR_TRIBBLE : false
13:23:38.876 INFO  IndexFeatureFile - Deflater: IntelDeflater
13:23:38.877 INFO  IndexFeatureFile - Inflater: IntelInflater
13:23:38.877 INFO  IndexFeatureFile - GCS max retries/reopens: 20
13:23:38.877 INFO  IndexFeatureFile - Requester pays: disabled
13:23:38.877 INFO  IndexFeatureFile - Initializing engine
13:23:38.877 INFO  IndexFeatureFile - Done initializing engine
13:23:39.092 INFO  FeatureManager - Using codec VCFCodec to read file file:///share_data/wujm/Config/reference/hg38/Mills_and_1000G_gold_standard.indels.hg38.vcf
13:23:39.138 INFO  ProgressMeter - Starting traversal
13:23:39.138 INFO  ProgressMeter -        Current Locus  Elapsed Minutes     Records Processed   Records/Minute
13:23:41.420 INFO  ProgressMeter -       chrX:152403569              0.0               1271815       33454142.9
13:23:41.421 INFO  ProgressMeter - Traversal complete. Processed 1271815 total records in 0.0 minutes.
13:23:41.594 INFO  IndexFeatureFile - Successfully wrote index to /share_data/wujm/Config/reference/hg38/Mills_and_1000G_gold_standard.indels.hg38.vcf.idx
13:23:41.594 INFO  IndexFeatureFile - Shutting down engine
[2021年12月6日 下午01时23分41秒] org.broadinstitute.hellbender.tools.IndexFeatureFile done. Elapsed time: 0.05 minutes.
Runtime.totalMemory()=2573729792
Tool returned:
/share_data/wujm/Config/reference/hg38/Mills_and_1000G_gold_standard.indels.hg38.vcf.idx

 ----------   如果这篇文章对您有帮助,请您随手点个赞,谢谢啦~~~~     ---------

评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

heianduck

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值