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

Delly vcf output #358

Open
Ana1335 opened this issue Dec 10, 2023 · 3 comments
Open

Delly vcf output #358

Ana1335 opened this issue Dec 10, 2023 · 3 comments
Assignees
Labels

Comments

@Ana1335
Copy link

Ana1335 commented Dec 10, 2023

Hello, I am running delay on three cancer cell lines to find TVs.
This is my running command (I found it is very time consuming to run all TVs together so running them separately)

for name in DEL DUP INV TRA INS; 
    do delly call -t ${name} -g ${FASTA} ${INPUT_DIR}/${SAMPLE_ID}.rmdup1.bam  > ${OUT_DIR}/${SAMPLE_ID}_${name}_delly.vcf
done

but the output vcfs I get have only these with no variants:

0%   10   20   30   40   50   60   70   80   90   100%
|----|----|----|----|----|----|----|----|----|----|
***************************************************
[2023-Dec-10 05:32:43] Split-read alignment

0%   10   20   30   40   50   60   70   80   90   100%
|----|----|----|----|----|----|----|----|----|----|
***************************************************
[2023-Dec-10 06:28:07] Junction read annotation

0%   10   20   30   40   50   60   70   80   90   100%
|----|----|----|----|----|----|----|----|----|----|
***************************************************
[2023-Dec-10 06:47:33] Breakpoint spanning coverage annotation

0%   10   20   30   40   50   60   70   80   90   100%
|----|----|----|----|----|----|----|----|----|----|
***************************************************
[2023-Dec-10 07:19:49] Read-depth annotation

0%   10   20   30   40   50   60   70   80   90   100%
|----|----|----|----|----|----|----|----|----|----|
***************************************************
[2023-Dec-10 08:25:23] Genotyping

0%   10   20   30   40   50   60   70   80   90   100%
|----|----|----|----|----|----|----|----|----|----|
***************************************************
[2023-Dec-10 08:25:25] Library statistics
Sample: RD.rmdup1


Am I doing anything wrong?
Thanks

@tobiasrausch
Copy link
Member

Probably easier to output the BCF:

delly call -o ${SAMPLE_ID}.bcf -g ${FASTA} ${INPUT_DIR}/${SAMPLE_ID}.rmdup1.bam

@tobiasrausch tobiasrausch self-assigned this Dec 12, 2023
@Ana1335
Copy link
Author

Ana1335 commented Dec 12, 2023

Thanks @tobiasrausch. Yes, it is running now but my files are so big and I am doing in a for loop, ie looping through each SV type

for name in DEL DUP INV TRA INS; 
    do delly call -t ${name} -o ${OUT_DIR}/${SAMPLE_ID}_${name}_delly.bcf -g ${FASTA} ${INPUT_DIR}/${SAMPLE_ID}.rmdup1.bam   
done


It is running since yesterday but there is no output file yet. Should I expect this? (I mean should I expect to see the output file when the run is compelete?). Thanks a lot

@tobiasrausch
Copy link
Member

Sure, once the run is complete there should be a BCF output file for each SV type.

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

No branches or pull requests

2 participants