Question: bedtools probelm of syntax
0
gravatar for Korsocius
4.2 years ago by
Korsocius110
Korsocius110 wrote:

Hi I would like to use this syntax, but bedtools has a problem with -f function, do you know where could be a problem pls? ("$z" is a bam file.....)

 bedtools map -a /home/filip/Desktop/binary_NIFTY.bed -b $z -c 10,10 -o count,concat -f | awk -v OFS="\t" 'BEGIN {print "CHR\tSTART\tSTOP.............................

Thank you for help.
Filip

dna map bedtools • 1.1k views
ADD COMMENTlink modified 4.2 years ago • written 4.2 years ago by Korsocius110
3
gravatar for Aaronquinlan
4.2 years ago by
Aaronquinlan10k
United States
Aaronquinlan10k wrote:

You need to supply a value for the -f argument.  For example, if you require 50% overlap, provide "-f 0.05".

ADD COMMENTlink written 4.2 years ago by Aaronquinlan10k

But there is a deafult value or not?  (for one base 1bp)

ADD REPLYlink written 4.2 years ago by Korsocius110
1

Yep, default (without providing -f) will detect one base pair.  If you use, -f you are trying to override that default, and as such, must provide a value to use as the override.

    -f    Minimum overlap required as a fraction of A.
        - Default is 1E-9 (i.e., 1bp).
ADD REPLYlink written 4.2 years ago by Aaronquinlan10k

It works......

ADD REPLYlink written 4.2 years ago by Korsocius110
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1348 users visited in the last hour