Generating multiple species consensus alignments for tree topology tests
0
0
Entering edit mode
3 months ago
mglasena ▴ 30

I have generated a VCF file with 9 samples, each being a separate species aligned to a common reference. I am hoping to generate multiple sequence alignment files in fasta format for each gene for use with iqtree tree topology tests. For example, for each gene in the reference, I'd like to produce something like the following:

>sample/species_1
AATT
>sample/species_2
ATTT
>sample_species_3
TTTT

I am assuming bcftools consensus is the tool to use for this. I read a related tutorial using this program. However, bcftools consensus does not accommodate multi-sample vcf files. What would be the most efficient way to create a multiple sequence alignment file for each region in a BED file?

bcftools alignment fasta iqtree vcf • 166 views
ADD COMMENT
0
Entering edit mode

Have you tried play with bedtools?

ADD REPLY

Login before adding your answer.

Traffic: 2604 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6