forked from pinetree1/crispr-dav
-
Notifications
You must be signed in to change notification settings - Fork 0
/
conf.txt.template
64 lines (45 loc) · 1.79 KB
/
conf.txt.template
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
### Case sensitive for all section names, parameters, and values
## Add a genome section to reflect the genome of interest.
## If there are multiple genomes, each genome must have its own section.
## When amplicon sequence is used as reference, this section is ignored.
#
[genomex]
ref_fasta = genome/genomex.fa
bwa_idx = genome/genomex.fa
refGene = genome/refgenex.txt
## Required tools
[app]
# The paths of abra jar file.
abra = /path/to/ABRA/abra-0.97-SNAPSHOT-jar-with-dependencies.jar
# Path to prinseq-lite.pl. Make sure prinseq path is executable.
prinseq = /path/to/prinseq-lite-0.20.4/prinseq-lite.pl
# bwa must be added to PATH
# Path of samtools. By default use samtools in PATH.
samtools = /path/to/samtools-1.3.1/samtools
# Path of flash. By default use flash in PATH.
flash = /path/to/flash2
# bedtools (e.g. v2.25) must support: bedtools intersect -F
# Default is bedtools in PATH.
bedtools = /path/to/bedtools2.25/bin/bedtools
# Java must be 1.7 or later. Default is java in PATH
java = /path/to/jdk1.7.0_51/bin/java
# Path of pysamstats executable. By default use pysamstats in PATH.
pysamstats = /path/to/bin/pysamstats
# R must have ggplot2, reshape2, and naturalsort packages
# By default use Rscript in PATH.
rscript = /path/to/R-3.2.1/bin/Rscript
## Filtering parameters
[prinseq]
# remove reads with quality score mean below this value (default: 30)
min_qual_mean = 30
# remove reads with length less than this value (default: 50)
min_len = 50
# remove reads with percentage of Ns over this value (default: 3).
ns_max_p = 3
[other]
# realign for large indel detection with ABRA: Y(default) or N
realign_flag = Y
# require minimum BWA mapping quality score (default: 20)
min_mapq = 20
# number of bases on each side of sgRNA to view base changes (default: 40)
wing_length = 12