Skip to content

Commit

Permalink
Merge pull request #24 from MattBixley/patch-3
Browse files Browse the repository at this point in the history
Update 3_trimmingfiltering.md
  • Loading branch information
DininduSenanayake authored Mar 11, 2024
2 parents b1aebb2 + a9b3b62 commit 578b2f4
Showing 1 changed file with 15 additions and 21 deletions.
36 changes: 15 additions & 21 deletions docs/3_trimmingfiltering.md
Original file line number Diff line number Diff line change
Expand Up @@ -35,32 +35,30 @@ from FASTQ files.
- The following syntax will remove the adapter sequence AACCGGTT from the file SRR014335-chr1.fastq, create a new file called SRR014335-chr1_trimmed.fastq, and write a summary to the log file SRR014335-chr1.log:

```bash

$ pwd
/home/$USER/RNA_seq

$ mkdir Trimmed

$ module load cutadapt/4.1-gimkl-2022a-Python-3.10.5

$ cutadapt -q 20 -a AACCGGTT -o Trimmed/SRR014335-chr1_cutadapt.fastq Raw/SRR014335-chr1.fastq > Trimmed/SRR014335-chr1.log

pwd
```
```
/home/$USER/RNA_seq
```
```bash
mkdir Trimmed
module load cutadapt/4.1-gimkl-2022a-Python-3.10.5
cutadapt -q 20 -a AACCGGTT -o Trimmed/SRR014335-chr1_cutadapt.fastq Raw/SRR014335-chr1.fastq > Trimmed/SRR014335-chr1.log
```
!!! quote ""
* **-q** (`--quality-cutoff`) parameter can be used to trim low-quality ends from reads. If you specify a single cutoff value, the 3’ end of each read is trimmed.

```bash

$ less Trimmed/SRR014335-chr1.log

less Trimmed/SRR014335-chr1.log
```

Now we should trim all samples.

```bash
$ cd Raw

$ ls
cd Raw
ls
```
```
SRR014335-chr1.fastq SRR014336-chr1.fastq SRR014337-chr1.fastq SRR014339-chr1.fastq SRR014340-chr1.fastq SRR014341-chr1.fastq
```
```bash
Expand All @@ -79,16 +77,12 @@ done

```bash
cd ../MultiQC
```
```bash
cp ../Trimmed/*log .
```
```bash
multiqc .
```
![image](./Images/MQC2.png)

- - -


<p align="center"><b><a class="btn" href="https://genomicsaotearoa.github.io/RNA-seq-workshop/" style="background: var(--bs-dark);font-weight:bold">Back to homepage</a></b></p>
<p align="center"><b><a class="btn" href="https://genomicsaotearoa.github.io/RNA-seq-workshop/" style="background: var(--bs-dark);font-weight:bold">Back to homepage</a></b></p>

0 comments on commit 578b2f4

Please sign in to comment.