This is a test version of Biostars. For the public version, visit https://www.biostars.org.
convert *.bdg file to *.bw file

Hello,

I tried to convert *.bdg file to *.bw file.

I have bedClip software in the folder.

Thanks in advance for any help!

Best,

Yue

administrator@ACB:~/bedtools2-2.26.0$ sh bdg2bw.sh flag1_FE.bdg hg38.len
bdg2bw.sh: 19: bdg2bw.sh: bedClip: not found
/usr/bin/bedtools
sort: cannot read: flag1_FE.bdg.clip: No such file or directory
bdg2bw.sh: 23: bdg2bw.sh: Bad substitution
chip-seq
• 3,942 views • link
written by yueli7 yueli7

1 answer

You need to have software on $PATH, not in the current working folder. You seem to be running into this error a lot. Follow the same procedure I gave you earlier, but with the path to bedClip and bedtools

• 0 views • link
written by Ram Ram

Hello, RamRs

Thank you so much for your great help!

I have all the file in my folder.

Thank you again!

Best,

Yue Li

administrator@ACB-HuangLab-Ubuntu:~/bedtools2-2.26.0/bin$ ./bedClip
bedClip - Remove lines from bed file that refer to off-chromosome locations.
usage:
bedClip [options] input.bed chrom.sizes output.bed
chrom.sizes is a two-column file/URL: <chromosome name=""> <size in="" bases="">
If the assembly <db> is hosted by UCSC, chrom.sizes can be a URL like
http://hgdownload.soe.ucsc.edu/goldenPath/<db>/bigZips/<db>.chrom.sizes
or you may use the script fetchChromSizes to download the chrom.sizes file.
If not hosted by UCSC, a chrom.sizes file can be generated by running
twoBitInfo on the assembly .2bit file.
options:
-truncate  - truncate items that span ends of chrom instead of the
            default of dropping the items
-verbose=2 - set to get list of lines clipped and why
administrator@ACB-HuangLab-Ubuntu:~/bedtools2-2.26.0/bin$ sh bdg2bw.sh
Need 2 parameters! <bedgraph> <chrom info="">
/usr/bin/bedtools
administrator@ACB-HuangLab-Ubuntu:~/bedtools2-2.26.0/bin$ sh bdg2bw.sh    flag1_FE.bdg hg38.len
bdg2bw.sh: 19: bdg2bw.sh: bedClip: not found
/usr/bin/bedtools
bdg2bw.sh: 21: bdg2bw.sh: cannot create flag1_FE.bdg.sort.clip: Permission denied
bdg2bw.sh: 23: bdg2bw.sh: Bad substitution
• 0 views • link
written by yueli7 yueli7

What is the output to:

which -a bedClip
which -a bedtools

If either of those gives you an error, you need to add the location of those binaries to your $PATH. Computers cannot see what you see, so just because the executable in your directory, you shell won't recognize them.

• 0 views • link
written by Ram Ram

Hello, RamRs,

Thank you so much for your great help!

Thank you again!

Best,

Yue 

administrator@ACB-HuangLab-Ubuntu:~/bedtools2-2.26.0/bin$ PATH=$PATH:~/bedtools2-2.26.0/bin sh bdg2bw.sh flag1_FE.bdg chromInfo.txt
• 0 views • link
written by yueli7 yueli7

Also, don't use sh to execute shell scripts unless you're sure the script is POSIX compliant. Use either . or bash. Certain features (such as process substitution) will not work when executed with sh.

• 0 views • link
written by Ram Ram

Hello, RamRs,

Thank you so much for your remaining.

It finally work!

Thank you again!

Best,

Yue li

administrator@Ab-Ubuntu:~/bedtools2-2.26.0/bin$ PATH=$PATH:~/bedtools2-2.26.0/bin ./bdg2bw.sh xaf.txt chromInfo.txt
• 0 views • link
updated by Ram Ram • written by yueli7 yueli7

Log in to answer this question.

More posts like this
  • Nextflow - I/O and index issues
    written by bhumm • bhumm
    This is a crosspost from stack overflow ([Nextflow - I/O and index issues][1]). Based on this thread ([how to handle exact duplicates...][2]) crossposting is allowed …
  • KCCG /ClinSV issue with hg38
    written by Ishak • Ishak
    Dear Biostars I work on KCCG /ClinSV [here][1] to extract the structural wgs with hg38 I run this code as shown clinsv -r all -p …
  • Problem with Unicode while using BioPython on Windows
    written by williamsam • williamsam
    Hi all- I am a graduate student working with an undergraduate student (remotely) on a phylogenetics project. I was showing her how to use the …
  • Build signal track
    written by yueli7 yueli7
    Hello, I followed the post https://github.com/taoliu/MACS/wiki/Build-Signal-Track "Fix the bedGraph and convert them to bigWig files And you will have these bigwig files: flag1_FE_01.bw and flag1_logLR.bw" …
  • Problem installing fcgene using cygwin64.
    written by paulwjeffries • paulwjeffries
    I am trying to install fcgene so that I can convert vcf dosage files to plink files. It seems that I have successfully configured, but …
  • error in star
    written by yueli7 yueli7
    Hello, When I ran STAR, it happens like this: Anyone has suggestion? Thanks in advance for great help! yue administrator@ACB-HuangLab-Ubuntu:~/STAR-2.6.0a/source$ ./STAR --runThreadN 30 --runMode genomeGenerate …
  • GBS pipeline POD Error: does not exist or is unreadable
    written by csjlee08 • csjlee08
    Hello, I am trying to use GBS-POD pipeline, having trouble with the system to find my files although I did specify the file paths in …
  • dependency plyr in R 2.15.1
    written by yueli7 yueli7
    Hi, I try to install qvalue. It seems dependency 'plyr' is not available. My R version is R 2.15.1 ``` > source("http://bioconductor.org/biocLite.R") A new version …
  • SIFT ncbi executables setup incorrectly?
    written by arronslacey arronslacey
    EDIT: solved. whilst I had properly defined the path to the NCBI executables in `SIFT_for_submitting_fasta_seq.csh`, I had not properly set the path in `seqs_chosen_via_median_info.csh` It …