HTSeq-Count

"Julia [guest]" <[email protected]> Fri, 5 Sep 2014 07:41:05 -0700 (PDT)
Newsgroups gmane.science.biology.informatics.conductor
Message-ID <[email protected]>
Hi all,
I am new to the field of seq and performed a RIP-Seq experiment using HTSeq count as counter.
I get now the following (using union, but doesn´t look better for interesection_strict):
__no_feature	1503377
__ambiguous	490772
__too_low_aQual	0
__not_aligned	0
__alignment_not_unique	5277314
	
When I sum up counts for all genes, I get 3227845.

The number for __no_feature, __ambiguous, __alignment_not_unique look very high.

Does somebody have an idea for that?

(Additional info: We did random priming and mapped with STAR and masked rRNA loci)

Best wishes
Julia

 -- output of sessionInfo(): 

.

--
Sent via the guest posting facility at bioconductor.org.

_______________________________________________
Bioconductor mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/bioconductor
Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor