Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • intersectBed gives empty file

    Hi All,

    Stumbled on a strange error when running intersectBed (BEDtools) on two samples that I know are partially overlapping. Here are the first five rows of the two files:

    File1.gtf, first 5 rows:
    Code:
    chr1	bed2gff	exon	4541	4636	0.0	*	.	 gene_id "chr1:4540-4636_1"; transcript_id "chr1:4540-4636_1"; exon_number "1";
    chr1	bed2gff	exon	4640	4692	15.0	-	.	 gene_id "chr1:4639-6522_1"; transcript_id "chr1:4639-6522_1"; exon_number "1";
    chr1	bed2gff	exon	4793	4901	15.0	-	.	 gene_id "chr1:4639-6522_1"; transcript_id "chr1:4639-6522_1"; exon_number "2";
    chr1	bed2gff	exon	6470	6522	15.0	-	.	 gene_id "chr1:4639-6522_1"; transcript_id "chr1:4639-6522_1"; exon_number "3";
    chr1	bed2gff	exon	4640	4692	8.0	-	.	 gene_id "chr1:4639-6522_2"; transcript_id "chr1:4639-6522_2"; exon_number "1";
    File2.gtf, first 5 rows:
    Code:
    chr1	bed2gff	exon	4541	4636	0.0	*	.	 gene_id "chr1:4540-4636"; transcript_id "chr1:4540-4636"; exon_number "1";
    chr1	bed2gff	exon	4640	4692	0.0	-	.	 gene_id "chr1:4639-6522"; transcript_id "chr1:4639-6522"; exon_number "1";
    chr1	bed2gff	exon	4793	4901	0.0	-	.	 gene_id "chr1:4639-6522"; transcript_id "chr1:4639-6522"; exon_number "2";
    chr1	bed2gff	exon	5610	5810	0.0	-	.	 gene_id "chr1:4639-6522"; transcript_id "chr1:4639-6522"; exon_number "3";
    chr1	bed2gff	exon	5864	5890	0.0	-	.	 gene_id "chr1:4639-6522"; transcript_id "chr1:4639-6522"; exon_number "4";
    I ran
    Code:
    intersectBed -a File2.gtf -b File1.gtf -wa -wb -f 1.0 -s -r  > overlap_newSegments_internal
    and end up with an empty result file.

    Any idea why this might be? Is intersectBed sensitive to the order in the files?

    Thanks,
    Boel
    Last edited by Boel; 12-29-2010, 08:28 AM.

  • #2
    What version are you using? I copied your examples and it worked with v2.10.0:
    Code:
    chr1	bed2gff	exon	4541	4636	0.0	*	.	 gene_id "chr1:4540-4636"; transcript_id "chr1:4540-4636"; exon_number "1";	chr1	bed2gff	exon	4541	4636	0.0	*	.	 gene_id "chr1:4540-4636_1"; transcript_id "chr1:4540-4636_1"; exon_number "1";
    chr1	bed2gff	exon	4640	4692	0.0	-	.	 gene_id "chr1:4639-6522"; transcript_id "chr1:4639-6522"; exon_number "1";	chr1	bed2gff	exon	4640	4692	15.0	-	.	 gene_id "chr1:4639-6522_1"; transcript_id "chr1:4639-6522_1"; exon_number "1";
    chr1	bed2gff	exon	4640	4692	0.0	-	.	 gene_id "chr1:4639-6522"; transcript_id "chr1:4639-6522"; exon_number "1";	chr1	bed2gff	exon	4640	4692	8.0	-	.	 gene_id "chr1:4639-6522_2"; transcript_id "chr1:4639-6522_2"; exon_number "1";
    chr1	bed2gff	exon	4793	4901	0.0	-	.	 gene_id "chr1:4639-6522"; transcript_id "chr1:4639-6522"; exon_number "2";	chr1	bed2gff	exon	4793	4901	15.0	-	.	 gene_id "chr1:4639-6522_1"; transcript_id "chr1:4639-6522_1"; exon_number "2";

    Comment


    • #3
      Yes there is nothing wrong with my files, and I was using an up to data version of BEDtools. But! The second file was HUGE, and the command worked when I changed the order of the files. So problems solved now.

      Comment

      Latest Articles

      Collapse

      • seqadmin
        The Impact of AI in Genomic Medicine
        by seqadmin



        Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
        02-26-2024, 02:07 PM
      • seqadmin
        Multiomics Techniques Advancing Disease Research
        by seqadmin


        New and advanced multiomics tools and technologies have opened new avenues of research and markedly enhanced various disciplines such as disease research and precision medicine1. The practice of merging diverse data from various ‘omes increasingly provides a more holistic understanding of biological systems. As Maddison Masaeli, Co-Founder and CEO at Deepcell, aptly noted, “You can't explain biology in its complex form with one modality.”

        A major leap in the field has
        ...
        02-08-2024, 06:33 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, Yesterday, 06:12 AM
      0 responses
      19 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 02-23-2024, 04:11 PM
      0 responses
      67 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 02-21-2024, 08:52 AM
      0 responses
      74 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 02-20-2024, 08:57 AM
      0 responses
      66 views
      0 likes
      Last Post seqadmin  
      Working...
      X