RADseq chromosomes unknown
0
0
Entering edit mode
@clivedarwell-9195
Last seen 5.9 years ago
Japan

Hi all

 

I have a RADseq data set (from Stacks) which has no chromosome identifiers for the SNP calls. When importing the VCF file I get an error message due to this. So, I converted the "unk" chromosome information in the VCF into 1s. I can now load the file. When I filter out the linked SNPs I get around 400 SNPs remaining out of 40000. My PCA analysis shows virtually no structure when we were expecting strong signal and the top principal components are low and do not drop significantly after the first two (they are all around 5%). If I run the data without filtering out the linked SNPs I get a similar result which seems even more surprising.

 

Does anyone have any thoughts on the source of these odd results?

 

Thanks

 

Clive

 

snprelate RADseq Stacks • 582 views
ADD COMMENT

Login before adding your answer.

Traffic: 193 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6