SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
total number of counts with HTseq oliviera Bioinformatics 17 07-26-2013 07:33 AM
Relationship between kmer size and read length and total no.of reads. ? vaibhavvsk Introductions 2 03-19-2013 06:25 AM
Handmade mRNAseq library prep protocol w/ 10ng total RNA lukas1848 Sample Prep / Library Generation 0 11-02-2011 04:31 AM
total rna quality for library construction wayland Illumina/Solexa 3 10-13-2011 05:10 AM
total RNA quality for library construction neveaire Sample Prep / Library Generation 4 04-05-2011 11:39 AM

Reply
 
Thread Tools
Old 11-07-2012, 02:05 PM   #1
Anna Nad
Junior Member
 
Location: Basel

Join Date: Nov 2012
Posts: 4
Default more counts than total library size ?

Hello,
I have mRNASeq data displaying counts per splice variants of mRNA, and I was just wondering how it is possible that the total number of counts in the sample is much higher than the total number of mapped reads (and of total number of reads generated by Illumina).
thank you in advance for your answer,
anna
Anna Nad is offline   Reply With Quote
Old 11-07-2012, 02:16 PM   #2
aliceb
Member
 
Location: Switzerland

Join Date: Jan 2010
Posts: 18
Default

Hi Anna,

It sounds like some of your reads are mapping to more than one place. This way they would be reported more than once in your counts.

How are you doing the mapping?

-Alice
aliceb is offline   Reply With Quote
Old 11-07-2012, 02:37 PM   #3
Anna Nad
Junior Member
 
Location: Basel

Join Date: Nov 2012
Posts: 4
Default

Hello Alice,
The tool BWA is used to map the reads with a maximum set at 2 mismatches in the first 32 bases of the sequences, and a maximum of n mismatches in total.
In my sample, library size was about 277 million reads and mapped reads were about 192 million (69,5%). But total counts is about 861 million reads.
So mismatches were the problem? so 192 million are mapped, but it does not mean uniquely mapped ?
thank you Alice,
anna
Anna Nad is offline   Reply With Quote
Old 11-07-2012, 04:20 PM   #4
aliceb
Member
 
Location: Switzerland

Join Date: Jan 2010
Posts: 18
Default

Hi Anna,

I'm afraid I haven't used bwa, but we see similar results with bowtie2 (although not quite to the same extent as you).

When we recover counts from our sam file (we're using ht-seq-count), the number of non-unique reads is reported, which we interpret as reads that are counted more than once because they are mapping to more than one place.

I would try your alignments allowing for fewer mismatches and see what happens.

Good luck.
aliceb is offline   Reply With Quote
Old 11-07-2012, 05:04 PM   #5
timydaley
Member
 
Location: Los Angeles

Join Date: Jun 2010
Posts: 26
Default

BWA does report secondary alignments. This is likely the cause, you can remove all but the primary alignment with samtools view -F 0x0100.

Last edited by timydaley; 11-07-2012 at 05:42 PM.
timydaley is offline   Reply With Quote
Old 11-13-2012, 02:04 AM   #6
Anna Nad
Junior Member
 
Location: Basel

Join Date: Nov 2012
Posts: 4
Default

thank you Timydaley
Anna Nad is offline   Reply With Quote
Reply

Tags
counts, library, mapped, size

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:57 PM.


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