Question: HTseq or featureCounts without gff file
0
gravatar for arta
2.0 years ago by
arta540
Sweden
arta540 wrote:

Hi,

I am doing RNA-seq analysis for allele-specific gene expression. To do that, two mouse strains were hybrid. I do have reference fasta files for both but do no not have gff file for one strain.

So far, I have done the mapping and alignment using STAR and I got Aligned.sortedByCoord.out.bam. Now I would like to compare read counts that are mapped same exons or genes to find allele-specific gene expression. But I am not sure which tool to use. I know both HTseq and featureCounts require gff file. Any ideas ??

sequencing snp rna-seq next-gen • 1.2k views
ADD COMMENTlink modified 2.0 years ago by Devon Ryan90k • written 2.0 years ago by arta540
2
gravatar for Devon Ryan
2.0 years ago by
Devon Ryan90k
Freiburg, Germany
Devon Ryan90k wrote:

You should be using one of the many available ASE programs rather than trying to roll your own solution.

Anyway, this is what lapels and suspenders are for. You then only need an annotation for one parental line.

ADD COMMENTlink written 2.0 years ago by Devon Ryan90k

Thanks Devon, I will look them.

ADD REPLYlink written 2.0 years ago by arta540
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: 934 users visited in the last hour