forked from brianjohnhaas/indrops
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Early_med_bud.yaml
44 lines (42 loc) · 1.63 KB
/
Early_med_bud.yaml
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
project_name : 'AxoSC4'
project_dir : '/n/scratch2/nl91/indrops/inDrops4_06_04_18'
paths :
bowtie_dir : ''
rsem_dir : ''
python_dir : ''
java_dir : ''
samtools_dir : ''
bowtie_index : '/n/scratch2/nl91/indrops/indrops_index/03_01_18_indrops'
sequencing_runs :
- name : 'Seq3'
version: 'v3'
dir : '/n/scratch2/nl91/indrops/inDrops4_06_04_18'
fastq_path : '20180531_pool_NL5679_S1_{read}_001.fastq.gz'
libraries:
- {library_name: "S_1", library_index: "CTTAATAG"}
- {library_name: "S_2", library_index: "ATAGCCTT"}
- {library_name: "S_3", library_index: "TAAGGCTC"}
- {library_name: "S_4", library_index: "TCGCATAA"}
- {library_name: "S_5", library_index: "TTACCTCC"}
parameters :
umi_quantification_arguments:
m : 10 #Ignore reads with more than M alignments, after filtering on distance from transcript end.
u : 2 #Ignore counts from UMI that should be split among more than U genes.
d : 600 #Maximal distance from transcript end, NOT INCLUDING THE POLYA TAIL
split-ambigs: False #If umi is assigned to m genes, add 1/m to each gene's count (instead of 1)
min_non_polyA: 15 #Require reads to align to this much non-polyA sequence. (Set to 0 to disable filtering on this parameter.)
output_arguments:
output_unaligned_reads_to_other_fastq: False
filter_alignments_to_softmasked_regions: False
bowtie_arguments:
m : 200
n : 1
l : 15
e : 100
trimmomatic_arguments:
LEADING: "28"
SLIDINGWINDOW: "4:20"
MINLEN: "16"
argument_order: ['LEADING', 'SLIDINGWINDOW', 'MINLEN']
low_complexity_filter_arguments:
max_low_complexity_fraction : 0.50