-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathconfig_file.txt
57 lines (57 loc) · 2.33 KB
/
config_file.txt
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
sample_name=sample
sample_fastq_dir=/mnt/hadoop-cnvrf-public/data/sample
control_name=control
control_fastq_dir=/mnt/hadoop-cnvrf-public/data/control
final_results_folder=/mnt/hadoop-cnvrf-public/sample_results
# Copy of genome indices
local_hg19_index=/mnt/hg19/hg19_index.tar.gz
# tso_reference.bed must have the fourth column (name column)
ref_pileup=/mnt/hadoop-cnvrf-public/tso_reference_aug2018.bed
exons_bed=/mnt/hadoop-cnvrf-public/subset.bed
# Number of references to use
num_in_sample_references=4
# Min number reference that need to contain a variant for it to be called
min_num_variant_references=2
# Ordered genes (comma delimited)
ordered_genes=FANCA:GJB6:FANCC
training_data=/mnt/hadoop-cnvrf-public/het_del_training_aug2018.txt
username=HadoopCnvRF
# ######################################################
# Advanced users
cnv_home=/mnt/hadoop-cnvrf-public
hdfs_dir=/user/hadoop/cnv
hdfs_hg19_compressed=/user/hadoop/cnv/hg19_index.tar.gz
genome_extention=.fa
hdfs_picard_jar_path=/user/hadoop/cnv/picard.jar
# NOTE: The temp folder you specify below has to be present and accessible
# in slave nodes. /tmp is usually present and writablein slave nodes. If you
# want to use a different folder, or /tmp is not available, you will need to
# create one using a bootstrap script
temp_dir=/tmp
R_libs=/R_libs
R_scripts=/mnt/hadoop-cnvrf-public
# See http://www.acgt.me/blog/2015/3/17/more-madness-with-mapq-scores-aka-why-bioinformaticians-hate-poor-and-incomplete-software-documentation
# 0 = maps to 5 or more locations, 1 = maps to 3-4 locations, 3 = maps to 2 locations. We are using 3 as the filter because we
# want reads that map to 2 locations. We are using the samtools view -q filter, (-q INT only include reads with mapping quality >= INT [0])
mapping_quality_filter=3
min_control_coverage=20
debug_mode=0
# NOTE: There is a bug in the prog that needs to be fixed. If you specify
# only one hadoop node, and user has just one pair of fastq files, you end
# up with just one BAM and samtools merge command fails
number_partitions=2
SubSampleNumber=2000
r1_signature=R1
r2_signature=R2
fastq_extension=.fastq
window_length=60
het_del_min_filter=0.25
het_del_max_filter=0.75
gain_min_filter=1.2
het_lengh_threshold=180
hom_lengh_threshold=180
gain_lengh_threshold=800
number_threads=4
num_default_reducers=40
mergebam_timeout=21600000
inputfile_timeout=21600000