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!
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!
Comment