unique hits

arrchi · 05-26-2011, 12:55 PM

Hi there,

I am looking for a way to find unique hits for my RNA-seq data. After searching online and in this community, I still can't find a good way to find unique hits from a sam file. Any input will be welcome.

I used tophat to align the data to the HG19 and samtools -bq -1 to generated reliable hits.

Here is part of the output:

[Tophat_out]$ grep -w SRR087416.97659 accepted_hits_realiable.sam
SRR087416.97659 0 chr1 11356 1 36M * 0 0 CAGCTAGGGACATTGCAGGGTCCTCTTGCTCAAGGT BBBCBCB9CBBBC@:+>3A97AACABA9@CCCCB9# NM:i:0 NH:i:4 CC:Z:chr12 CP:i:94218
SRR087416.97659 16 chr12 94218 1 36M * 0 0 ACCTTGAGCAAGAGGACCCTGCAATGTCCCTAGCTG #9BCCCC@9ABACAA79A3>+:@CBBBC9BCBCBBB NM:i:0 NH:i:4 CC:Z:chr15 CP:i:102519779
SRR087416.97659 16 chr15 102519779 1 36M * 0 0 ACCTTGAGCAAGAGGACCCTGCAATGTCCCTAGCTG #9BCCCC@9ABACAA79A3>+:@CBBBC9BCBCBBB NM:i:0 NH:i:4 CC:Z:chr2 CP:i:114359624

My questions are: Does this mean that the read SRR087416.97659 maps to chr1, chr12, and chr15?

If I understand the concept "unique-hit" correctly, then this read can not be counted as an unique hit. Am I right?

Thanks,

-A

rnaeye · 05-29-2011, 02:57 PM

Quote:

Originally Posted by arrchi

Hi there,

I am looking for a way to find unique hits for my RNA-seq data. After searching online and in this community, I still can't find a good way to find unique hits from a sam file. Any input will be welcome.

I used tophat to align the data to the HG19 and samtools -bq -1 to generated reliable hits.

Here is part of the output:

[Tophat_out]$ grep -w SRR087416.97659 accepted_hits_realiable.sam
SRR087416.97659 0 chr1 11356 1 36M * 0 0 CAGCTAGGGACATTGCAGGGTCCTCTTGCTCAAGGT BBBCBCB9CBBBC@:+>3A97AACABA9@CCCCB9# NM:i:0 NH:i:4 CC:Z:chr12 CP:i:94218
SRR087416.97659 16 chr12 94218 1 36M * 0 0 ACCTTGAGCAAGAGGACCCTGCAATGTCCCTAGCTG #9BCCCC@9ABACAA79A3>+:@CBBBC9BCBCBBB NM:i:0 NH:i:4 CC:Z:chr15 CP:i:102519779
SRR087416.97659 16 chr15 102519779 1 36M * 0 0 ACCTTGAGCAAGAGGACCCTGCAATGTCCCTAGCTG #9BCCCC@9ABACAA79A3>+:@CBBBC9BCBCBBB NM:i:0 NH:i:4 CC:Z:chr2 CP:i:114359624

My questions are: Does this mean that the read SRR087416.97659 maps to chr1, chr12, and chr15?

If I understand the concept "unique-hit" correctly, then this read can not be counted as an unique hit. Am I right?

Thanks,

-A

You are correct: SRR087416.97659 sequence maps more than one location in the genome. This means that you have no way of knowing where that sequence is coming from. You need to sort your output by sequence ID, then find uniq IDs. Does this help? Unique read should hit the genome only once at one specific location.

rnaeye · 05-29-2011, 02:59 PM

btw, what instrumentation this output is coming from? Thanks.

arrchi · 05-31-2011, 01:52 PM

Thanks for your message.

The result is generated by a Mac Pro with 8GB Memory and 1TB hard drive. Is that you want to know?

rnaeye · 05-31-2011, 01:56 PM

