Comparing loci across catalogs
Entering edit mode
17 months ago
Austin • 0

Hi all,

Context: I'm using the Stacks (v2.59) command to analyze garden and wild samples of Quercus acerifolia, and have aligned my paired-end NextRAD reads (filtered using process_radtags) to a reference genome (Q. robur) using GSNAP. The goal is to measure how many wild alleles are observed (captured) within garden samples. To this end, I want to generate a catalog of wild loci, a separate catalog of garden loci, and then measure the proportion of wild loci observed in the garden catalog.

My question is: how can I match RAD loci that were generated in separate catalogs? For instance, can I use adegenet to read in the files generated by Stacks and match the names of loci, since both garden and wild samples are aligned to the same reference (using the same parameters)?

My preference is to do this in R (ideally with adegenet, but open to other approaches). If not, is there a way I could just compare across 2 loci catalogs (i.e. find matches between two .fa.gz files)?

Thanks for any tips!

Stacks reference alignment RADseq adegenet • 332 views
Entering edit mode

Hi @Austin, Did you ever figure this out?


Login before adding your answer.

Traffic: 1626 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6