Add all non overlapping genes/transcripts from one gtf to another with bedtools intersect
If trying to add all non-overlapping transcripts from one gtf to another gtf would would the appropriate bedtools intersect command be.
I am starting with:
intersectBed –wa –a augustus.gtf –b genemark.gtf > mynew.gtf
But I think this will only write the original entry in A for each overlap but not the non overlapping parts?
• 1,554 views
•
link
written
by
robert.murphy •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Bedtools give different output
written by anikcropscience •
-
Determine the significance of peak overlaps (bedtools)
written by kstangline •
-
bedtools intersect and find genes
written by hosin •
-
Combining two databases in gtf file, bedtools gives wierd result
written by marina.v.yurieva
-
Bedtools intersect -u or -wa causing issues
written by alex.v.nesta •
-
Beedtools intersect: more peaks intersected than intersected portions of those peaks?
written by salamandra
-
bedtools closest - result for only the first query
written by yjyj •
-
Bedtools intersect question
written by xiaoyonf •
-
Get completely overlapping features
written by iraun
-
Fraction of bases covered by each element in bed
written by avinash.sahu •
Option
-vis what you want.So use
-vto get non overlapping then later on join the two together?That could be an option, yes.
Is there a better way you would suggest?