Would it be possible to add an outgroup masking feature to snptoolkit combine? For example, the outgroup is processed along with the study samples and included in the combine command like the study samples via VCF and bam files, but it is specified with an --outgroup flag and only SNP positions called in the study samples get called in the outgroup sample for the multiple sequence SNP alignment. The alignment size can get very large in studies with highly divergent outgroups, so this would help with tree building efficiency.
This is not an issue, but a request/suggestion.
Would it be possible to add an outgroup masking feature to snptoolkit combine? For example, the outgroup is processed along with the study samples and included in the combine command like the study samples via VCF and bam files, but it is specified with an --outgroup flag and only SNP positions called in the study samples get called in the outgroup sample for the multiple sequence SNP alignment. The alignment size can get very large in studies with highly divergent outgroups, so this would help with tree building efficiency.
Cheers, Kelly