GetPileupSummaries intervals-list with Targeted Sequencing?
0
1
Entering edit mode
10 weeks ago

Hi!

I am applying the GetPileUpSummaries, for somatic variant calling starting from targeted sequencing .fasta. I aligned the file with the GrCh38 reference. And currently I am at the GetPileUpSummariesStep.

  gatk --java-options -Xmx200G GetPileupSummaries \
    -I $RECBAM \
    -L ???? \
    -O $OUTPUT \
    -V $GERM # Germline reference --> small_exac_common_3.hg38.vcf.gz

For the -L argument i need an intervals.bed file. Given the fact that my fasta is targeted sequencing, which file should I use here?

Best Regards, Manuel

GetPileUpSummaries GATK Variant-Calling • 587 views
ADD COMMENT
0
Entering edit mode

how about using the bed file of your targets....

ADD REPLY
0
Entering edit mode

thanks for the suggestion. how can i obtain the bedfile of my targets (kit is proprietary - can I obtain bed from fasta? if yes how?)

ADD REPLY
0
Entering edit mode

kit is proprietary

I can't imagine you get a targeted sequencing without having a BED file....

ADD REPLY
0
Entering edit mode

they only provided me the sorted.bam files aligned with hg19. I converted bam into fastq and started variant calling from there. Asking for bed file is not a suitable option in the current scenario.

ADD REPLY
0
Entering edit mode

Asking for bed file is not a suitable option in the current scenario.

???????

ADD REPLY
1
Entering edit mode

SOLVED by using bedtools -> bamToBed tool.

bedtools bamtobed -i  $SORTED_BAM > intervals.bed
ADD REPLY

Login before adding your answer.

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