SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Ion Torrent



Similar Threads
Thread Thread Starter Forum Replies Last Post
unique mapped reads from Mapsplice yuhj RNA Sequencing 1 09-25-2011 04:40 PM
Regarding Unique reads, Unique alignments sridharacharya RNA Sequencing 2 09-20-2010 05:39 AM
SHRiMP - how to obtain unique mapped reads? idonaldson Bioinformatics 10 07-02-2010 02:48 PM
Unique VS Non-Unique read analysis samt Bioinformatics 2 09-29-2009 09:44 AM
Unique mapped reads definition confusing... edge Bioinformatics 6 09-10-2009 12:44 AM

Reply
 
Thread Tools
Old 12-03-2013, 06:10 PM   #1
gigigou
Member
 
Location: Nanjing,CHINA

Join Date: May 2012
Posts: 31
Default How to get unique mapped with TAMP?

Hi ,all.

I'm using TMAP to align Ion Torrent reads. I wonder how I can get unique mapped reads from the TMAP output SAM files.

I find the -a parameter, it seems that to set -a 0 will produce unique map. But from the description I think that -a 0 means if there are multiple alignments TMAP will chose the one with best score, obviously this is not unique.

So if anyone can help me with this?

Many thanks!
gigigou is offline   Reply With Quote
Old 12-04-2013, 12:21 AM   #2
Michael.Ante
Senior Member
 
Location: Vienna

Join Date: Oct 2011
Posts: 123
Default

Hi gigigou,

After running TMAP you have to count the mapped reads which mapped only once.
One straight-forward (but maybe inefficient) approach is to count them with awk on the command line:
Code:
samtools view tmap.sort.bam | awk '!match($3,/\*/){t[$1]++}END{for(i in t){if(t[i]==1){print i}}}'
This will give you a list of read-IDs which aligned to a chromosomal location (field 3 must not be '*') and occured only once. You can use this list to filter your bam/sam-file.

I hope this will help you.
Michael.Ante is offline   Reply With Quote
Old 12-04-2013, 07:09 PM   #3
gigigou
Member
 
Location: Nanjing,CHINA

Join Date: May 2012
Posts: 31
Default

Quote:
Originally Posted by Michael.Ante View Post
Hi gigigou,

After running TMAP you have to count the mapped reads which mapped only once.
One straight-forward (but maybe inefficient) approach is to count them with awk on the command line:
Code:
samtools view tmap.sort.bam | awk '!match($3,/\*/){t[$1]++}END{for(i in t){if(t[i]==1){print i}}}'
This will give you a list of read-IDs which aligned to a chromosomal location (field 3 must not be '*') and occured only once. You can use this list to filter your bam/sam-file.

I hope this will help you.
Thank you!

I'll try it and tell you if it works!

Many thanks.
gigigou is offline   Reply With Quote
Reply

Tags
tmap

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 01:20 AM.


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