Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • krapulaxdoctor
    Member
    • May 2015
    • 22

    #1

    picard metrics error

    Dear all,
    I would like to ask for help to solve this issue with Picard tools:

    I already looked through forums, picard FAQ, but I id not find similar problems.

    I would like to run the CollectOxoGMetrics command to test my datafiles.
    For this i used the following command:


    Code:
    user@work:~$ java -jar -Xmx6g /picard/dist/picard.jar CollectOxoGMetrics R=/mm10-Genome/genome.fa I=rnaseq.bam O=picard-oxo-output.txt
    the program starts up and writes some line:

    [Fri Oct 23 13:15:24 CEST 2015] picard.analysis.CollectOxoGMetrics INPUT=rnaseq.bam OUTPUT=picard-oxo-output.txt REFERENCE_SEQUENCE=/mm10-Genome/genome.fa MINIMUM_QUALITY_SCORE=20 MINIMUM_MAPPING_QUALITY=30 MINIMUM_INSERT_SIZE=60 MAXIMUM_INSERT_SIZE=600 USE_OQ=true CONTEXT_SIZE=1 STOP_AFTER=2147483647 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
    [Fri Oct 23 13:15:24 CEST 2015] Executing as user@work on Linux 3.13.0-65-generic amd64; OpenJDK 64-Bit Server VM 1.7.0_79-b14; Picard version: 1.135(a61ab80f95b753621cb252b17f2234d8bcb77096_1436675286) JdkDeflater
    INFO 2015-10-23 13:15:24 CollectOxoGMetrics Generated 16 context strings.
    INFO 2015-10-23 13:15:24 CollectOxoGMetrics Loading dbSNP File: null
    INFO 2015-10-23 13:15:24 CollectOxoGMetrics Starting iteration.

    , but after 2 minutes it stops, writing:

    [Fri Oct 23 13:17:57 CEST 2015] picard.analysis.CollectOxoGMetrics done. Elapsed time: 2.54 minutes.
    Runtime.totalMemory()=229638144
    And the text file it generates as output contains a few line:

    ## htsjdk.samtools.metrics.StringHeader
    # picard.analysis.CollectOxoGMetrics INPUT=rnaseq.bam OUTPUT=picard-oxo-output.txt REFERENCE_SEQUENCE=/mm10-Genome/genome.fa MINIMUM_QUALITY_SCORE=20 MINIMUM_MAPPING_QUALITY=30 MINIMUM_INSERT_SIZE=60 MAXIMUM_INSERT_SIZE=600 USE_OQ=true CONTEXT_SIZE=1 STOP_AFTER=2147483647 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
    ## htsjdk.samtools.metrics.StringHeader
    # Started on: Fri Oct 23 13:15:24 CEST 2015
    And nothing else.

    Then , i wanted to test whether my bam file has a problem, for this i ran

    Code:
    java -jar -Xmx6g /picard/dist/picard.jar ValidateSamFile I=rnaseq.bam O=picard-validate-output.txt
    the output files has error messages, (100 lines from which only one is different) such as:

    ERROR: Read groups is empty
    WARNING: Read name BIOMICS-HISEQLN:157:HMLNVADXX:2:1205:5080:8933, A record is missing a read group
    WARNING: Read name BIOMICS-HISEQHI:555:HMLGGADXX:2:1211:16446:43292, A record is missing a read group
    WARNING: Read name BIOMICS-HISEQLN:159:HMJFYADXX:1:1205:18060:52304, A record is missing a read group
    WARNING: Read name BIOMICS-HISEQHI:555:HMLGGADXX:2:2210:16735:34742, A record is missing a read group
    WARNING: Read name BIOMICS-HISEQLN:157:HMLNVADXX:2:1111:19948:49163, A record is missing a read group
    WARNING: Read name BIOMICS-HISEQHI:555:HMLGGADXX:2:1111:9860:29396, A record is missing a read group
    WARNING: Read name BIOMICS-HISEQLN:157:HMLNVADXX:1:1103:12058:44655, A record is missing a read group
    WARNING: Read name BIOMICS-HISEQLN:159:HMJFYADXX:1:2115:4058:6432, A record is missing a read group
    WARNING: Read name BIOMICS-HISEQLN:157:HMLNVADXX:1:2208:5115:74838, A record is missing a read group

    Could you please suggest me an idea how to deal with this error?
    Does picard need more than 8GB memory in general?

    Thank you very much for your help. And for your time.


    Regards,

Latest Articles

Collapse

  • SEQadmin2
    Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
    by SEQadmin2



    Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

    Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
    ...
    10-14-2026, 09:55 AM
  • SEQadmin2
    New Genomics Technologies Take Aim at Long-Standing Limits
    by SEQadmin2


    Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

    We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
    ...
    09-28-2026, 10:25 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, 10-22-2026, 11:13 AM
0 responses
29 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-22-2026, 10:48 AM
0 responses
22 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-15-2026, 10:52 AM
0 responses
39 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-14-2026, 05:39 AM
0 responses
45 views
0 reactions
Last Post SEQadmin2  
Working...