-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathRun_FastQC.sh
34 lines (26 loc) · 1.04 KB
/
Run_FastQC.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
# References:
# https://www.bioinformatics.babraham.ac.uk/projects/fastqc/INSTALL.txt
# https://manpages.ubuntu.com/manpages/trusty/man1/fastqc.1.html
# Locate to the desired directory
cd scRNA/data
# View the first 6 entries
zcat HumanBrain_S1_L001_R1_001.fastq.gz | head -6
# View the number of lines in the FASTQ file
zcat HumanBrain_S1_L001_R1_001.fastq.gz | wc -l
# Before trimming: perform quality control analysis on particular or all FASTQ files stored in compressed (gzipped) format
zcat HumanBrain_S1_L001_R1_001.fastq.gz | fastqc stdin --outdir=../report/
zcat *fastq.gz | fastqc stdin --outdir=../report/
#!/bin/bash
cd ~/single_cell_data/raw_fastq
for filename in *.fastq.gz; do
echo "Processing $filename"
zcat $filename | fastqc stdin --outdir=./report/
done
# After trimming: do the FastQC again
cat *trimmed.fastq | fastqc stdin --outdir=../report/
#!/bin/bash
cd ~/single_cell_data/trim_fastq
for filename in *trimmed.fastq; do
echo "Processing trimmed $filename"
cat $filename | fastqc stdin --outdir=./report/
done