add: spades de novo

RGBEPP.d

@@ -16,6 +16,7 @@ void show_help(string pkgver) {
 	    Version: ", pkgver, "
 	    License: GPL-2.0-only
 	    Author: Guoyi Zhang
+	    -c\t--config\tconfig file for software path (optional)
 	    -g\t--genes\t\tgene file path (optional, if -r is specified)
 	    -f\t--functions\tfunctions type (optional): all clean map 
 	      \t           \tpostmap varcall consen align
@@ -24,7 +25,13 @@ void show_help(string pkgver) {
 	    -m\t--memory\tmemory settings (optional, default 16 GB)
 	    -r\t--reference\treference genome path
 	    -t\t--threads\tthreads setting (optional, default 8 threads)
-	    --macse\t\tMacse jarfile path
+	    --fastp\t\tFastp path (optional)
+	    --bowtie2\t\tBowtie2 path (optional)
+	    --samtools\t\tSamtools path (optional)
+	    --bcftools\t\tBcftools path (optional)
+	    --macse\t\tMacse jarfile path (optional)
+	    --trimal\t\tTrimal path (optional)
+	    --spades\t\tSpades python path (optional)
 	    for example: ./RGBEPP -f all -l list -t 8 -r reference.fasta \n");
 }
 
@@ -108,7 +115,22 @@ string[] getARG_G(string ARG_R){
     return ARG_G;
 }
 
-void processQualityControl(string[] ARG_L, int ARG_T, string DirRaw, string DirQcTrim, string PathFastp) {
+string getValueFromConfig(string file, string key) {
+    string content = readText(file);
+    string value;
+    auto regex = regex(key ~ r"\s*=\s*(.+)");
+
+    foreach (line; content.splitter("\n")) {
+        if (auto match = matchFirst(line, regex)) {
+            value = match.captures[1];
+            break;
+        }
+    }
+
+    return value;
+}
+
+void processQcTrim(string[] ARG_L, int ARG_T, string DirRaw, string DirQcTrim, string PathFastp) {
     // Prepare directory
     createDir(DirQcTrim);
     writeln("QcTrimming::Start");
@@ -313,22 +335,19 @@ void processAlign(string[] ARG_G, string DirConsensus, string DirAlign, string P
 
 }
 
-string getValueFromConfig(string file, string key) {
-    string content = readText(file);
-    string value;
-    auto regex = regex(key ~ r"\s*=\s*(.+)");
-
-    foreach (line; content.splitter("\n")) {
-        if (auto match = matchFirst(line, regex)) {
-            value = match.captures[1];
-            break;
+void processAssembly(string[] ARG_L, int ARG_M, int ARG_T, string DirQcTrim, string DirAssembly, string PathSpades){
+    createDir(DirAssembly);
+    foreach (string file; ARG_L) {
+       string baseName = getBaseName(file);
+       string DirAss = DirAssembly ~ "/" ~ baseName;
+       createDir(DirAss);
+       string inputFileR1 = DirQcTrim ~ "/" ~ baseName ~ "_R1.fastq.gz";
+       string inputFileR2 = DirQcTrim ~ "/" ~ baseName ~ "_R2.fastq.gz";
+       string[] cmdAssembly = [PathSpades, "--pe1-1", inputFileR1, "--pe1-2", inputFileR2, "-t", ARG_T.to!string, "-m", ARG_M.to!string, "--careful", "--phred-offset", "33", "-o", DirAss];
+    	executeCommand(cmdAssembly);
     }
 }
 
-    return value;
-}
-
-
 void main(string[] args) {
     string pkgver = "0.0.3";
 
@@ -336,6 +355,7 @@ void main(string[] args) {
     string DirRaw = DirHome ~ "/00_raw";
     string DirQcTrim = DirHome ~ "/01_fastp";
     string DirMap = DirHome ~ "/02_bowtie2";
+    string DirAssembly = DirHome ~ "/02_spades";
     string DirBam = DirHome ~ "/03_bam";
     string DirVcf = DirHome ~ "/04_vcf";
     string DirConsensus = DirHome ~ "/05_consen";
@@ -347,8 +367,10 @@ void main(string[] args) {
     string PathBcftools = "/usr/bin/bcftools";
     string PathMacse = "/usr/share/java/macse.jar";
     string PathTrimal = "/usr/bin/trimal";
+    string PathSpades = "/usr/bin/spades.py";
 
     int ARG_T = 8;
+    int ARG_M = 16;
     string[] ARG_G;
     string[] ARG_L;
     string ARG_C;
@@ -385,10 +407,34 @@ void main(string[] args) {
 		    i++;
                     ARG_T = args[i].to!int;
                     break;
+                case "--fastp":
+		    i++;
+                    PathFastp = args[i];
+                    break;
+                case "--bowtie2":
+		    i++;
+                    PathBowtie2 = args[i];
+                    break;
+                case "--samtools":
+		    i++;
+                    PathSamtools = args[i];
+                    break;
+                case "--bcftools":
+		    i++;
+                    PathBcftools = args[i];
+                    break;
                 case "--macse":
 		    i++;
                     PathMacse = args[i];
                     break;
+                case "--trimal":
+		    i++;
+                    PathTrimal = args[i];
+                    break;
+                case "--spades":
+		    i++;
+                    PathSpades = args[i];
+                    break;
                 default:
                     break;
             }
@@ -418,7 +464,11 @@ void main(string[] args) {
     writeln("RGBEPP::Start");
     // Perform steps based on provided function argument
     if (ARG_F == "all" || ARG_F == "clean") {
-        processQualityControl(ARG_L, ARG_T, DirRaw, DirQcTrim, PathFastp);
+        processQcTrim(ARG_L, ARG_T, DirRaw, DirQcTrim, PathFastp);
+    }
+
+    if (ARG_F == "assembly") {
+	processAlign(ARG_G, DirConsensus, DirAlign, PathMacse);
     }
 
     if (ARG_F == "all" || ARG_F == "map") {
@@ -442,6 +492,8 @@ void main(string[] args) {
 	processAlign(ARG_G, DirConsensus, DirAlign, PathMacse);
     }
 
+
+
     writeln("RGBEPP::End");
 }