Question: DELLY produces empty output
0
gravatar for liorglic
12 months ago by
liorglic340
liorglic340 wrote:

Hello,
I am trying DELLY v0.8.2 for the first time, using ~x40 coverage data on a plant genome. The bam file was obtained by running bwa mem, then sorting and indexing the result.
I am using default params for DELLY. My run results in a bcf file with no data (i.e. header is present, but no variants detected) and prints the following message to stdout:

Sample:SRR1945819_map_to_pan.sort,ReadSize=51,Median=407,MAD=11,UniqueDiscordantPairs=230596

I find it very unlikely that there are no SVs between my sample and the reference, so I guess something is wrong with my data or the way I use DELLY.
Any ideas on what might be going wrong?
Thanks!

cnv delly sv • 288 views
ADD COMMENTlink modified 12 months ago • written 12 months ago by liorglic340

You may consider posting on the GitHub issues page: https://github.com/dellytools/delly/issues

If you take that option, link back to this thread.

ADD REPLYlink written 12 months ago by Kevin Blighe71k

DELLY has a restriction on mapping quality - and your reads are 51bp (if I am correct), so some problems with multi-mappers and thus "hidden" SVs are expected, but surely not empty output. I would try other callers such as Manta and if nothing works, write to the author of DELLY.

ADD REPLYlink written 12 months ago by German.M.Demidov1.9k

Tanks @German.M.Demidov - I am getting the same result with another bam file which has lower depth (~x10) but longer reads (101). Interestingly, I've just tried Manta on the same data and got an empty vcf again! I am starting to wonder what's wrong with this data. Perhaps PE/SR approaches are incapable of dealing with them for some reason. Any idea on what to try next?

ADD REPLYlink written 12 months ago by liorglic340
4
gravatar for liorglic
12 months ago by
liorglic340
liorglic340 wrote:

OK, this is a bit embarrassing, but turns out the problem was not with DELLY (or Manta), but rather with bcftools. Turns out I was using a very old version (0.1.19) that apparently can't read the new bcf. Once upgraded to 1.10.2, I am now able to read the output bcf, which is not empty at all.

ADD COMMENTlink written 12 months ago by liorglic340
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1692 users visited in the last hour
_