Closed:Picard tools: error with CollectInsertSizeMetrics
0
0
Entering edit mode
6.2 years ago
anais1396 ▴ 30

Hi everyone !!

I try to make a quality control of my ATAC-seq datas with the Picard tool CollectInserSizeMetrics. As input, I provided a BAM files which is sorted and with removed duplicates (this file is the output of the Picard tool MarkDuplicates)

I ran the following command line (as we can see here : https://broadinstitute.github.io/picard/command-line-overview.html) :

java -jar /gpfs/work/amontois/Programmes/picard/build/libs/picard.jar CollectInsertSizeMetrics I=alignSRR75_marked_duplicates.bam O=alignSRR75_insert_size_metrics.txt H=alignSRR75_insert_size_histogram.pdf M=0.5

But no file.txt and no histogram were generated...

I had the following error message :

12:06:51.065 INFO  NativeLibraryLoader - Loading libgkl_compression.so from jar:file:/gpfs/work/amontois/Programmes/picard/build/libs/picard.jar!/com/intel/gkl/native/libgkl_compression.so
[Sun Feb 18 12:06:51 CET 2018] CollectInsertSizeMetrics HISTOGRAM_FILE=alignSRR75_insert_size_histogram.pdf INPUT=alignSRR75_marked_duplicates.bam OUTPUT=alignSRR75_insert_size_metrics.txt    DEVIATIONS=10.0 MINIMUM_PCT=0.05 METRIC_ACCUMULATION_LEVEL=[ALL_READS] INCLUDE_DUPLICATES=false ASSUME_SORTED=true STOP_AFTER=0 VERBOSITY=INFO QUIET=false VALIDATION_STRINGENCY=STRICT COMPRESSION_LEVEL=5 MAX_RECORDS_IN_RAM=500000 CREATE_INDEX=false CREATE_MD5_FILE=false GA4GH_CLIENT_SECRETS=client_secrets.json USE_JDK_DEFLATER=false USE_JDK_INFLATER=false
[Sun Feb 18 12:06:51 CET 2018] Executing as amontois@nic65 on Linux 3.10.0-693.11.6.el7.x86_64 amd64; Java HotSpot(TM) 64-Bit Server VM 1.8.0_152-b16; Deflater: Intel; Inflater: Intel; Picard version: 2.17.6-8-ga8222cd-SNAPSHOT
WARNING 2018-02-18 12:06:51     CollectInsertSizeMetrics        All data categories were discarded because they contained < 0.05 of the total aligned paired data.
WARNING 2018-02-18 12:06:51     CollectInsertSizeMetrics        Total mapped pairs in all categories: 0.0
[Sun Feb 18 12:06:51 CET 2018] picard.analysis.CollectInsertSizeMetrics done. Elapsed time: 0.00 minutes.
Runtime.totalMemory()=2075918336

I don't understand what is wrong... Can someone help me to fix that ??

Thank you in advance !

Anais

sequencing ATAC-seq • 211 views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 1968 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6