my questions was what DNA sequencing platform this output is from, such as Illumina, ABI SOLiD, etc. thanks.

arrchi · 06-01-2011, 06:11 AM

Oh. Sorry.

This is Illumina RNA sequencing data.

arrchi · 06-01-2011, 08:17 AM

Does anybody knows that the value 0 in "SRR087416.97659 0" means? I checked samtools menu, it only says if 0x1 is unset, no assumptions can be made about .....

arrchi · 06-01-2011, 08:53 AM

I found an old post saying that

Quote:

Flag 0 means "the read is not paired and mapped, forward strand".

Hope it is true.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
how to extract unique hits from a sam file	gfmgfm	Bioinformatics	10	07-07-2019 08:33 AM
Bowtie call to get unique, multi-hits and nonmatching reads	PFS	Bioinformatics	3	07-07-2019 08:18 AM
Regarding Unique reads, Unique alignments	sridharacharya	RNA Sequencing	2	09-20-2010 06:39 AM
Unique VS Non-Unique read analysis	samt	Bioinformatics	2	09-29-2009 10:44 AM
Unique and repeat Hits	pfranchini	Bioinformatics	1	07-15-2009 09:51 AM

05-26-2011, 12:55 PM	#1
arrchi Member Location: ma Join Date: Mar 2011 Posts: 46	unique hits Hi there, I am looking for a way to find unique hits for my RNA-seq data. After searching online and in this community, I still can't find a good way to find unique hits from a sam file. Any input will be welcome. I used tophat to align the data to the HG19 and samtools -bq -1 to generated reliable hits. Here is part of the output: [Tophat_out]$ grep -w SRR087416.97659 accepted_hits_realiable.sam SRR087416.97659 0 chr1 11356 1 36M * 0 0 CAGCTAGGGACATTGCAGGGTCCTCTTGCTCAAGGT BBBCBCB9CBBBC@:+>3A97AACABA9@CCCCB9# NM:i:0 NH:i:4 CC:Z:chr12 CP:i:94218 SRR087416.97659 16 chr12 94218 1 36M * 0 0 ACCTTGAGCAAGAGGACCCTGCAATGTCCCTAGCTG #9BCCCC@9ABACAA79A3>+:@CBBBC9BCBCBBB NM:i:0 NH:i:4 CC:Z:chr15 CP:i:102519779 SRR087416.97659 16 chr15 102519779 1 36M * 0 0 ACCTTGAGCAAGAGGACCCTGCAATGTCCCTAGCTG #9BCCCC@9ABACAA79A3>+:@CBBBC9BCBCBBB NM:i:0 NH:i:4 CC:Z:chr2 CP:i:114359624 My questions are: Does this mean that the read SRR087416.97659 maps to chr1, chr12, and chr15? If I understand the concept "unique-hit" correctly, then this read can not be counted as an unique hit. Am I right? Thanks, -A

05-29-2011, 02:59 PM	#3
rnaeye Member Location: East Cost Join Date: May 2011 Posts: 79	btw, what instrumentation this output is coming from? Thanks.

05-31-2011, 01:52 PM	#4
arrchi Member Location: ma Join Date: Mar 2011 Posts: 46	Thanks for your message. The result is generated by a Mac Pro with 8GB Memory and 1TB hard drive. Is that you want to know?

05-31-2011, 01:56 PM	#5
rnaeye Member Location: East Cost Join Date: May 2011 Posts: 79	my questions was what DNA sequencing platform this output is from, such as Illumina, ABI SOLiD, etc. thanks.

06-01-2011, 06:11 AM	#6
arrchi Member Location: ma Join Date: Mar 2011 Posts: 46	Oh. Sorry. This is Illumina RNA sequencing data.

06-01-2011, 08:17 AM	#7
arrchi Member Location: ma Join Date: Mar 2011 Posts: 46	Does anybody knows that the value 0 in "SRR087416.97659 0" means? I checked samtools menu, it only says if 0x1 is unset, no assumptions can be made about .....