SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Samtools Flag to identify outward facing (mate-pair not paired end) reads lesander Bioinformatics 5 10-28-2011 01:24 PM
Finding reads where the mate is unmapped Ash Bioinformatics 6 09-15-2011 06:03 AM
Improperly paired mates spikesd17 Bioinformatics 2 02-07-2011 11:17 AM
Finding unmapped reads using samtools Ash Bioinformatics 2 10-28-2010 07:20 AM
matching unmapped paired SOLiD reads smarkel General 6 10-09-2009 11:49 AM

Reply
 
Thread Tools
Old 10-31-2011, 01:29 PM   #1
newguy
Junior Member
 
Location: usa

Join Date: Oct 2011
Posts: 1
Default Finding unmapped paired-end mates with samtools

I believe the reference sequence is incorrect in my region of interest. I have access to alignment files for several whole genomes and I'd like to use Velvet to do a targeted de novo assembly of just this region. In order to take full advantage of the paired-end reads, I want to compile a file with all reads for which one of the paired-end mates is mapped to the region (I would like the mapped reads, plus any unmapped paired-end mates that don't align to the faulty reference plus any paired-end mates that align elsewhere in the genome). I can use samtools view to identify the mapped reads, but the mates have the same read name, and are only differentiated by 64 vs 132 in the flag code. What is the best way to get the mates that are unmapped or mapped elsewhere?
newguy is offline   Reply With Quote
Old 10-31-2011, 02:14 PM   #2
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

SAMtools specs call for the unmapped mate of a mapped read to be given the position of the mapped end.

So if you think that the area from 10,000-12,000 is not right, get all the reads whose coordiantes are 10,000-12,000, and that'll include all the reads that really do map there, as well as the unmapped mates of any of those reads.

As for reads where one end maps there, and the other maps elsewhere, I don't think you can do that directly with samtools. But once you get a list of all the names of reads that map to your region, you can get all of the reads with that name out of your fastq, and that's really the set of reads you want anyway.

Last edited by swbarnes2; 10-31-2011 at 02:19 PM.
swbarnes2 is offline   Reply With Quote
Old 11-01-2011, 04:04 PM   #3
kjlee
Member
 
Location: Atlanta, GA

Join Date: Jun 2011
Posts: 12
Default

Hey newguy,

I am having the same issue. For me, I am looking for insertions (specifically mobile element insertions) that have occurred in the re-sequenced genome with respect to the reference.

I ran RetroSeq, and it gave me a bunch of points that it has called as points of insertion of mobile elements. What I want to do is extract all of the reads around that (potential) break-point and do a local assembly. I appreciate swbarnes2's input, but I, perhaps like newguy, was hoping that there was a tool already written for this purpose.

swbarnes2: Thank you for your suggestion; you have guided my coding. (It probably can be done with a simple grep.)

Cheers,

Kevin
kjlee is offline   Reply With Quote
Old 09-11-2012, 11:04 AM   #4
liuseq
Junior Member
 
Location: NY

Join Date: Nov 2010
Posts: 2
Default Mobile element insertion detection question

me too, if you have the experiance of how to deal with mobile element insertion or known scripts to detect that, please share with me, I am interested it too. I have data from complete genomics, they have generated MEI files for sequencing data, but I dont know how to summary those novel MEI.

thanks



Quote:
Originally Posted by kjlee View Post
Hey newguy,

I am having the same issue. For me, I am looking for insertions (specifically mobile element insertions) that have occurred in the re-sequenced genome with respect to the reference.

I ran RetroSeq, and it gave me a bunch of points that it has called as points of insertion of mobile elements. What I want to do is extract all of the reads around that (potential) break-point and do a local assembly. I appreciate swbarnes2's input, but I, perhaps like newguy, was hoping that there was a tool already written for this purpose.

swbarnes2: Thank you for your suggestion; you have guided my coding. (It probably can be done with a simple grep.)

Cheers,

Kevin
liuseq is offline   Reply With Quote
Reply

Tags
paired-end, samtools

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 02:38 PM.


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