cd /work/home/acwnw4bl7y/HT22/HT22
ls *lite.1|while read id;
do
mv $id ${id%%.*}.sra;
done
cd /work/home/acwnw4bl7y/HT22/HT22
export PATH=/work/home/acwnw4bl7y/sratoolkit.3.1.1-centos_linux64/bin:$PATH
ls *.sra|while read id;
do
fastq-dump --gzip --split-3 $id;
done
cd /work/home/acwnw4bl7y/HT22/HT22
source ~/miniconda3/bin/activate
conda activate pyht
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
cd /work/home/acwnw4bl7y/HT22/HT22
source ~/miniconda3/bin/activate
conda activate pyht
ls *1.clean.fq.gz|while read id;
do
hisat2 -t -p 16 -x /work/home/acwnw4bl7y/GRCm39/GRCm39 \
-1 $id -2 ${id%_*}_2.clean.fq.gz \
2>${id%%_*}.hisat2.log \
|samtools sort -@ 16 -o ${id%_*}_ht2p.bam
done
cd /work/home/acwnw4bl7y/HT22/HT22
source ~/miniconda3/bin/activate
conda activate pyht
mkdir -p htseq
ls *.bam |while read id;
do
htseq-count -f bam -s no -i gene_id $id \
/work/home/acwnw4bl7y/Mus_musculus.GRCm39.106.gtf \
1>./htseq/${id%%_*}.txt 2>./htseq/${id%%_*}.HTseq.log;
done
最后編輯于 ：2024.10.18 17:03:32
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
