gravatar for newbie

4 hours ago by

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";

link

modified 4 hours ago

written
4 hours ago
by

newbie90



Source link