Entering edit mode
3.6 years ago
robert.murphy
▴
90
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?
Option
-v
is what you want.So use
-v
to get non overlapping then later on join the two together?That could be an option, yes.
Is there a better way you would suggest?