I am having trouble running a .sh file. The code 'x=${file_name:0:$z-11}' is giving me a bad substitution error. However when I run in BASH it works. Thing is when this goes to production the .sh will not be running in BASH. Is there a way to substring a string not in BASH or a way to invoke... (2 Replies)
Hi all-
I'm trying to search through some .gz log files to verify certain feeds have passed through our app.
I have a small script that I wrote in hopes that I could automate the checking but haven't been able to get the zgrep to work. When I copy it to the command line directly it works... (2 Replies)
Ubuntu 9.10 is my linux distro
Based on forums they say that the problem is with environment .
here is my case:
login as user, then sudo -s
using this command: s3cmd put file s3://bucket >>worked!
now here is the simple script intended for testing:
#! /bin/bash
env >/tmp/cronjob.log... (1 Reply)
Hey guys. Hopefully this is an easy one but having reference similar problems on the web I still can't fix it.
I am doing a recursive find and replace from a script. Of course I could just run the damn thing from the command line but it's bugging me now and want to get it working.
grep -rl... (4 Replies)
I'm reading about debugging aids in bash and have come across the set command. It says in my little book that an addition to typing
set
you can also use them "on the command line when running a script..." and it lists this in a small table:
set -o option Command Line... (5 Replies)
Hi All,
This command works when I type it on but when I run the batch file it doesn't..any ideas why?
attrib.exe * | find /c /v "" >filecount.txt (1 Reply)
I am working with a sh script on a solaris 9 zone (sol 10 host) that grabs information to build the configuration command line. the variables Build64, SSLopt, CONFIGopt, and CC are populated in the script. the script includes
CC=`which gcc`
CONFIGopt=' --prefix=/ --exec-prefix=/usr... (8 Replies)
OSX 10.9
I am building a script that evaluates the difference between 2 files. Here is a command that does not work transparently.
Running this command in Terminal yields great results; however when I put that line in a .sh script, I get the errors shown below. Am I doing something silly?
... (1 Reply)
Oracle Linux 6.
Trying to set up a simple monitoring of memory usage. This command does exactly what I want at the command line:
echo $(date +%Y-%m-%d" "%H:%M:%S) $(grep PageTables /proc/meminfo) >> /home/oracle/meminfo.logBut when I put it in my crontab:
* * * * * echo $(date +%Y-%m-%d"... (2 Replies)
The below command works in the terminal interactively but not as part of a bash script. I though maybe I needed to escape the "$dir" so it isn't interpreted literally, but that's not it. Thank you :).
interactively in terminal
dir=/path/to
new=$(ls "$dir"/*.csv -tr | tail -n 1) && echo... (6 Replies)
Discussion started by: cmccabe
6 Replies
LEARN ABOUT DEBIAN
vcf-annotate
VCF-ANNOTATE(1) User Commands VCF-ANNOTATE(1)NAME
vcf-annotate - annotate VCF file, add filters or custom annotations
SYNOPSIS
cat in.vcf | vcf-annotate [OPTIONS] > out.vcf
DESCRIPTION
About: Annotates VCF file, adding filters or custom annotations. Requires tabix indexed file with annotations.
Currently annotates only the INFO column, but it will be extended on demand.
OPTIONS -a, --annotations <file.gz>
The tabix indexed file with the annotations: CHR FROM[ TO][ VALUE]+.
-c, --columns <list>
The list of columns in the annotation file, e.g. CHROM,FROM,TO,-,INFO/STR,INFO/GN. The dash in this example indicates that the third
column should be ignored. If TO is not present, it is assumed that TO equals to FROM.
-d, --description <file|string>
Header annotation, e.g. key=INFO,ID=HM2,Number=0,Type=Flag,Description='HapMap2 membership'. The descriptions can be read from a
file, one annotation per line.
-f, --filter <list>
Apply filters, list is in the format flt1=value/flt2/flt3=value/etc.
-h, -?, --help
This help message.
Filters:
+ Apply all filters with default values (can be overridden, see the example below).
-X Exclude the filter X
1, StrandBias
FLOAT Min P-value for strand bias (given PV4) [0.0001]
2, BaseQualBias
FLOAT Min P-value for baseQ bias [1e-100]
3, MapQualBias
FLOAT Min P-value for mapQ bias [0]
4, EndDistBias
FLOAT Min P-value for end distance bias [0.0001]
a, MinAB
INT Minimum number of alternate bases [2]
c, SnpCluster
INT1,INT2 Filters clusters of 'INT1' or more SNPs within a run of 'INT2' bases []
D, MaxDP
INT Maximum read depth [10000000]
d, MinDP
INT Minimum read depth [2]
q, MinMQ
INT Minimum RMS mapping quality for SNPs [10]
Q, Qual
INT Minimum value of the QUAL field [10]
r, RefN
Reference base is N []
W, GapWin
INT Window size for filtering adjacent gaps [10]
w, SnpGap
INT SNP within INT bp around a gap to be filtered [10]
Example:
zcat in.vcf.gz | vcf-annotate -a annotations.gz -d descriptions.txt | bgzip -c >out.vcf.gz zcat in.vcf.gz | vcf-annotate -f
+/-a/c=3,10/q=3/d=5/-D -a annotations.gz -d descriptions.txt | bgzip -c >out.vcf.gz
Where descriptions.txt contains:
key=INFO,ID=GN,Number=1,Type=String,Description='Gene Name' key=INFO,ID=STR,Number=1,Type=Integer,Description='Strand'
vcf-annotate 0.1.5 July 2011 VCF-ANNOTATE(1)