Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

about somatic SV calling #364

Open
dongjuleem opened this issue Jan 22, 2024 · 0 comments
Open

about somatic SV calling #364

dongjuleem opened this issue Jan 22, 2024 · 0 comments

Comments

@dongjuleem
Copy link

dongjuleem commented Jan 22, 2024

Sorry for the basic question. I am currently conducting somatic SV calling using WGS paired data from 30 tumors and controls. I have used the pipeline to call each paired data as follows:

delly call -x hg19.excl -o t1.bcf -g hg19.fa tumor1.bam control1.bam
delly filter -f somatic -o t1.pre.bcf -s samples.tsv t1.bcf

As a first question, I noticed that some candidates with filter=pass in delly call disappear after the delly filter process.
I would like to inquire about the reason for this and whether this process is correct.

As a final question, in the README, it recommends the following step:
delly call -g hg19.fa -v t1.pre.bcf -o geno.bcf -x hg19.excl tumor1.bam control1.bam ... controlN.bam

Does this mean I have to repeat this process 30 times for each of the 30 WGS paired data?"

always thanks

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant