SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   HTSeq - paired reads count (http://seqanswers.com/forums/showthread.php?t=47132)

nofar 10-02-2014 04:02 AM

HTSeq - paired reads count
 
Hello,
I was wondering if there were any changes in the latest version on HTSeq (after Version 0.5.4) that could affect the count of paired reads (other than what's mentioned in the documentation)?
I have a test dataset which I tested a few months ago (before the version update) and I see that the new version gives slightly different results (about 10% of the genes's counts have changed).
Based on the documentation, the only change that could possibly affect the counts is the default for parameter minaqual, but I ran the new version with a=0 (which was the default in the previous version) and there was still a difference in the counts. What else could be affecting this? I haven't changed anything in the dataset or in the parameters/versions of other programs (same version of tophat, same GTF file etc.)
Thanks


All times are GMT -8. The time now is 07:16 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.