[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