-
Notifications
You must be signed in to change notification settings - Fork 3
/
config_file.sh
executable file
·24 lines (18 loc) · 2.12 KB
/
config_file.sh
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
##########################################################################################
run_he_script="run_hyper_editing.sh" # if needed insert the proper path before the script name
bwa_aln_soft="bwa" # if needed insert the proper path before the tool name
bwa_mem_soft="bwa" # if needed insert the proper path before the tool name
SamToFastq_soft="bam2fastx" # if needed insert the proper path before the tool name
SamTools_soft="samtools" # if needed insert the proper path before the tool name
source_file="file_list" # path+full_name of the input files to run: fastq_file_path+name /TAB/ out_prefix (if the input files are of paired-end reads, each of the paired files should appear in separate line).
genome_bwa_ind="genome_prefix" # path+prefix of the index genome expected 5 files like: prefix.amb, prefix.ann, prefix.bwt, prefix.pac, prefix.sa
genome_trans_bwa_ind="transformed_genome_prefix" # path+prefix of the index transformed genome: for each of the 6 transformed (a2c a2g a2t g2c t2c t2g)=>tt: 5 files: prefix.tt.amb, prefix.tt.ann, prefix.tt.bwt, prefix.tt.pac, prefix.tt.sa + 1 fasta file prefix.tt.fa => tot 36 files
genome_fasta="genome_prefix.fa" # path+full_name of the fasta file of the original genome
Q="33" # PHRED score offset of base quality (33 or 64)
PE="1" # 0 => single end, 1 => paired end
GAP=500000 # gap max size between the pairs
dir_pre="output_folder_prefix" # path+full_name of the output directory
bwa_run="1" # 0 => pipeline will not run BWA mapping if analyseMM file exists, 1 => first run BWA mapping of the source files
ue_detect_args="0.05 0.6 30 0.6 0.1 0.8 0.2" # args meaning: -Min of edit sites at Ultra-Edit read -Min fraction of edit sites/mm sites -Min sequence quality for counting editing event -Max fraction of same letter in cluster -Min of cluster length -Max initiate index of cluster -Min ending index of cluster
################################################################################################
$run_he_script $genome_bwa_ind $genome_trans_bwa_ind $genome_fasta $Q $PE $GAP $dir_pre $bwa_run $ue_detect_args $source_file $bwa_aln_soft $bwa_mem_soft $SamToFastq_soft $SamTools_soft