Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
ddRAD, STACKS, and Trimming Carcharodon Bioinformatics 0 07-12-2017 11:09 AM
ddRADSeq data analysis for population structure (Stacks, dDocent...) yarinka Bioinformatics 8 03-23-2015 02:01 PM
Gviz problem---too many stacks to draw angel-sakura Bioinformatics 1 01-04-2015 01:20 AM
Stacks exce_velvetg not present bio_jit Bioinformatics 0 07-02-2014 09:43 PM
Stacks process_radtags N in overhang CGO Bioinformatics 0 06-18-2014 06:17 AM

Thread Tools
Old 07-13-2017, 05:09 PM   #1
Location: Honolulu, HI

Join Date: Jul 2015
Posts: 40
Default STACKS and PE Illumina Data (ddRAD)

I have a dataset from a ddRAD library (Illumina HiSeq4K, 150PE). I'm playing around with STACKS, and noticed that it has some limitations when dealing with paired-end data from ddRAD datasets.

Namely, it treats them as separate/independent loci.

This seems like a problem from a population genomic point of view, since there is a base assumption (admittedly often violated) of independence between loci. If we have paired data (single-reads and paired-end ("PE") reads), we know that those two loci aren't at all likely to be independent.

What I have done so far is demultiplex individuals' data based on inline barcode sequence(s) and give a rough quality filter (sliding-window 15%, min quality score = 10). So I'm left with four files for each individual:

One file of SE reads and one file of PE reads, in-phase (where both reads from the same fragment/cluster were kept and not discarded).

One file of SE reads, whose PE counterpart has been discarded.

One file of PE reads, whose SE counterpart has been discarded.

Only one or the other read - SE or PE - needs to have a SNP. The other read can be discarded. Can STACKS keep track of this? Can it go through assembly and SNP-calling and keep track of header-titles, and use that information to figure out what those sequence pairs are?

Or is that information lost? (I suspect that for the remainders - the latter two files described - it would be quite difficult to recover that information...)

Alternatively, is it worth throwing caution to the wind and using the SE and PE reads and throwing all of those data together at the end? Is there another approach here?

Many thanks,
Carcharodon is offline   Reply With Quote

paired end, paired reads, stacks

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 07:27 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO