Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
How can i get the number of read aligned to one scaffold youngjune29bhak Bioinformatics 0 08-03-2015 11:13 PM
bowtie number of mismatches and multiple aligned reads carolW Bioinformatics 17 12-23-2014 05:03 PM
How to count aligned RNA-seq reads after sequenced and aligned by Illumina? IceWater Illumina/Solexa 5 04-05-2012 09:18 AM
TopHat/Bowtie - number of reads aligned mgibson Bioinformatics 7 10-22-2011 08:04 PM
the number of aligned reads by maq and bwa totalnew Bioinformatics 3 12-14-2009 11:18 AM

Thread Tools
Old 10-21-2018, 11:06 PM   #1
Junior Member
Location: Jeddah

Join Date: Oct 2018
Posts: 1
Default Htseq_processing double number of aligned reads

Hi All,

I am trying to generate count data from the bam files using Htseq. Our data was produced using Illumina stranded library and sequenced using Hiseq 3000. After mapping using Hisat2, I know that there are about 37 million aligned read pairs. However, when I am generating count data using HTseq I can see that it has processed ~ 78 million reads and the total counts from all features is coming around 64 million which I feel is not correct. If we have 37 million aligned reads how can we get double number of count. I have used -s reverse option as the library was prepared using Illumina stranded protocol. I am not able to guess where things are going wrong. It would be great if anyone can give me a probable reason for this.

Best regards, Amit
amit123 is offline   Reply With Quote

expression analysis, hisat2, htseq count, rnaseq alignment

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 01:17 PM.

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