From 93b01501a8b5804fbc12eef0380fc3695e5bdbfd Mon Sep 17 00:00:00 2001
From: YuuMJ <76215523+YuuMJ@users.noreply.github.com>
Date: Sun, 1 Jun 2025 16:27:48 +0800
Subject: [PATCH] update

---
 .../4_Sacc_cerevisiae_dec_cov1_ymj.pbs        | 20 +++++++++++++++++++
 ...��先检查路径和文件名是否修改正确！！！解决不了再来问我 |  0
 2 files changed, 20 insertions(+)
 create mode 100644 3_decontamination/4_Sacc_cerevisiae_dec_cov1_ymj.pbs
 create mode 100644 注意！！！如果有报错，请先检查路径和文件名是否修改正确！！！解决不了再来问我

diff --git a/3_decontamination/4_Sacc_cerevisiae_dec_cov1_ymj.pbs b/3_decontamination/4_Sacc_cerevisiae_dec_cov1_ymj.pbs
new file mode 100644
index 0000000..54b3423
--- /dev/null
+++ b/3_decontamination/4_Sacc_cerevisiae_dec_cov1_ymj.pbs
@@ -0,0 +1,20 @@
+#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`
+
+blastn -db /home/Data_01/share_databases/human_genome_new/GCF_000001405.39_GRCh38.p13_rna_db \
+-query 2_Sacc_cerevisiae_dec_bac_hum.fasta \
+-out 3_Sacc_cerevisiae_dec_bac_hum_geno_rna.tab \
+-outfmt "6 qseqid qlen sseqid slen pident qcovs qstart qend sstart send mismatch gapopen evalue bitscore" \
+-evalue 1e-5 -max_target_seqs 1 -num_threads 48
+
+sed -i '1i\tquery ID\tquery sequence length\tsubject ID\tsubject qequence length\tpercent identities\tcoverage\tqstart\tqend\tsstart\tsend\tmismatch\tgapopen\tevalue\tbitscores ' 3_Sacc_cerevisiae_dec_bac_hum_geno_rna.tab
+
+echo End Time is `date`
diff --git a/注意！！！如果有报错，请先检查路径和文件名是否修改正确！！！解决不了再来问我 b/注意！！！如果有报错，请先检查路径和文件名是否修改正确！！！解决不了再来问我
new file mode 100644
index 0000000..e69de29