Question: Where to find a valid GTF file for stringtie asssembly?
1
gravatar for lakhujanivijay
21 months ago by
lakhujanivijay5.1k
India
lakhujanivijay5.1k wrote:

Hi all

I am trying to locate a valid GTF file for transcriptome assembly of Cenchrus americanus (Pearl Millet) for the genome information avaialble at this NCBI link here.

Assembly with Stringtie the GFF file downloaded from here does not work and throws following error.

$ stringtie -p 55 -c 10 -G cenchrus_americanus.gff -l A -o A.gtf -C A.ref.transcript.covered.txt A.sorted.bam

Error: could not any valid reference transcripts in cenchrus_americanus.gff (invalid GTF/GFF file?)

It is obvious as the GFF file does not contain information for transcripts, just a single type of field in the 3rd column

$ awk -F "\t" '{print $3}' cenchrus_americanus.gff | sort | uniq

region

There is only this paper which I could find (referenced to the genome provided by NCBI)

There is one other link which I found here and which have a couple of files, however, not sure what all do they mean. Can anybody help me with finding the proper GTF file ? Lastly, I also looked at Ensembl plants database with no luck!

Thanks

Vijay

transcriptome stringtie gtf • 1.1k views
ADD COMMENTlink modified 21 months ago by Damian Kao15k • written 21 months ago by lakhujanivijay5.1k

hey, did you find the GTF file you were looking for? I am also looking for it

ADD REPLYlink written 13 months ago by c_u250
0
gravatar for Damian Kao
21 months ago by
Damian Kao15k
USA
Damian Kao15k wrote:

The .gff file you downloaded doesn't describe gene annotations. It just describes the genome scaffolds/contigs. There doesn't appear to be any gene annotations in the NCBI ftp you posted.

You do not need to provide StringTie with a set of reference annotations. You can perform an assembly just with your transcriptome reads.

ADD COMMENTlink written 21 months ago by Damian Kao15k
1

The .gff file you downloaded doesn't describe gene annotations.

Yes , exactly!

You do not need to provide StringTie with a set of reference annotations.

I am attempting reference-based transcriptome assembly. Any clues for me?

ADD REPLYlink modified 21 months ago • written 21 months ago by lakhujanivijay5.1k
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: 957 users visited in the last hour