This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Identification of break points for a balanced translocation

Hi,

From the karyotype results we could see a balanced translocation in one of the sample and we have whole genome data available for the same sample. It would be really great if someone can give suggestion in this regard. We have extracted the clipped reads aligned to both the chromosome and part of the reads which are clipped have been realigned to reference genome. We would like to cluster the reads which aligned to specific region. Any tools are available to perform this?

wgs breakpoint analysis balanced translocation
• 2,334 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by nkausthu • nkausthu

2 answers

If you have paired-end whole genome data you can use Delly.

• 1 views • link
written by trausch trausch

Thank you and I ran Delly for the paired end data. But its taking too much of time to give the results. As I have the karyotype results can I specify the specific chromosomes and run Delly?

• 1 views • link
written by nkausthu • nkausthu

You can just extract the 2 translocated chromosomes of interest:

samtools view -b <input.bam> chrA chrB > subset.bam
samtools index subset.bam

Or if you know the approximate regions:

samtools view -b <input.bam> chrA:10000000-20000000 chrB:50000000-60000000 > subset.bam
samtools index subset.bam
• 1 views • link
written by trausch trausch

last year, I found a balanced translocation just by looking at the discordant pairs Identifying discordantly mapped reads .

There was a large number of pairs that were obviously clustered in the same region

• 1 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this