Skip to content

Commit

Permalink
Apply suggestions from code review
Browse files Browse the repository at this point in the history
Co-authored-by: Kirill Tsukanov <tskir@users.noreply.github.com>
  • Loading branch information
tcezard and tskir authored Jun 6, 2020
1 parent 8bfb9e1 commit 3a081cb
Show file tree
Hide file tree
Showing 2 changed files with 7 additions and 3 deletions.
6 changes: 5 additions & 1 deletion README_nextflow.md
Original file line number Diff line number Diff line change
Expand Up @@ -24,5 +24,9 @@ tests/test_pipeline.sh

Or with real data:
```bash
nextflow run remapping_pipeline.nf --oldgenome <genome.fa> --newgenome <new_genome.fa> --vcffile <source.vcf> --outfile <remap.vcf>
nextflow run remapping_pipeline.nf \
--oldgenome <genome.fa> \
--newgenome <new_genome.fa> \
--vcffile <source.vcf> \
--outfile <remap.vcf>
```
4 changes: 2 additions & 2 deletions remapping_pipeline.nf
Original file line number Diff line number Diff line change
Expand Up @@ -12,7 +12,7 @@ def helpMessage() {
--outfile final vcf containing variants mapped to new genome [required]
Parameters:
--flankingseq The length of the flanking sequences that generates reads (default 50)
--scorecutoff Percentage of the flanking sequences that should be used as the alignment Score cut-off threshold (default 0.6)
--scorecutoff Percentage of the flanking sequences that should be used as the alignment score cut-off threshold (default 0.6)
--diffcutoff Percentage of the flanking sequences that should be used as the AS-XS difference cut-off threshold (default 0.04)
"""
}
Expand Down Expand Up @@ -53,7 +53,7 @@ outfile_dir = file(params.outfile).getParent()
/*
* Index the new reference genome using bowtie_build
*/
if ( !file("${newgenome_dir}/${newgenome_basename}.1.bt2").exists()){
if (!file("${newgenome_dir}/${newgenome_basename}.1.bt2").exists()){

process bowtieGenomeIndex {
// Memory required is 10 times the size of the fasta in Bytes or at least 1GB
Expand Down

0 comments on commit 3a081cb

Please sign in to comment.