Add all non overlapping genes/transcripts from one gtf to another with bedtools intersect
0
0
Entering edit mode
3.6 years ago

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?

bedtools annotation gtf • 1.0k views
ADD COMMENT
0
Entering edit mode

Option -v is what you want.

ADD REPLY
0
Entering edit mode

So use -v to get non overlapping then later on join the two together?

ADD REPLY
0
Entering edit mode

That could be an option, yes.

ADD REPLY
0
Entering edit mode

Is there a better way you would suggest?

ADD REPLY

Login before adding your answer.

Traffic: 2703 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6