Go Back   SEQanswers > Applications Forums > RNA Sequencing

Similar Threads
Thread Thread Starter Forum Replies Last Post
MiSeq gDNA reads still fail "Kmer content" and "per base seq content" after trimming" ysnapus Illumina/Solexa 4 11-12-2014 07:25 AM
DEXSeq error in estimateDispersions: match.arg(start.method, c("log(y)", "mean")) fpadilla Bioinformatics 14 07-03-2013 02:11 PM
"allele balance ratio" and "quality by depth" in VCF files efoss Bioinformatics 2 10-25-2011 11:13 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

Thread Tools
Old 11-18-2016, 02:59 AM   #1
Junior Member
Location: SPAIN

Join Date: Nov 2016
Posts: 1
Default Too much "no feature" counts??

Iím having some issues analysing RNASeq data on Chipster. I have single-end rna samples (brain) sequenced on a illumina NextSeq platform from rattus norvergicus.
When i run the TopHat (Rattus6.0.83/fr-firststrand/20 hits allowed per read/Ö) i obtain:

Input : 29084606
Mapped : 27577865 (94.8% of input)
of these: 1284892 ( 4.7%) have multiple alignments (17957 have>20)
94.8% overall read mapping rate.

When running the HTSeq (Rattus6.0.83/"reverse"/no paired end/chromosome names in the BAM "1"/ union / exon / gene_id) for counting the reads per gene i found a lot of no_feature reads.

__no_feature 7413760
__ambiguous 111934
__too_low_aQual 0
__not_aligned 0
__alignment_not_unique 4203664

Is it normal?

Then in the differential expression using the DESeq2 i get extremely few or even no significant padj values. Only get 12 signicant padj values even compairing MALES VH and FEMALES VH which seems quite shocking to us.
Iím i doing something wrong?
Thank you very much in advance for the help.
JOM... is offline   Reply With Quote

rnaseq alignment, rnaseq analysis, rnaseq deseq padj, rnaseq htseq

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 05:43 AM.

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