How to extracts variants from multiple sequence alignment files?
0
0
Entering edit mode
15 months ago

I have a clustalw file in .aln format I want to get the variants only from this file

clustalw biopython fasta fastq • 740 views
0
Entering edit mode
0
Entering edit mode

I try it but I got this: Alignment sequences.fasta contains sequences of unequal length. Expected length is 712 but got 688 in sequence MT364213.1

0
Entering edit mode

Did you use Pierre's script or https://github.com/sanger-pathogens/snp-sites ?

0
Entering edit mode
0
Entering edit mode

That link requires an aligned fasta format file which is not the same as .aln file you have. I suggest you use Pierre's script in this answer: Getting A Vcf File From A Fasta Alignment

0
Entering edit mode

I found that the sequences must have the same length

0
Entering edit mode

If your sequences aren't the same length, they aren't aligned, so you need to check your input data properly.