How to prep MAFFT alignment files for use in R packages ape and pegas?
Hello everyone,
I have a single alignment file that I would like to analyze using the Tajima's D test in pegas. The alignment file is in Clustal format, and was generated using MAFFT. Within my file, I have alignments for every gene in the reference genome, aligned to the rest of the genomes in the population I'm studying.
Is there any file processing I have to do in order to import the alignment file into R? I know I have to use the ape function DNAbin in order to use the Tajima's D test, but I'm not sure what needs to be done in terms of processing my alignment file to be ready for import.
If anyone could help, I would be very grateful!
• 1,600 views
•
link
written
by
juliuspasion1994 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Tajima's D calculation for specific genes in multi-sample VCF
written by Muhammad •
-
Tajima's D test applied to tandemly repeated genes
written by jaer.wees •
-
Need help with Codeml analysis
written by Manoj
-
Population Genetics from SNP data- no VCF
written by vidyavuru •
-
Pairwise alignment from Multiple alignment MAFFT file
written by brimaloney24 •
-
Calculating Tajima's D for an individual genome vs. for population of genomes vs. for genes
written by severalorks
-
Trimming an alignment in R
written by Gon •
-
Pegas import fasta and calculate theta
written by Adrian Pelin
-
detection of positive selection in population using Tajima's n Fu Li test
written by evo_genomics •
-
Tajima'S D Test-Ngs Data
written by jackuser1979
