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!
1 answer
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.
Log in to answer this question.
-
FilterMutectCalls producing empty vcf file
written by ChumBucket2024 •
-
VG giraffe autoindex and alignement mapping warning "Refusing to perform too-large rescue alignment"
written by s_135 •
-
How to make Galaxy identify an adaptor when it doesn't? What is the possible source of errors?
written by Sanjukta •
-
Affymatrix microarray data do not have normal distribution after normalization (left skewed)
written by firstorthopedicdoctor •
-
Determining reason for low quality assembly
written by liorglic
-
Encountering errors with PAGIT icorn
written by Ace
-
Picard tools: error with CollectInsertSizeMetrics
written by anais1396 •
-
HT-seq Counts Very High Number of No Feature
written by Krish Srinivasan •
-
Illumina sequencing fail
written by gerberd1990 •
-
Help With Bioperl Remote Blast Script (Bio::Tools::Run::Remoteblast)
written by campbell.malachy •
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.
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.
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?