鸡（Gallus_gallus）转录组实战

1 软件安装

https://www.jianshu.com/p/eb89ab4af035
linux平台下需要安装的软件：fastqc，fastp，hisat2，samtools，htseq

2下载基因组序列和基因组注释文件

wget -c http://ftp.ensembl.org/pub/release-104/fasta/gallus_gallus/dna/Gallus_gallus.GRCg6a.dna_rm.toplevel.fa.gz
wget -c http://ftp.ensembl.org/pub/release-104/gtf/gallus_gallus/Gallus_gallus.GRCg6a.104.gtf.gz

解压Gallus_gallus.GRCg6a.dna_rm.toplevel.fa.gz

3构建索引文件

hisat2-build -p 2 Gallus_gallus.GRCg6a.dna_rm.toplevel.fa Gallus_gallus

4过滤raw reads

mkdir -p fastp
ls *1.fastq.gz|while read id;
do
fastp -5 20 -i ${id%_*}_1.fastq.gz -I ${id%_*}_2.fastq.gz \
-o ${id%_*}_1.clean.fq.gz -O ${id%_*}_2.clean.fq.gz \
-j ./fastp/${id%_*}.json -h ./fastp/${id%_*}.html;
done

5比对

ls *1.clean.fq.gz|while read id;
do
hisat2 -t -p 2 -x /media/lzx/0000678400004823/Indexs/Hisat2/Chicken/Gallus_gallus \
-1 $id -2 ${id%_*}_2.clean.fq.gz \
2>${id%%_*}.hisat2.log \
|samtools sort -@ 2 -o ${id%_*}_ht2p.bam
done