SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
target mapped read percentage eren Illumina/Solexa 0 08-13-2011 10:06 PM
expected percentage of mapped reads in chip-seq experiment andreiafonseca General 16 03-25-2011 05:14 AM
Paired-end reads mapped 50% while each pair reads mapped 80% zack80.liu Bioinformatics 3 03-03-2011 01:06 AM
low percentage of reads mapped rahilsethi SOLiD 3 09-13-2010 06:01 AM
SOLiD SAGE: low percentage of reads mapped rahilsethi SOLiD 0 09-09-2010 11:04 AM

Reply
 
Thread Tools
Old 02-28-2011, 10:15 AM   #1
zack80.liu
Junior Member
 
Location: TX

Join Date: Feb 2011
Posts: 9
Default Percentage of mapped reads ?

Hi all,

I am new to the RNA seq world. I have processed bunch of illumina paired end reads. The reported alignment is only 49%.. I was wondering if this is normal. Here's the command and output. Thanks.

bowtie -n 3 -p 10 --best -e 200 --trim5 15 --trim3 25 --sam BowtieIndexes/mm9 -1 Raw_files/sample-1_export.fq -2 sample-2_export.fq sample.sam


# reads processed: 21482797
# reads with at least one reported alignment: 10637665 (49.52%)
# reads that failed to align: 10845132 (50.48%)
Reported 10637665 paired-end alignments to 1 output stream(s)
zack80.liu is offline   Reply With Quote
Old 02-28-2011, 01:11 PM   #2
plassaaw
Junior Member
 
Location: Cincinnati, OH

Join Date: Feb 2011
Posts: 7
Default

If you could provide us with a little more information about the data you are analyzing that would be helpful. What are the reads from?
I have had some samples that have had 90+% accuracy, but with things that are likely high in indels like cancer, I have gotten as low as 60ish %..
plassaaw is offline   Reply With Quote
Old 02-28-2011, 01:37 PM   #3
zack80.liu
Junior Member
 
Location: TX

Join Date: Feb 2011
Posts: 9
Default

This is a mouse cancer dataset. However not a lot indels are expected. Some of the samples are normal controls. Any more information needed?
zack80.liu is offline   Reply With Quote
Old 02-28-2011, 02:30 PM   #4
plassaaw
Junior Member
 
Location: Cincinnati, OH

Join Date: Feb 2011
Posts: 7
Default

That's fairly low then. Try rerunning a set of the samples instead of as paired-end, as separate single reads. Also, try trimming less off of the reads, because there is a possibility that that decreases the number of possible locations for the paired ends to match up.

Let me know the results.
plassaaw is offline   Reply With Quote
Old 02-28-2011, 10:03 PM   #5
Jon_Keats
Senior Member
 
Location: Phoenix, AZ

Join Date: Mar 2010
Posts: 279
Default

If you are analyzing RNAseq data why are you using bowtie and not tophat? Your percent aligned will be low with bowtie alone as only reads aligning to an exons will map, while all reads crossing a exon-exon junction will not align. Depending on the read length this can be a large percentage or your reads.
Jon_Keats is offline   Reply With Quote
Old 03-01-2011, 06:43 AM   #6
zack80.liu
Junior Member
 
Location: TX

Join Date: Feb 2011
Posts: 9
Default

Plassaaw,

I used mm9_with_junctions index which supposedly takes care of the exon-exon junctions .

When I run single reads on these two paired_end reads _1 and _2, 80% and 78% of the reads are aligned.

However, when I run paired-end alignment, only 50% of the reads are aligned... Any idea?



# reads processed: 21482797
# reads with at least one reported alignment: 17370381 (80.86%)
# reads that failed to align: 4112416 (19.14%)
Reported 17370381 alignments to 1 output stream(s)
zack80.liu is offline   Reply With Quote
Old 03-01-2011, 09:08 AM   #7
zack80.liu
Junior Member
 
Location: TX

Join Date: Feb 2011
Posts: 9
Default

single read mapping can reach 80%

# reads processed: 34439174
# reads with at least one reported alignment: 27167675 (78.89%)
# reads that failed to align: 7271499 (21.11%)

Well paired mapping is only 52%

Any idea? Thanks!
zack80.liu 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 04:39 PM.


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