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

Writing saving generated gtf class object #1

Open
BrunoGuillotin opened this issue Sep 14, 2022 · 0 comments
Open

Writing saving generated gtf class object #1

BrunoGuillotin opened this issue Sep 14, 2022 · 0 comments

Comments

@BrunoGuillotin
Copy link

BrunoGuillotin commented Sep 14, 2022

Hi,
Thanks for creating this super simple tool to open and modify gtfs. I used it to add the introns not annotated in my gtf (from Zea mays aka corn). Worked like a charm.

gtf_file <- file.path("Zea_mays.B73_RefGen_v4.47.gtf")
gtf <- read_format(gtf_file)
update gtf with intron coordinates from the exon gaps
and return the updated object
ans <- construct_introns(gtf, update=TRUE)[]default
same as above, but returns only constructed introns
introns <- construct_introns(gtf, update=FALSE)

But how do I save the generated gtf (introns and ans) ?
I would like to use it with bedtools but every ways i tried to save the gtf object give me a file with an excessive amount of column that is not recognise as a gtf.

Thanks in advance for your great help,
B.

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

1 participant