Skip to content
Open
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
1 change: 0 additions & 1 deletion input/sampledata.txt
Original file line number Diff line number Diff line change
@@ -1,5 +1,4 @@
SampleName SampleID AssemblyGroup FastqFile
HN008_subsampled HN008 samplegroup1 HN008/HN008_subsampled
HN009_subsampled HN009 samplegroup2 HN009_subsampled
HN016_subsampled HN016 samplegroup1 HN016_subsampled
HN036_subsampled HN036 samplegroup1 HN036_subsampled
Expand Down
4 changes: 2 additions & 2 deletions rules/02-trim.smk
Original file line number Diff line number Diff line change
Expand Up @@ -71,8 +71,8 @@ rule trimmomatic:
"02-trim",\
"{sample}_2.unpaired.fastq.gz")
log:
err = os.path.join(OUTPUTDIR, "logs",\
"01-setup", "02-trim",\
err = os.path.join(OUTPUTDIR, "logs",\
"01-setup", "02-trim",\
"PE_{sample}_err.log"),
out = os.path.join(OUTPUTDIR, "logs",\
"01-setup", "02-trim",\
Expand Down
4 changes: 2 additions & 2 deletions rules/11-SWAM.smk
Original file line number Diff line number Diff line change
Expand Up @@ -41,9 +41,9 @@ rule merge_all:
'''
echo "" > {output}
for assemblyfile in {params.assemblyfiles}; do
currvar=$(echo $assemblyfile | rev | cut -f1 -d\/ | rev)
currvar=$(echo $assemblyfile | rev | cut -f1 -d\\/ | rev)
sed 's/>.*/&_'"$currvar"'/' $assemblyfile >> {output}
done
#cat {params.assemblyfiles} > {output}
'''


2 changes: 1 addition & 1 deletion rules/12-MAD.smk
Original file line number Diff line number Diff line change
Expand Up @@ -116,7 +116,7 @@ rule replace_mad_ids:
"rename","MAD.{filter_workflow}.fasta")
run:
sequences=[curr for curr in SeqIO.parse(input.mad,"fasta")]
concordance=pd.read_csv(input.concordance,sep="\s+")
concordance=pd.read_csv(input.concordance,sep="\\s+")
concord_dict=dict(zip(concordance["MAD_id"],concordance["Concordance_Short_ID"]))
for i,curr in enumerate(sequences):
curr.id = concord_dict[curr.id]
Expand Down