Question: FeatureCounts on HG38
1
gravatar for elb
2.1 years ago by
elb180
Torino
elb180 wrote:

Hi guys, I'm trying to move from HG19 to HG38. While running FeatureCounts using Refseq*.gtf I see some differences between HG19 and HG38. This I suppose is expected but I only would like to know if I'm doing wrong since I'm a newbie in the RNA-Seq HG38 analysis. Here some outputs of FeatureCounts for one sample (single end) as example:

Sample1.

HG19:
Features : 460394
Meta-features : 24103
Chromosomes/contigs : 49
Total reads : 24944206.
Successfully assigned reads : 11116093 (44.6%)

The same sample with HG38:
Features : 536042
Meta-features : 50502
Chromosomes/contigs : 256
Total reads : 27605790
Successfully assigned reads : 6766947 (24.5%)

The alignment was performed using STAR and in both cases the Uniquely mapped reads % was around 93% or in any case it was comparable between the HG19 and HG38. Moreover, FeatureCounts was used as follows:

Fcounts <- featureCounts(files=fileName, 
                           annot.ext="/.../hg38.refseq.gtf", 
                           isGTFAnnotationFile=TRUE, GTF.attrType="gene_id", 
                           GTF.featureType="exon", nthreads=4, minMQS=3, isPairedEnd=FALSE)

Can anyone give me some feedback?

Thank you in advance

hg19 rna-seq star refseq • 1.1k views
ADD COMMENTlink modified 2.1 years ago by h.mon30k • written 2.1 years ago by elb180
1

I changed your post from Job to Question.

ADD REPLYlink modified 2.1 years ago • written 2.1 years ago by Medhat8.7k

At least if you run featureCounts on the command line you get a file with summary metrics, if R produces that as well then have a look at it. If R doesn't produce that file then run featureCounts directly rather than via the R wrapper.

ADD REPLYlink written 2.1 years ago by Devon Ryan96k
2

the stats are in Fcounts$stat

ADD REPLYlink written 2.1 years ago by WouterDeCoster44k

How the number of reads is different, if it is the same sample?

Total reads : 24944206

And

Total reads : 27605790

ADD REPLYlink written 2.1 years ago by h.mon30k

Successfully assigned reads is different

ADD REPLYlink written 2.1 years ago by elb180
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 877 users visited in the last hour