#PBS -N 4_Sacc_cerevisiae_dec_cov1_ymj
#PBS -l nodes=1:ppn=48
#PBS -r y
#PBS -q master_431

#转录组不需要做这一步，直接跳到下一步

cd /home/Data_01/yizhenzhen06/public_share/Saccharomyces_cerevisiae_DNA/3_decontamination/

echo Start Time is `date`

seqkit seq -n -i ../2_assembly/spades/scaffolds.fasta > 4_Sacc_cerevisiae_dec_all.txt
awk '{if($1~"_0.")print$1}' 4_Sacc_cerevisiae_dec_all.txt > 4_Sacc_cerevisiae_dec_all_cov1.id
seqkit grep -v -f 4_Sacc_cerevisiae_dec_all_cov1.id 4_Sacc_cerevisiae_dec_all.fasta > 4_Sacc_cerevisiae_dec_all_cov1.fasta 

echo End Time is `date`