Issues with BedToIntervalList Picard tool
1
0
Entering edit mode
19 months ago
Elisa • 0

Hi everyone,

I'm trying to convert a bed file into interval list format. I'm not sure about which file I should choose as SD. I've tried using 'genome.dict' (genome.fa is the reference genome I've used to align) or the BAM file of the aligned reads; in both cases I get the following error:

Exception in thread "main" picard.PicardException: Sequence 'chr1_KI270763v1_alt' was not found in the sequence dictionary

How can I solve this issue?

Thanks in advance

Picard • 741 views
ADD COMMENT
0
Entering edit mode
19 months ago

Sequence 'chr1_KI270763v1_alt' was not found in the sequence dictionary

your bed file contains a chromosome "chr1_KI270763v1_alt" that is not defined in the REFERENCE/DICT file. if you're not interested with that contig, you can just remove it from your bed with grep -v

ADD COMMENT
0
Entering edit mode

I'm interested with that contig so I can't remove it. Is it allowed to add that contig to the reference genome?

ADD REPLY
0
Entering edit mode

it depends of your downstream analysis. For example if you're using a bam and gatk, the internal dictionary should be the very same than the reference sequence.

ADD REPLY

Login before adding your answer.

Traffic: 1281 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