瀏覽代碼

added question 6

nicolas-zimmermann 5 年之前
父節點
當前提交
c1f3f77edf
共有 1 個文件被更改,包括 3 次插入3 次删除
  1. 3 3
      wgs.sh

+ 3 - 3
wgs.sh 查看文件

@@ -28,12 +28,12 @@ soft/megahit --k-list 21 --mem-flag 0 -1 EchE_R1.fastq.gz -2 EchE_R2.fastq.gz -o
28 28
 
29 29
 # Question 4
30 30
 
31
-soft/prodigal -i contigs/final.contigs.fa -d genes.fna
31
+soft/prodigal -i contigs/final.contigs.fa -d results/genes.fna
32 32
 
33 33
 # Question 5
34 34
 
35
-sed "s:>:*\n>:g" genes.fna | sed -n "/partial=00/,/*/p"|grep -v "*" > genes_full.fna
35
+sed "s:>:*\n>:g" results/genes.fna | sed -n "/partial=00/,/*/p"|grep -v "*" > results/genes_full.fna
36 36
 
37 37
 # Question 6
38 38
 
39
-
39
+blastn -query results/genes_full.fna -db databases/resfinder.fna -perc_identity 0.8 -qcov_hsp_perc 0.8 -evalue 1e-3 -outfmt 6 -out results/annoted_genes.dat