Hi,
Apparently I need to use a gtf or gff3 file for annotation with easyRNAseq now when processing hg19 bamfiles now - as method biomaRt uses GRCh39 and I get an error. So far so bad.
However, I seem not to be able to download a gtf or gff3 file for hg19 that easyRNAseq would accept:
> DGE<-simpleRNASeq(bamFiles=getBamFileList(bamfiles), param=param, nnodes=4, verbose=T, override=T)
...
==========================
Processing the annotation
==========================
Validating the annotation source
Read 993 records
Error in .validate(obj, verbose = verbose) :
The provided gff3 contains no annotation of type 'mRNA' and/or 'exon' in the first 1000 lines.
I've tried files the following annotation files but none seems to work:
ref_GRCh37.p13_top_level.gff3.gz
ref_GRCh37.p13_scaffolds.gff3.gz
gencode.v19.annotation.gtf
gencode.v19.annotation.gff3
Homo_sapiens.GRCh37.70.gtf
will try the Illumina iGenome next. What am I missing here?!?
Apparently I need to use a gtf or gff3 file for annotation with easyRNAseq now when processing hg19 bamfiles now - as method biomaRt uses GRCh39 and I get an error. So far so bad.
However, I seem not to be able to download a gtf or gff3 file for hg19 that easyRNAseq would accept:
> DGE<-simpleRNASeq(bamFiles=getBamFileList(bamfiles), param=param, nnodes=4, verbose=T, override=T)
...
==========================
Processing the annotation
==========================
Validating the annotation source
Read 993 records
Error in .validate(obj, verbose = verbose) :
The provided gff3 contains no annotation of type 'mRNA' and/or 'exon' in the first 1000 lines.
I've tried files the following annotation files but none seems to work:
ref_GRCh37.p13_top_level.gff3.gz
ref_GRCh37.p13_scaffolds.gff3.gz
gencode.v19.annotation.gtf
gencode.v19.annotation.gff3
Homo_sapiens.GRCh37.70.gtf
will try the Illumina iGenome next. What am I missing here?!?
Comment