forked from omics-school/analyse-rna-seq
-
Notifications
You must be signed in to change notification settings - Fork 1
/
script2.sh
43 lines (36 loc) · 2.06 KB
/
script2.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
# Référence de l'échantillon à analyser.
sample="SRR2960338"
# Chemin relatif et nom du fichier contenant le génome de référence.
genome="genome/GCF_000214015.3_version_140606.fna"
# Chemin relatif et nom du fichier contenant les annotations.
annotations="genome/GCF_000214015.3_version_140606.gff"
echo "=============================================================="
echo "Contrôle qualité"
echo "=============================================================="
mkdir -p reads_qc
fastqc "reads/${sample}.fastq.gz" --outdir reads_qc
echo "=============================================================="
echo "Indexation du génome de référence"
echo "=============================================================="
mkdir -p index
bowtie2-build "${genome}" "index/O_tauri"
echo "=============================================================="
echo "Alignement des reads sur le génome de référence"
echo "=============================================================="
mkdir -p map
bowtie2 -p 2 -x "index/O_tauri" -U "reads/${sample}.fastq.gz" -S "map/bowtie-${sample}.sam" 2> "map/bowtie-${sample}.out"
echo "=============================================================="
echo "Conversion en binaire, tri et indexation des reads alignés"
echo "=============================================================="
samtools view -@ 2 -b "map/bowtie-${sample}.sam" -o "map/bowtie-${sample}.bam"
samtools sort -@ 2 "map/bowtie-${sample}.bam" -o "map/bowtie-${sample}.sorted.bam"
samtools index "map/bowtie-${sample}.sorted.bam"
echo "=============================================================="
echo "Comptage"
echo "=============================================================="
mkdir -p count
htseq-count --stranded=no --type="gene" --idattr="ID" --order=name --format=bam "map/bowtie-${sample}.sorted.bam" "${annotations}" > "count/count-${sample}.txt"
echo "=============================================================="
echo "Nettoyage des fichiers inutiles"
echo "=============================================================="
rm -f "map/bowtie-${sample}.sam" "map/bowtie-${sample}.bam"