bedtools not working with gff from flybase

Hi,

I got a gff file from flybase that i am trying to use with bedtools. However, it is throwing an error:

(0)> bedtools intersect -a my.vcf -b dmel-2L-r5.41.gff
Error: Invalid record in file dmel-2L-r5.41.gff. Record is 
2L      FlyBase chromosome_band -204333 1326937 .       +       .       ID=band-21_chromosome_band;Name=band-21

I downloaded the gff from here

I don't know why this format would be incorrect. From what I know of gff files, it seems correct.

Any help would be appreciated.

Thank you,
Best,

asif

bedtools gff flybase vcf

• 2,432 views • link

written 10.4 years ago by asif.zubair •

Why's your start coordinate negative?

• 0 views • link

written 10.4 years ago by Sam

1 answer

Yes, it turns out there is a problem with flybase's gff file. It seems flybase people fixed that in a newer version.

Thank you for pointing it out!

• 0 views • link

written 10.4 years ago by asif.zubair •

More posts like this

how to get sequences by location when pyfaidx.Fasta(read_long_names=True) creates keys from FastaRe…
written 2.6 years ago by mk

I'm attempting to use pyfaidx to index the dmel r6 genome build so that I can get actual sequences from tuples like `(chromosome, start,end)`. In …
bedtools intersect VCF with BED Invalid record error.
written 3.8 years ago by nicole.zeltser •

I am attempting to intersect a VCF file with coordinates from a BED file. e.g. bedtools intersect -a /path/to/vcf -b /path/to/.bed The command results in …
STAR indexing error: terminate called after throwing an instance of 'std::out_of_range'
written 4.3 years ago by khf0008 •

Hi all, This is my first time posting on Biostars, and first time doing any sort of bioinformatics project, so please be gentle! I have …
Error in gff file
written 4.8 years ago by sherafzalk769 •

I am trying to use gff with bedtools intersect to get the reads count to the gene. However, it is throwing an error: ``` Error: …
Bedtools intersect error in record file
written 6.2 years ago by jamie.pike

I keep getting an error with bedtools intersect when comparing a .bed file to a .gff, e.g. Error: Invalid record in file /home/u1983390/Fusarium_data/MIMPS/MIMP_Searches/IntersectBED/GFFS_SameNameFormula/AGND.SNF.gff. Record is …
How to retrieve VCF of C. elegans and D. melanogaster
written 7.0 years ago by Fadel •

Hi, As simple as this, I couldn't know how to find VCF variation files of C. elegans nor D. melanogaster. Since dbSNP does only host …
Problem with .vcf annovar
written 7.3 years ago by zizigolu

Hi, I faced a problem I am trying to annotate my `.vcf` files by `annovar` to get `maf` as by maf I am able to …
Using xdformat on fasta from REPBASE returns erorr.
written 10.5 years ago by asif.zubair •

I am trying to use `xdformat` to make `xdf` files from `fasta` files downloaded from REPBASE - http://www.girinst.org/server/archive/RepBase14.04. I concatenated all the references: cat RepBase14.04.fasta/*.ref …
Bedtools Intersect Error
written 10.8 years ago by dally

I am trying to intersect two bed files using the intersectBed command of bedtools. But I am generating an error when I haven't run into …