-
Notifications
You must be signed in to change notification settings - Fork 241
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
bcftools mpileup call with empty output #2069
Comments
Check if there are any reads in that region
A common problem is a mismatch in chromosome naming convention ( |
The next thing to check then is what raw mpileup looks like. Maybe all positions in that region are non-variant, therefore
or
|
In the first screenshot all bases in all reads are identical to the reference, hence there are no variants, nothing to call. The second screenshot shows several bases with no coverage. There seem to be no overlapping reads, therefore the five positions with zero coverage are not called as deletion. It is possible there was a read with an indel and was filtered by mpileup (eg because of the --min-ireads option), but this I cannot tell just from seeing the screenshot. |
This requires a test case to reproduce and debug the problem. This script can be used to create a small slice of the bam and the reference https://github.com/pd3/mpileup-tests/blob/main/misc/create-bam-test |
Hello, I'm interesting in identify if my sample has a insertion, deletion or if it's heterozygote (in a specific region). To do this I ran the next code
bcftools mpileup -Ou -r chrX:start-end -f genome.fa $Input_Path/input.bam | bcftools call -Ou -mv -o test_option3
When I run it I don't get any error messages. The following are the terminal messages.
The problem is my output file
test_option3
contain only the header. The last rows are:I also tried to execute the same analysis in Galaxy and I got the same output file, Is it possible that the error was on my bam file? How can I test it?
The text was updated successfully, but these errors were encountered: