This is a test version of Biostars. For the public version, visit https://www.biostars.org.
GetPileupSummaries intervals-list with Targeted Sequencing?

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
• 2,032 views • link
updated by Ram Ram • written by Manuel Sokolov Ravasqueira • Manuel Sokolov Ravasqueira

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

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

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?)

• 0 views • link
written by Manuel Sokolov Ravasqueira • Manuel Sokolov Ravasqueira

kit is proprietary

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

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

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.

• 0 views • link
written by Manuel Sokolov Ravasqueira • Manuel Sokolov Ravasqueira

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

???????

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

0 answers

No answers yet.

Log in to answer this question.

More posts like this