This is a test version of Biostars. For the public version, visit https://www.biostars.org.
BedTools can't find "hg19.genome"

Hello,

I am running the following command:

bedtools random -n 100000 -l 118 -g hg19.genome >random_h19.bed

and get this error:

Error: The requested genome file (hg19.genome) could not be opened. Exiting!

Any ides how to specify the path to "hg19.genome" and where to find it? (We have an old version of bedtools installed on a large cluster).

Thanks!

chip-seq
• 6,105 views • link
updated by GenoMax GenoMax • written by biostart biostart

1 answer

It is a file with the chromosome sizes: https://github.com/arq5x/bedtools/blob/master/genomes/human.hg19.genome

• 0 views • link
written by GenoMax GenoMax

Ok, I downloaded the genome file from the UCSC web site and manually put it in my folder, now it works

• 0 views • link
written by biostart biostart

Yes, above comment by @biostart works.

• 0 views • link
written by cheyujlee • cheyujlee

Log in to answer this question.

More posts like this
  • chromosome mm10 size with genomecoveragebed
    written by Jerry • Jerry
    Hi, I would to convert my bed file in bedgraph using bedtools. I use this command: genomeCoverageBed -bg -i file.bed -g mm10_canon_chr.sizes > file.bg but …
  • How to calculate enrichment P-value?
    written by biostart biostart
    <p>Hello,</p> <p>Quick question: How to calculate the P value for the enrichment of my dataset in a certain feature?</p> <p>I have calculated (using bedtools), that …
  • [MACS] IndexError: list index out of range
    written by bio_zhangxl • bio_zhangxl
    <p>histone modify data mapping to reference genome</p> <p><em><strong>.sra--&gt;fastq</strong></em></p> <pre> fastq-dump -I --split-3 SRR948816.sra</pre> <p><em><strong>fastq--bowtie2--sam</strong></em></p> <pre> bowtie2 -p 8 -x genome SRR948816.fastq &gt; SRR948816.sam bowtie2 -p …
  • Galaxy "Stitch Maf Blocks" Produced Null Bases
    written by lam • lam
    <p>problem: "stitch MAF blocks" produced NULL bases </p> <p>steps: From <a href='http://genome.ucsc.edu/cgi-bin/das/hg19/dna?segment=chrY:23310464,23331995'>http://genome.ucsc.edu/cgi-bin/das/hg19/dna?segment=chrY:23310464,23331995</a></p> <p>we got the following bases:</p> <pre><code>agtttacccaactcacaggtatttgaggatacaaacccatgactgggcttggctttaaaaattcctatctgaaatttctagt.......... </code></pre> <p>We want to do multiple-species alignment …
  • Any Quick Sort Method For Large Bed File As 20G In Size?
    written by Xianjun Xianjun
    <p>We are often suggested to sort the input bed file by "sort -k1,1 -k2,2n" in order to invokes a memory-efficient algorithm designed for large files, …
  • Tophat Reporting Output Trakcs Failed!!
    written by allqwdd • allqwdd
    <p>2012-10-17 01:01:25] Reporting output tracks [FAILED] </p> <h2><strong><em>Error: bammerge failed to open BAM file tophatSample/tmp/leftkeptr</em></strong><em>eads.mapped.bam</em></h2> <p>that was the error message Is anyone knows why this …
  • Intersectbed/Coveragebed -Split Purify Exon?
    written by Puriney Puriney
    <p>all.reads.bam file records mapped RNA-seq reads data, including: </p> <ol> <li>exon:exon junction</li> <li>exon body</li> <li>intron body</li> <li>exon:intron junction</li> </ol> <p><strong>Q1:</strong> When calculating RPKM for given …
  • Correlation Of Fpkm And Length Normalized Transcript Mapped Read Count
    written by johannes.helmuth • johannes.helmuth
    <p>Hello,</p> <p>in the process of estimating expression for a 16 human tissue dataset (<a href="http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE30611">"Human Body Map 2.0 GSE30611"</a>) I used different methods to estimate …
  • Vcftools,Vcf-Merge Error
    written by Joanne Lim • Joanne Lim
    <p>I had run vcftools to merge some vcf files with no success. As mentioned, the *.vcf files was pre-processed with bgzip and tabix.</p> <p>Here are …
  • Bed File Bedpe Format
    written by Assa Yeroslaviz Assa Yeroslaviz
    <p>Hi,</p> <p>I'm having trouble with converting the bam file into bed -bedpe using the bedtools.</p> <pre><code>workflow: samtools sort -n mut.bam mut.Namesorted bamTobed -i mut.Namesorted.bam -bedpe …