SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > SOLiD



Similar Threads
Thread Thread Starter Forum Replies Last Post
target mapped read percentage eren Illumina/Solexa 0 08-13-2011 10:06 PM
GAII low number of mapped reads aligenie Bioinformatics 3 06-21-2011 11:55 PM
expected percentage of mapped reads in chip-seq experiment andreiafonseca General 16 03-25-2011 05:14 AM
Percentage of mapped reads ? zack80.liu Bioinformatics 6 03-01-2011 09:08 AM
SOLiD SAGE: low percentage of reads mapped rahilsethi SOLiD 0 09-09-2010 11:04 AM

Reply
 
Thread Tools
Old 09-09-2010, 10:41 AM   #1
rahilsethi
Member
 
Location: Pittsburgh, PA

Join Date: May 2010
Posts: 22
Post low percentage of reads mapped

Hello,

I was analyzing the SOLiD SAGE data of human sequence using SOLiD SAGE Analysis tool. I performed Mapping using 27 bp length with 1 mismatch. The reference was the complete set of human mRNA sequence from Refseq db. I then calculated % of reads that mapped to reference using the results file that gives the list of tags and their corresponding read files. I ran the analysis for 4 SAGE data and I got the following percentage:
SAGE A : 15 % apporx.; SAGE B 16% approx.; SAGE C 17% and SAGE D 20 %

What can be the reason for such low percentage of reads mapping to human mRNA reference? Is it a general result for most of the SOLiD SAGE experiments?

Thanks

Last edited by rahilsethi; 09-09-2010 at 10:46 AM.
rahilsethi is offline   Reply With Quote
Old 09-09-2010, 11:22 AM   #2
NextGenSeq
Senior Member
 
Location: USA

Join Date: Apr 2009
Posts: 482
Default

BLAT some of the unmapped reads to see what they are. If they are actually real mRNA, its your analysis. If they are genomic DNA or something else, it's your library.
NextGenSeq is offline   Reply With Quote
Old 09-09-2010, 12:06 PM   #3
rahilsethi
Member
 
Location: Pittsburgh, PA

Join Date: May 2010
Posts: 22
Default

What it has to do with my analysis? The result is straight from the software SOLiD SAGE Analysis tool. In the result.tab file it produces read ids for the tags are mentioned. I counted the unique set from those read ids and divided it by the unique set of read ids from the read file to get the percentage of reads mapped to human mRNA obtained from Refseq database.
I will still BLAT some of them and see where they are mapping
If they are not mapping to mRNA then something should be with the reads generated by SOLiD SAGE run
rahilsethi is offline   Reply With Quote
Old 09-13-2010, 06:01 AM   #4
pmiguel
Senior Member
 
Location: Purdue University, West Lafayette, Indiana

Join Date: Aug 2008
Posts: 2,317
Default

Hi Rahlisethi,

NextGenSeq is giving you a "sanity check" to help with your troubleshooting. The SAGE might be working fine but your RNA might contain lots of transcripts from repetitive elements, for example. Then their position is not uniquely mappable in the genome. Or you might get hits to E. coli, or some other unexpected species, and suspect contamination from some source.

--
Phillip
pmiguel 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 09:20 AM.


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