SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
chip-seq analysis against E.coli tujchl Bioinformatics 3 10-27-2015 05:47 AM
E coli de novo sequencing SHB General 5 10-28-2011 03:47 AM
Escherichia coli O104:H4 + 11 strains assemblies What_Da_Seq Pacific Biosciences 0 07-06-2011 06:20 PM
Illumina for resequencing E. coli jkersh Illumina/Solexa 16 11-11-2010 12:05 PM
Sequencing of E. coli strains pDNA 454 Pyrosequencing 2 05-12-2010 12:45 PM

Reply
 
Thread Tools
Old 05-14-2012, 06:15 AM   #1
aggp11
Member
 
Location: Wisconsin

Join Date: Jun 2011
Posts: 87
Default Mapping E.coli PE reads

Hi,

I have to align PE reads from the E.coli DH10B sequencing. I am thinking of using BWA sampe. Could E.coli being a circular genome cause any alignment issues? If yes, then I don't see an option to for aln/sampe that could handle this. So is anybody aware of doing PE E.coli data alignment using BWA? Also, are there any better alignment tools for E.coli PE data?

Any help regarding this would be appreciated.

Thanks,
Praful

Last edited by aggp11; 05-14-2012 at 07:47 AM.
aggp11 is offline   Reply With Quote
Old 05-14-2012, 07:09 AM   #2
maubp
Peter (Biopython etc)
 
Location: Dundee, Scotland, UK

Join Date: Jul 2009
Posts: 1,543
Default

Quote:
Originally Posted by aggp11 View Post
Could E.coli being a circular genome cause any alignment issues?
Yes, because BWA won't map any reads over the origin, and could potentially reject paired reads spanning the origin.
maubp is offline   Reply With Quote
Old 05-14-2012, 08:49 AM   #3
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

You will have some reads that cross from the "end" of your genome to the origin of your genome, and those won't map. And you will have reads where one read is at the "end", and the other is after the origin, and these will be marked as improperly paired.

But it's not going to make much difference in the grand scheme of things.
swbarnes2 is offline   Reply With Quote
Old 05-14-2012, 09:48 AM   #4
aggp11
Member
 
Location: Wisconsin

Join Date: Jun 2011
Posts: 87
Default

Thanks a lot for your prompt replies. This makes sense and this is what I thought might happen, but just wasn't sure. So based on this, I would assume BWA aln/sampe should still work for most of the reads and I'll just have a slightly bigger number of improperly paired reads.

Thanks again guys.
Praful
aggp11 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:47 AM.


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