featureCounts results has many psedugenes(ENSG***) together with genes

Dear all,

Do you know how to generate the results to all genes(-t gene already used)?

featureCounts -T 10 -s 2 -p -t gene  -a gencode.v40.primary_assembly.annotation.gtf -o gene_count.txt -g gene_name *bam

enter image description here

featurecounts rna-seq

• 1,468 views • link

updated 4.3 years ago by iraun

• written 4.3 years ago by dew •

1 answer

What is the issue? featureCounts will generate counts for the genes you are providing in your GTF (regardless of being genes or pseudogenes). If you don't want to generate counts for the pseudogenes you can remove them from the GTF or from the output counts file from featureCounts.

• 0 views • link

written 4.3 years ago by iraun

Thanks a lot!

So should the pseudogenes be used in the downstream DEGs finding? Thanks! Best!

• 0 views • link

written 4.3 years ago by dew •

No problem :).

I am afraid that I can't answer to that question. It depends on your hypothesis and experiment. Are you interested in the expression patterns of pseudogenes? Having differentially expressed pseudogenes is something that you are looking for in your experimental design?

• 0 views • link

written 4.3 years ago by iraun

Thanks! Got it. Thank you very much!Best!

• 0 views • link

written 4.3 years ago by dew •

More posts like this

BCFTOOLS VIEW parameter?
written 1.9 years ago by charlie •

I am trying to split my multiple samples vcf files to single one with `Bcftools view`, and `-O z` which can generate the compressed VCF. …
SNP and Genes position matching dataframes in R
written 2.4 years ago by QX •

Hi all, I have 2 datasets that look like this: data1: ![enter image description here][1] data2 ![enter image description here][2] Now I want to loop …
Increasing per base G content in QC sequencing files
written 2.5 years ago by robertsr •

Hello, I conducted some metagenomic sequencing as follows: - Metagenomic sequencing from human stool samples - PCR-free library prep using NEB kit (450bp insert) - …
Annotation GTF/GFF Arabidopsis thaliana
written 2.7 years ago by san96 •

Hello, this is my first time working with Arabidopsis and I am quantifying with featureCounts as follows: featureCounts -p --countReadPairs -t exon -g gene_id -a …
Why did I get more unique genes than unique transcripts w/ STAR -> featureCounts?
written 3.4 years ago by O.rka

Here's my STAR command: ``` STAR --genomeDir human/GRCh38.p13_GENCODE.39_ERCC92 --readFilesIn veba_output/preprocess/31O0P_S1/output/trimmed_1.fastq.gz veba_output/preprocess/31O0P_S1/output/trimmed_2.fastq.gz --outFileNamePrefix star_output/31O0P_S1/intermediate/1__star/ --runThreadN 1 --outReadsUnmapped Fastx --readFilesCommand zcat && samtools sort -@ 1 -T …
featureCounts
written 4.3 years ago by kcarey •

Hey Guys, I have been trying to run this code to generate counts data using featureCounts Subread package on my server: featureCounts -p -T 20 …
strange QQ plots
written 4.4 years ago by dew •

Dear all, May I know if you have ever come across the QQ plot which does not display as close to a line, thank you …
Subread FeatureCounts produces zero percent successfull alignment
written 4.6 years ago by mathavanbioinfo

Dear All, I working on rice genome transcriptome analysis. I have done alignment using the Hisat2 tool process greater than 80% s score. Then I …
Can not plot degree distribution in Cytoscape Network Analyzer on logscale
written 4.8 years ago by bhodai •

I am using Cytoscape 3.9.0. When I tried - Tools > Analyze Network, the following window poped-up ![enter image description here][1] When I clicked on …
Low assignment rate with featureCounts
written 5.0 years ago by gt •

I used STAR to align my reads (brain samples) to human reference genome. Getting good unique mapping rates (~70-90%). ![enter image description here][1] However, when …