DNAbin2genind and ambiguity codes from Fasta

[burbrink]

Hello,

It seems that when I read fasta files with ambiguity codes for heterozygous states (e.g., w = A/T) using DNAbin2genind this is coded as NA in the genind object. Unfortunately, other functions doing this at the genome scale that also rely on this (e.g., multidna2genind) also pass this as NA to the genind object. Do you know of any fix for this?

Thanks so much!

Frank

[zkamvar]

Hello,

I'm not sure why, but DNAbin2genind always assumed haploid sequences for FASTA, which you can see here:

adegenet/R/sequences.R

Line 65 in 78be588

res <- genind(temp, ploidy=1, pop=pop, type="codom")

I believe that modifying DNAbin2genind with a flag that allows for heterozygous sites is possible. Would you like to make a pull request for this functionality?