SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Breakdancer_max output is empty inou13 Bioinformatics 4 04-25-2012 03:54 PM
GATK IndelRealigner returns empty output file bjfan Bioinformatics 3 06-08-2011 10:43 AM
Bedtools intersectBed with vcf format epigen Bioinformatics 2 03-10-2011 08:56 AM
intersectBed gives empty file Boel Bioinformatics 2 01-04-2011 06:29 AM
Bfast output and "Empty Sequence Dictionary" in .sam output aiden Bioinformatics 1 05-28-2010 06:50 PM

Reply
 
Thread Tools
Old 08-28-2012, 08:23 AM   #1
palc
Senior Scientist
 
Location: Auckland, New Zealand

Join Date: Jun 2012
Posts: 6
Default intersectBed (BEDtools) generating empty output file

I have used the Bedtools command intersectBed to check the overlapping between the two bed files. A is my INDEL file and B is my Reference file. But it is producing empty output file. I thought the problem was that the file B is much larger than file A. But I tried with changing the file order still its not creating any output.

Here is the reference B file (larger):

gff_seqname 0 1395 gene 0 +
gff_seqname 0 1395 exon 0 +
gff_seqname 1397 2498 gene 0 +
gff_seqname 1397 2498 exon 0 +
gff_seqname 2524 3619 gene 0 +

Here is my A file with just 51 INDELS:

CHROMOSOME 174708 174713 -GCCGG:2/6
CHROMOSOME 1078686 1078686 +A:105/112
CHROMOSOME 1229123 1229125 -CT:800/870
CHROMOSOME 1234830 1234830 +AT:134/134
CHROMOSOME 1234833 1234834 -A:134/134

here is my command:

intersectBed -a raw.bed -b sodalis_galaxy.bed -wa -wb >test13.bed
palc is offline   Reply With Quote
Old 08-28-2012, 08:36 AM   #2
palc
Senior Scientist
 
Location: Auckland, New Zealand

Join Date: Jun 2012
Posts: 6
Default

I can see the chromosome in column 1 doesn't match..!!
palc is offline   Reply With Quote
Reply

Tags
bedtools galaxy

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 09:13 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO