[BioC] Where is DEXSeq bioconductor package?
fabrice.ciup at gmail.com
Thu Aug 4 19:17:50 CEST 2011
I am curios what is the potential problem we will have if the -s
parameters set wrong? If there is a way to know the -s parameters from
When I run dexseq_count.py on my sam file from bwa, it always has this warning:
UserWarning: Malformed SAM line: RNAME != '*' although flag bit
algnt = SAM_Alignment.from_SAM_line( line )
UserWarning: Read GA2-EAS337_0017_FC:1:8:11202:16260#0 claims to have
an aligned mate which could not be found. (Is the SAM file properly
"which could not be found. (Is the SAM file properly sorted?)" )
but I can found GA2-EAS337_0017_FC:1:8:11202:16260#0 mate in the data,
as show in the attachments a.sam. And the a.sam is sorted by the name
(default by bwa).
In htseq-count, it also have this problem.
htseq-count a.sam Homo_sapiens.GRCh37.63/Homo_sapiens.GRCh37.63.gtf
100000 GFF lines processed.
200000 GFF lines processed.
2000000 GFF lines processed.
2064769 GFF lines processed.
Warning: Malformed SAM line: RNAME != '*' although flag bit &0x0004 set
Warning: Read GA2-EAS337_0017_FC:1:8:11202:16260#0 claims to have an
aligned mate which could not be found. (Is the SAM file properly
4 sam line pairs processed.
Why here report 4 line pairs? In fact, there are 3 sam line pairs in a.sam.
BTW: in your opinions, to do this kind of count using BWA or tophat
mapping, which one is better? (Illumina pair-end human RNA-seq data)
More information about the Bioconductor