# HG changeset patch # User mingchen0919 # Date 1523298042 14400 # Node ID 38435c74e795c003d6d9eeeb1228bc8279d6df10 # Parent ecf34bd5298790e28f0972c8f9489af621353d75 working version diff -r ecf34bd52987 -r 38435c74e795 split.pl --- a/split.pl Mon Apr 09 13:32:55 2018 -0400 +++ b/split.pl Mon Apr 09 14:20:42 2018 -0400 @@ -24,7 +24,7 @@ for($i = 1; $i <= $numOfFiles; $i++){ print "$i\n"; - open FILE, ">".$homfile.".".$i || die "Can't open file"; + open FILE, ">".$homfile.".".$i.".fasta" || die "Can't open file"; print FILE $lineIn; $seqs = 1; $lineIn = ; @@ -40,7 +40,7 @@ close FILE; } $i = $i -1; - open FILE, ">>".$homfile.".".$i; + open FILE, ">>".$homfile.".".$i.".fasta"; while ($lineIn = ){ print FILE $lineIn; } diff -r ecf34bd52987 -r 38435c74e795 split_fasta.Rmd --- a/split_fasta.Rmd Mon Apr 09 13:32:55 2018 -0400 +++ b/split_fasta.Rmd Mon Apr 09 14:20:42 2018 -0400 @@ -16,9 +16,9 @@ # single-end.sh cat <${X_d}/job-script.sh -cd ${WORKING_DIR}/fasta_files && ${X_t}/split.pl \\ - --input_file=${X_A} \\ - --seqs_per_file=${X_B} \\ > ${X_d}/fasta_splitter-log.txt 2>&1 +cd ${WORKING_DIR}/fasta_files && \\ + cp ${X_A} input.fasta && \\ + ${X_t}/split.pl input.fasta ${X_B} > ${X_d}/fasta_splitter-log.txt 2>&1 EOF ``` @@ -28,6 +28,12 @@ sh ${X_d}/job-script.sh ``` +```{bash} +# remove original input fasta +rm ${WORKING_DIR}/fasta_files/input.fasta +mv ${WORKING_DIR}/fasta_files/out ${X_d}/split-out.txt +``` + ```{r, 'display output directory contents', results='asis', echo=FALSE} ## after the job is done, we list all files from the output directory. ## full relative path to the output directory needs to be displayed. diff -r ecf34bd52987 -r 38435c74e795 split_fasta.xml --- a/split_fasta.xml Mon Apr 09 13:32:55 2018 -0400 +++ b/split_fasta.xml Mon Apr 09 14:20:42 2018 -0400 @@ -16,9 +16,9 @@ -B '$number' ]]> - + - +