[adegenet-forum] different number of individuals in fastas
Katharine Coykendall
dollykc at gmail.com
Thu Feb 22 21:34:15 CET 2018
Hello,
I've been following the instructions found at
http://popgen.nescent.org/PopDiffSequenceData.html to do a hierarchical F
analysis on my data. I have two fasta files with the same 117 individuals
in each file. The files are different genes. When I put them together
using the read.multiFASTA command, it tells me that there are 120
sequences. When I turn it into a genid object it says there are 120
individuals and 19 alleles. I checked both files to make sure that the
sample names are the same between them and they are. I was wondering if it
had to do with add.gaps=TRUE so I tried to set it to false, but I get an
error
Error in file(con, "rb") : cannot open the connection
In addition: Warning message:
In file(con, "rb") : cannot open file 'FALSE': No such file or directory
Is there an easy way to look at the multiFASTA or genid file to see where
the extra three individuals are coming from?
Thanks,
Katharine
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/adegenet-forum/attachments/20180222/0324e40c/attachment.html>
More information about the adegenet-forum
mailing list