1、配置环境
/mnt/d/zhongyunan/genomedata/wastewater
conda create -n meta_env
conda activate meta_env
conda install -c bioconda fastp cd-hit bwa samtools megahit -y
fastp --version
megahit --version
bwa
2、下载人类基因组GRCh38，建立索引
mkdir -p ref
cd ref
wget ftp://ftp.ensembl.org/pub/release-110/fasta/homo_sapiens/dna/Homo_sapiens.GRCh38.dna.primary_assembly.fa.gz
gunzip Homo_sapiens.GRCh38.dna.primary_assembly.fa.gz
mv Homo_sapiens.GRCh38.dna.primary_assembly.fa GRCh38.fa
HUMAN_REF="/mnt/d/zhongyunan/genomedata/wastewater/ref/GRCh38.fa"
bwa index $HUMAN_REF
3、下载数据
略
查看文件数量
ls -1 | wc -l
查看文件夹大小
du -sh raw_data
