Skip to content

Problem preping stats from bam files ONCOCNV_getCounts.pl line 450 #5

@kmavrommatis

Description

@kmavrommatis

Hi,
I am using OncoCNV 6.6 to get copy number estimates on exome data.
There is an error that comes up when running the ONCOCNV_getCounts.pl

$perl $TOOLDIR/ONCOCNV_getCounts.pl getControlStats -m Exon -b $targetBed -c $controls -o $OUTDIR/Control.stats.txt

reading NA12891_Capture.MarkDuplicates.mdup..GATK.Recalibrate.bam
	sample name: NA12891_Capture.MarkDuplicates.mdup..GATK.Recalibrate

[ many lines counting reads ]
read 95700000 reads
	Total target length: 0
processed 2 controls, NA12878_Capture.MarkDuplicates.mdup..GATK.Recalibrate NA12891_Capture.MarkDuplicates.mdup..GATK.Recalibrate
Can't use an undefined value as an ARRAY reference at
	/opt/ONCOCNV-6.6/src//ONCOCNV_getCounts.pl line 450 (#1)
    (F) A value used as either a hard reference or a symbolic reference must
    be a defined value.  This helps to delurk some insidious errors.
    
Uncaught exception from user code:
	Can't use an undefined value as an ARRAY reference at /opt/ONCOCNV-6.6/src//ONCOCNV_getCounts.pl line 450.

the input files are :
bed file with the exome target regions

$head 1-target.bed 
chr1	12100	12258	AMPLENSG00000223972.4	0	ENSG00000223972.4
chr1	12555	12721	AMPLENSG00000223972.4.1	0	ENSG00000223972.4
chr1	13333	13701	AMPLENSG00000223972.4.2	0	ENSG00000223972.4
chr1	30336	30503	AMPLENSG00000243485.2	0	ENSG00000243485.2
chr1	35047	35544	AMPLENSG00000237613.2	0	ENSG00000237613.2
chr1	35620	35778	AMPLENSG00000237613.2.1	0	ENSG00000237613.2
chr1	69091	70008	AMPLENSG00000186092.4	0	ENSG00000186092.4
chr1	324296	324394	AMPLENSG00000237094.6	0	ENSG00000237094.6
chr1	324429	325605	AMPLENSG00000237094.6.1	0	ENSG00000237094.6
chr1	327736	328214	AMPLENSG00000237094.6.2	0	ENSG00000237094.6

and a bam file processed with Picard MarkDuplicates, GATK Realign and Recalibrate.
Before any processing since the BAM file contains aligned reads to unplaced scaffolds which are not present in the target.bed I am filtering with samtools -L <target.bed>

Any advice/help/suggestions are welcome
Thanks
Konstantinos

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions