-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathcellranger.sh
45 lines (36 loc) · 1.63 KB
/
cellranger.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
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
######################## In vitro model ########################
### cellranger command
counts_dir="/dataOS/rcalandrelli/MARGI/RNAseq_02152019/counts/"
cd $counts_dir
# R1 and R2 fastq files from each sample to be in separate folders under "fastq_dir"
fastq_dir="/dataOS/rcalandrelli/MARGI/RNAseq_02152019/fastq"
for i in $(find $fastq_dir -maxdepth 1 -mindepth 1 -type d -printf '%f\n'); do
cellranger count --id="count_""$i" \
--transcriptome=/dataOS/rcalandrelli/Software/refdata-cellranger-GRCh38-3.0.0 \
--fastqs=/dataOS/rcalandrelli/MARGI/RNAseq_02152019/fastq/"$i" \
--sample="$i" \
--localcores=24
done
### Unzip matrices files
for i in $(find $fastq_dir -maxdepth 1 -mindepth 1 -type d -printf '%f\n'); do
cd $counts_dir$i"/outs/filtered_feature_bc_matrix/"
gunzip -k *.gz
done
######################## In vivo model ########################
### cellranger command
counts_dir="/mnt/extraids/OceanStor-SysCmn-2/rcalandrelli/MARGI/RNAseq_human_vascular/counts/"
cd $counts_dir
# R1 and R2 fastq files from each sample to be in separate folders under "fastq_dir"
fastq_dir="/mnt/extraids/OceanStor-SysCmn-2/rcalandrelli/MARGI/RNAseq_human_vascular/fastq"
for i in $(find $fastq_dir -maxdepth 1 -mindepth 1 -type d -printf '%f\n'); do
cellranger count --id="count_""$i" \
--transcriptome=/dataOS/rcalandrelli/Software/refdata-cellranger-GRCh38-3.0.0 \
--fastqs=/dataOS/rcalandrelli/MARGI/RNAseq_human_vascular/fastq/"$i" \
--sample="$i" \
--localcores=24
done
### Unzip matrices files
for i in $(find $fastq_dir -maxdepth 1 -mindepth 1 -type d -printf '%f\n'); do
cd $counts_dir$i"/outs/filtered_feature_bc_matrix/"
gunzip -k *.gz
done