SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Illumina/Solexa



Similar Threads
Thread Thread Starter Forum Replies Last Post
Draft genome scaffolding with RNAseq paired-end reads JPZ Bioinformatics 7 02-10-2015 02:05 AM
Haplotype phasing from Miseq 150bp paired end reads? michael-steffen Bioinformatics 3 01-08-2015 11:41 PM
Assembly MiSeq Paired End PacMan Bioinformatics 12 07-14-2014 04:57 AM
Recommended aligner to use with a draft reference genome and paired end reads N311V Bioinformatics 2 07-16-2013 05:32 PM
250bp paired end Miseq reads alignment empyrean Bioinformatics 3 12-12-2012 12:36 PM

Reply
 
Thread Tools
Old 12-23-2015, 05:10 PM   #1
cktheory
Junior Member
 
Location: Pennsylvania

Join Date: Jan 2015
Posts: 3
Default MiSeq Paired End Reads for Draft Assembly

Hi Everyone,
I am trying to de novo assemble a draft bacterial genome with a single chromosome from another persons data. The insert size of the sequenced segments average to 580bp and the read length for the R1 and R2 files is 301. After filtering there is not enough overlap to merge the paired reads. Is there a way to still use this data to assemble a draft genome that can be annotated?

Last edited by cktheory; 12-23-2015 at 07:06 PM.
cktheory is offline   Reply With Quote
Old 12-23-2015, 10:00 PM   #2
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707
Default

Filtering does not change the read length. If the insert size is 580, then actually there is enough overlap to merge them - try BBMerge.

You don't need to merge them, though. You can just adapter-trim and then assemble with Spades. Merging and error-correction may give you a better assembly, but neither is required.
Brian Bushnell is offline   Reply With Quote
Reply

Tags
assembly, miseq, paired end reads

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 08:55 AM.


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