Commit 78d1317c authored by Laura Denies's avatar Laura Denies
Browse files

Update ownHMM_library.R

parent 0028bbdd
......@@ -7,11 +7,10 @@
### Input of the different files ###
####################################
gene_table_HMM<-read.delim(file = snakemake@input[[1]], header=FALSE, comment.char="#")
gene_table_HMM<-read.delim(file = snakemake@input[[1]], comment.char="#")
translation<-read.delim(file = snakemake@input[[2]], header=FALSE)
signalP<-read.table(file = snakemake@input[[3]], sep= "\t")
colnames(gene_table_HMM)<-c("Query_sequence","HMM_Name","Significance_Evalue","Score")
gene_table_HMM<-gene_table_HMM[,c(1,2,4,3)]
colnames(signalP)<-c("Query_sequence","SignalP")
......@@ -55,4 +54,3 @@ gene_table_Toxic<-merge.data.frame(gene_table_Toxic,signalP,by="Query_sequence")
gene_table_Toxic$Query_sequence<-sprintf("%010d",as.numeric(gene_table_Toxic$Query_sequence))
write.csv(gene_table_Toxic,file = snakemake@output[[2]])
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment