[adegenet-forum] different number of individuals in fastas

Katharine Coykendall dollykc at gmail.com
Thu Feb 22 23:02:55 CET 2018


Figured it out!  In one of my fasta files, I had three of the sample names
capitalized and the other file they weren't.

On Thu, Feb 22, 2018 at 3:34 PM, Katharine Coykendall <dollykc at gmail.com>
wrote:

> 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/698a93c0/attachment.html>


More information about the adegenet-forum mailing list