🧬 Advanced linux tutorial

The dataset for this exercise is available at

/data/PwaniGenomics/fastq_data

#create symbolic link to the data in your user directory

ln -s /data/PwaniGenomics/fastq_data

Note: FASTQ reads are in .gz format.

Genome File Descriptions (.fna Indexes)

• fna.amb – Contains info on ambiguous bases (e.g., N).
• fna.ann – Contains metadata such as sequence names and lengths.
• fna.bwt – Burrows-Wheeler transformed sequence (compressed format).
• fna.pac – Packed sequence data.
• fna.sa – Suffix array used to locate sequences.

Structure of a FASTQ File

FASTQ files (typically .fastq.gz) contain sequencing reads. Each read consists of four lines:

@SRR15369215.126490887       # Sequence identifier
GGACCTTCTGTCATTTCACT...      # Nucleotide sequence
+                             # Separator
AAFFFJJJJJJJJJJJJJJJJ...      # Base quality scores
  

Read Length Analysis

Count Total Reads for SRR836370 and SRR836361:

Reads Shorter than 100bp for SRR836370 and SRR836361:

Reads Longer than 100bp for SRR836370 and SRR836361:

Reads Not Equal to 100bp for SRR836370 and SRR836361:

Reads Contaminated with Illumina Adapter for SRR836370 and SRR836361:

Adapter sequence: CTGTCTCTTATACACATCT

Extra note

The Datasets can be accessed later too

Datasets are available at: https://zenodo.org/records/21871083