-
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
Absolutly All positions in mpileup (previously -aa) #2076
Comments
Hi, |
Samtools mpileup isn't going anywhere, but it no longer has VCF output and nor will this be added back as the decision was made to limit calling code to bcftools only for ease of maintenance. I can't comment on whether |
thank you for your response, it is reassuring to know that samtools mpileup will at least remain, as I build my pipeline around it. Perhaps |
What is the motivation for this feature request? |
Hi, thanks for getting back to me So I am working with custom scripts written by the group I am in, that "add" variant information generated with other methods (that output in fasta format) together in a VCF for downstream analysis, if this makes sense. For this, we would need a "starting" vcf that represent all positions. I realise this is probably something very few people would use, but since |
Dear developers,
I am looking to generate a VCF file in which all positions of the reference are represented - so also those where my input bam has no read coverage. In SAMtools mpileup, the option
-aa
exists which does what I need, but I understand SAMtools mpileup is depricated and shouldn't be used (?). I also am unable to pipe the output from SAMtools mpileup intoBCFtools call
for generating a VCF, have so far been unable to find another way to convert SAMtools' pileup format to a VCF. Is there a native option in BCFtools mpileup that would result the-aa
behaviour from SAMtools, ie resulting in a VCF where all bases of the reference are present?Any suggestions would be much appreciated! I'm attaching a small example file, in case that helps.
XMEN.zip
way to reproduce:
bcftools mpileup XMEN1_10.bam -f XMEN-reference_10.fasta
samtools mpileup -aa XMEN1_10.bam
Many thanks!
The text was updated successfully, but these errors were encountered: