SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Reads mapping of ends of chromosomes ETHANol Bioinformatics 0 11-03-2011 09:09 AM
Paired-end reads mapping to different chromosomes... rareaquaticbadger Bioinformatics 1 06-27-2011 04:24 AM
How many functionally-assigned reads is enough? idoerg Bioinformatics 0 04-25-2011 07:44 AM
What are the unmapped reads beelu Illumina/Solexa 1 09-09-2010 05:18 AM
tophat paired reads mapped to different chromosomes wenhuang Bioinformatics 3 06-23-2010 05:23 AM

Reply
 
Thread Tools
Old 12-01-2011, 10:54 AM   #1
efoss
Member
 
Location: Seattle

Join Date: Jul 2011
Posts: 98
Default unmapped reads assigned to chromosomes

If I use samtools' "idxstats" commands, it returns information with:

1. reference sequence name (chromosome name in my case)

2. length of reference sequence

3. number of mapped reads assigned to that reference sequence

4. number of unmapped reads assigned to that reference sequence

Number 4 makes no sense to me. If a read is unmapped then by definition you can't assign it to a genomic location. If I remove everything with map quality less than 1, most of these go away but not all of them, i.e. I still end up with unmapped reads assigned to chromosomes. Does anyone understand this?

Thank you.

Eric
efoss is offline   Reply With Quote
Old 12-01-2011, 11:08 AM   #2
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

Come on. If you aren't going to read the SAM format specs before you ask a question about the SAM format, at least search the message board before asking a question. I think I've answered this question three times this week alone.
swbarnes2 is offline   Reply With Quote
Old 12-01-2011, 11:16 AM   #3
efoss
Member
 
Location: Seattle

Join Date: Jul 2011
Posts: 98
Default

Quote:
Originally Posted by swbarnes2 View Post
Come on. If you aren't going to read the SAM format specs before you ask a question about the SAM format, at least search the message board before asking a question. I think I've answered this question three times this week alone.
Sorry. I had searched the message board unsuccessfully, but yes - it was stupid not to look at the sam documentation more closely:

4.1 For a unmapped paired-end or mate-pair read whose mate is mapped, the unmapped read should have RNAME and POS identical to its mate.

Eric
efoss is offline   Reply With Quote
Old 12-01-2011, 01:22 PM   #4
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

Yup. The other possibility is if you use bwa for alignments. bwa concatenates all the reference sequences together, so if a read hangs off of one reference onto another, it will be given the right chromosome and position, but the 4 flag will also be put in, as a sign that something is not quite right. I'm not sure if other aligners, like bowtie, will do that.
swbarnes2 is offline   Reply With Quote
Old 12-01-2011, 01:24 PM   #5
efoss
Member
 
Location: Seattle

Join Date: Jul 2011
Posts: 98
Default

Quote:
Originally Posted by swbarnes2 View Post
Yup. The other possibility is if you use bwa for alignments. bwa concatenates all the reference sequences together, so if a read hangs off of one reference onto another, it will be given the right chromosome and position, but the 4 flag will also be put in, as a sign that something is not quite right. I'm not sure if other aligners, like bowtie, will do that.
Thanks very much. That's good to know, since I am aligning with bwa.

Eric
efoss is offline   Reply With Quote
Reply

Tags
unmapped reads

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 11:33 AM.


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