Skip to content
This repository was archived by the owner on May 5, 2020. It is now read-only.

Speed up consensus iteration #24

Open
andreas-wilm opened this issue Aug 15, 2018 · 1 comment
Open

Speed up consensus iteration #24

andreas-wilm opened this issue Aug 15, 2018 · 1 comment
Assignees

Comments

@andreas-wilm
Copy link

Use plain consensus rules for variant calling in consensus update. This will be much faster and we can live with a few FP / FN, since a final LoFreq call is done eventually anyway

@andreas-wilm andreas-wilm self-assigned this Aug 15, 2018
@andreas-wilm
Copy link
Author

The following as taken from http://samtools.github.io/bcftools/howtos/consensus-sequence.html, with --ploidy 1 added to call should work:

bcftools mpileup -Ou -f reference.fa alignments.bam | bcftools call -mv -Oz -o calls.vcf.gz
tabix calls.vcf.gz
cat reference.fa | bcftools consensus calls.vcf.gz > consensus.fa

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

No branches or pull requests

1 participant