vcf文件新建索引(idx)

原文件如下:

vcf文件新建索引(idx)

 直接使用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

vcf文件新建索引(idx)

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

上一篇:GFPGAN源码分析—第四篇


下一篇:【路径规划】基于RRT算法机器人最短路径规划matlab代码