SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Assembler for PE-300 genomic reads AnjaF Illumina/Solexa 0 10-20-2014 02:51 AM
de novo assembler for se 36 bp illumina reads srividyanathan De novo discovery 2 05-28-2013 07:52 AM
HiSeq 2000 50PE or 100PE? luckylove Illumina/Solexa 2 01-20-2013 09:14 PM
What's the best assembler for 454 reads? AFBrito 454 Pyrosequencing 3 05-31-2012 04:38 AM
Looking for assembler to align very short reads andylai Bioinformatics 0 06-02-2011 08:58 PM

Reply
 
Thread Tools
Old 02-27-2015, 08:46 AM   #1
horvathdp
Member
 
Location: Fargo

Join Date: Dec 2011
Posts: 66
Default best assembler for 160PE and 100PE reads

Hi all,
I have about 60X combined coverage of my genome, but half of my reads are 100 base paired end reads and the other half is 160 base paired end reads. Any suggestions for which assembler would be the best one to use for a combined file of both read sets?

I have removed the high and low kmer read sets (since I really am primarily interested in single and low-copy regions of my genome) and total I have about 20GB of sequence for each end. I ran the Trinity Kmer reduction program (max of 30 X coverage/kmer) on my PE100 file and only reduced it from about 8GB to 5GB per end. This is about what I would have expected given my initial kmer selection scheme. The attempt to do this with the PE160 files failed for an unknown reason, but I would expect it to also would give me about a 40% reduction. Thus, I could probably -with some effort, get my combined files down to about 12Gb each before attempting the assembly.

Any advice about the best assembler or other strategies to assist in assembling my low-copy genome space (which should contain the bulk of my genes) is welcome.
horvathdp is offline   Reply With Quote
Reply

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 03:01 PM.


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