Call2plink sub-workflow

modules/local/call2plink/bedify.nf

@@ -0,0 +1,16 @@
+process bedify {
+
+   input:
+     set file(ped), file(map)
+   output:
+     file("${base}.bed"), emit: bed_ch
+     file ("${base}.bim"), emit: bim_ch
+     file("${base}.fam"), emit: fam_ch
+   script:
+      def base = ped.baseName
+      """
+      echo $base
+
+      plink --file $base --no-fid --make-bed --out $base
+      """
+  }

modules/local/call2plink/combineBed.nf

@@ -0,0 +1,19 @@
+// Combine  all plinks and possibly remove errprs
+ process combineBed {
+   input:
+     file(bed) 
+     file(bim) 
+     file(fam)
+   output:
+     tuple file("raw.bed"), file("raw.fam"), file("raw.log"), emit: plink_src
+     file("rawraw.bim"), emit: fill_in_bim_ch
+   script:
+    """
+    ls *.bed | sort > beds
+    ls *.bim | sort >  bims
+    ls *.fam | sort >  fams
+    paste beds bims fams >  mergelist
+    plink --merge-list mergelist  --make-bed --out raw
+    mv raw.bim rawraw.bim
+    """
+ } 

modules/local/call2plink/fillinbam.nf

@@ -0,0 +1,10 @@
+process fillInBim {  //  Deals with monomorphic or non-called SNPs
+    input:
+     file(inbim)
+     file(strand)
+     file(manifest)
+    output:
+     file("raw.bim"), emit: filled_bim_ch
+    script:
+       "fill_in_bim.py ${params.output_align} $strand $manifest $inbim  raw.bim"
+  }

modules/local/call2plink/illumina2lgen.nf

@@ -0,0 +1,15 @@
+process illumina2lgen {
+    input:
+    set file(report), file(array)
+
+    output:
+    tuple file("${output}.ped"), file("${output}.map"), emit: ped_ch
+    script:
+        def samplesize = params.samplesize
+	def idpat      = params.idpat
+        def output = report.baseName
+        """
+        hostname
+        topbottom.py $array $report $samplesize '$idpat' $output
+        """
+ }