bedtools probelm of syntax
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
1 answer
You need to supply a value for the -f argument. For example, if you require 50% overlap, provide -f 0.05.
• 0 views
•
link
updated
by
Ram
•
written
by
Aaronquinlan
•
written
by
Aaronquinlan
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).
• 0 views
•
link
updated
by
Ram
•
written
by
Aaronquinlan
•
written
by
Aaronquinlan
Log in to answer this question.
More posts like this
-
bamCoverage with DNAse-seq shift
written by JMR
-
How to get sample names and genotype for SNP in multi-sample VCF file with special occurence
written by User000
-
Running a code in batch mode on a list of files
written by zizigolu
-
Tutorial: Calculation of TSS Scores for ATAC-seq Data
written by Sael •
-
Merge series and recognize end of it - AWK
written by Korsocius
-
awk command in shell script
written by mittu1602
-
Reformat a table with awk
written by viv3kanand •
-
ERROR: illegal character '.' when running bedtools closest command
written by Matus •
-
coverageBed with -g option ?
written by bioguy24
-
Sorting to group
written by filipzembol
