Question: Problem with bedtools intersect with two gtfs
0
gravatar for newbie
14 days ago by
newbie90
newbie90 wrote:

I'm trying to use bedtools intersect on gencode annotation gtf and another gtf which is output from stringtie merge option.

bedtools intersect -v -b gencode.v27.annotation.gtf -a stringtiemerge.gtf > final.gtf

This is how gencode gtf looks:

enter image description here

And this is how the stringtie merged gtf looks:

enter image description here

But when I used bedtools intersect, I see there is warning. I can understand that it is because of the contigs which are not found in gencode gtf. So, the below warning, should I consider that or no need to care?

bedtools intersect -v -b gencode.v27.annotation.gtf -a stringtiemerge.gtf > final.gtf

***** WARNING: File stringtiemerge.gtf has inconsistent naming convention for record:
GL000008.2  StringTie   transcript  1577    8168    .   +   .   transcript_id "MSTRG.4.1"; gene_id "MSTRG.4"; xloc "XLOC_000001"; class_code "u"; tss_id "TSS1";

***** WARNING: File stringtiemerge.gtf has inconsistent naming convention for record:
GL000008.2  StringTie   transcript  1577    8168    .   +   .   transcript_id "MSTRG.4.1"; gene_id "MSTRG.4"; xloc "XLOC_000001"; class_code "u"; tss_id "TSS1";
rna-seq intersect bedtools gtf • 117 views
ADD COMMENTlink modified 14 days ago • written 14 days ago by newbie90
1

perhaps because the order of the attributes is unexpected? in the file you have gene_id followed by transcript_id, the warning shows it flipped. It shouldn't matter but it might.

(BTW including the a few lines from the data instead of images would make more sense)

ADD REPLYlink written 14 days ago by Istvan Albert ♦♦ 85k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1421 users visited in the last hour