Intrerval list from bed file
Hello
I am trying to convert my bed file as intervals list. My reference dictionary and bed file are from same genomic build GRCh38. But still it gives the below error
java -jar $EBROOTPICARD/picard.jar BedToIntervalList \
I=/scicore/home/cichon/GROUP/memory_optimization/Twist_ComprehensiveExome_targets_hg38.bed \
O=/scicore/home/cichon/GROUP/memory_optimization/list.interval_list \
SD=/scicore/home/cichon/GROUP/memory_optimization/data/reference/gch38.dict
To execute picard run: java -jar $EBROOTPICARD/picard.jar
[Mon Nov 08 09:58:00 CET 2021] picard.util.BedToIntervalList INPUT=/scicore/home/cichon/GROUP/memory_optimization/Twist_ComprehensiveExome_targets_hg38.bed SEQUENCE_DICTIONARY=/scicore/home/cichon/GROUP/memory_optimization/data/reference/gch38.dict OUTPUT=/scicore/home/cichon/GROUP/memory_optimization/list.interval_list SORT=true UNIQUE=false 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
[Mon Nov 08 09:58:00 CET 2021] Executing as thirun0000@shi20.cluster.bc2.ch on Linux 3.10.0-1062.18.1.el7.x86_64 amd64; OpenJDK 64-Bit Server VM 1.8.0_212-b03; Picard version: 2.6.0-SNAPSHOT
[Mon Nov 08 09:58:01 CET 2021] picard.util.BedToIntervalList done. Elapsed time: 0.01 minutes.
Runtime.totalMemory()=32440320
To get help, see http://broadinstitute.github.io/picard/index.html#GettingHelp
Exception in thread "main" picard.PicardException: Sequence 'chr1' was not found in the sequence dictionary
at picard.util.BedToIntervalList.doWork(BedToIntervalList.java:146)
at picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:208)
at picard.cmdline.PicardCommandLine.instanceMain(PicardCommandLine.java:95)
at picard.cmdline.PicardCommandLine.main(PicardCommandLine.java:105)
Priya
-
Hi Priyadarshini Thirunavukkarasu,
It looks like this is most likely due to an incompatibility in the naming of the chromosomes between your bed file and dictionary file. Could you take a look at this document about incompatibilities in input files and check to make sure that the chromosome naming matches between your files?
Kind regards,
Pamela
Please sign in to leave a comment.
1 comment