SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
HTSeq dealing with "*" qualities kamsen Bioinformatics 39 10-15-2018 05:21 AM
DEXSeq error in estimateDispersions: match.arg(start.method, c("log(y)", "mean")) fpadilla Bioinformatics 14 07-03-2013 02:11 PM
HTseq count - all reads are "alignment_not_unique" dmsoanes Bioinformatics 2 07-12-2012 05:07 AM
The position file formats ".clocs" and "_pos.txt"? Ist there any difference? elgor Illumina/Solexa 0 06-27-2011 07:55 AM
"Systems biology and administration" & "Genome generation: no engineering allowed" seb567 Bioinformatics 0 05-25-2010 12:19 PM

Reply
 
Thread Tools
Old 09-01-2013, 10:17 PM   #1
flyingoyster
Member
 
Location: NJ

Join Date: Aug 2011
Posts: 10
Default "not_aligment" from the output of HTseq-count

What's the meaning of "not_alignment" from the output of HTseq-count? The manual says it is the "reads in the Sam/Bam file without alignment". Does this mean the reads that are not mapped to the reference? or anything else? Thanks!
flyingoyster is offline   Reply With Quote
Old 09-02-2013, 01:28 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

Yeah, these are just reads that didn't align. Whether the count for this is 0 or not will depend on the aligner that you're using (some don't put unaligned reads in the output SAM/BAM file, others do).
dpryan is offline   Reply With Quote
Old 09-02-2013, 05:13 PM   #3
flyingoyster
Member
 
Location: NJ

Join Date: Aug 2011
Posts: 10
Default

Assuming that the file I used for input of HTseq-count includes the unaligned reads in the output BAM file, But the number of "no_alignment" reads doesn't match the number of reads not mapped to the reference. Any idea about this? Thanks!
flyingoyster is offline   Reply With Quote
Old 09-02-2013, 10:57 PM   #4
Simon Anders
Senior Member
 
Location: Heidelberg, Germany

Join Date: Feb 2010
Posts: 994
Default

it should match. Maybe it's paired-end? (htseq-count counts read pairs, not reads)
Simon Anders is offline   Reply With Quote
Old 09-02-2013, 11:10 PM   #5
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

@flyingoyster: Try using the -o option to see what's going on in your case.
dpryan is offline   Reply With Quote
Old 09-03-2013, 12:21 AM   #6
flyingoyster
Member
 
Location: NJ

Join Date: Aug 2011
Posts: 10
Default

@Simon Anders, Yes, My data is paired-end. If HTseq-count counts read pairs, not reads, that makes sense. Thanks a lot!
@dpryan, I used Galaxy server, so not be able to choose options.

If I want to calculate the number of reads within genes, is it right to substract the number of no_feature reads from the total number of mapped reads? Thanks!
flyingoyster 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 10:17 AM.


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