Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Tiny bug - gtf output header #386

Open
SheepwormJM opened this issue Aug 1, 2023 · 2 comments
Open

Tiny bug - gtf output header #386

SheepwormJM opened this issue Aug 1, 2023 · 2 comments

Comments

@SheepwormJM
Copy link

SheepwormJM commented Aug 1, 2023

Using the agat_convert_sp_gff2gtf.pl to convert a gff3 file to a gtf file I found that the output gtf file had a slightly confusing header:

##gtf-version X
# GFF-like GTF i.e. not checked against any GTF specification. Conversion based on GFF input, standardised by AGAT.
##sequence-region hcontortus_chr4_Celeg_TT_arrow_pilon 1 51826579
# Gene gene:HCON_00118380
hcontortus_chrX_Celeg_TT_arrow_pilon    trf     tandem_repeat   1       2516    4529    -       .       gene_id "tandem_repeat-28261"; ID "tandem_repeat-28261";

Clearly the ##sequence-region.... and #Gene gene:HCON... are not needed, but because they're kept it makes the first gene line, a gene on chrX, not chr4, seem wrong at first glance.

General (please complete the following information):

  • AGAT version v1.2.0
  • AGAT installation/use A container I think
  • OS: Ubuntu 18.04

To Reproduce

agat_convert_sp_gff2gtf.pl --gff haemonchus_contortus.PRJEB506.WBPS18.annotations.gff3 --o haemonchus_contortus.PRJEB506.WBPS18.annotations.gtf

Input gff3 file here: https://parasite.wormbase.org/Haemonchus_contortus_prjeb506/Info/Index/

@Juke34
Copy link
Collaborator

Juke34 commented Sep 1, 2023

Thank you for your feedback

@Juke34
Copy link
Collaborator

Juke34 commented Oct 13, 2023

A solution would be to add a bolean option keep_header=true in the config yaml file. Does it sound fine to you? Or we keep cleaning this kind of stuff manually before running AGAT?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants