Correct GFF3 translation for FeatureCounts
0
0
Entering edit mode
rridley3 • 0
@0278a549
Last seen 3 months ago
United States

Hi, I am currently running FeatureCounts following mapping meta-transcriptomic data back to various reference genomes. I have predicted genes using prodigal, which outputs into GFF3 format:

JAACZP010000001.1_1 pyrodigal_v1.1.2 CDS 47 1186 127.9 - 0 ID=1_1;partial=00;[…]

I looked to run featurecounts using parameters -t CDS, -p, and -g ID for this file, and it did not map any reads. However, when I use the NCBI downloaded gtf formatted file, appearing as:

JAACZP010000001.1 Genbank region 1 242615 . + . ID=JAACZP010000001.1:1..242615;Dbxref=taxon[…]

SAM file line looks like: D00468:253:HLLHJBCX2:1:1101:9224:1989 77 0 0 * 0 0 NAGAAGGG[...]

I use the exact same parameters, and this file maps 52% of the reads. Is there an error in the way I am passing this file to featureCounts?

RNASeq FeatureCounts Rna • 99 views
ADD COMMENT

Login before adding your answer.

Traffic: 247 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6