Skip to content

Commit

Permalink
Fixes to GMAP parameters
Browse files Browse the repository at this point in the history
- set --max-intronlength-ends to 500kb, this is due to cases
  where GMAP soft-clips the terminal exon(s) if they are over
  10kb (default value)
- set chimeric alignment option to match min clip length
  • Loading branch information
mcmero committed Jun 24, 2020
1 parent bbde885 commit 9604ffb
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion MINTIE.groovy
Original file line number Diff line number Diff line change
Expand Up @@ -197,7 +197,7 @@ align_contigs_against_genome = {
output.dir = sample_name
produce('aligned_contigs_against_genome.sam'){
exec """
$gmap -D $gmap_refdir -d $gmap_genome -f samse -t $threads -n 0 $input.fasta > $output
$gmap -D $gmap_refdir -d $gmap_genome -f samse -t $threads -x $min_gap --max-intronlength-ends=500000 -n 0 $input.fasta > $output
""", "align_contigs_against_genome"
}
}
Expand Down

0 comments on commit 9604ffb

Please sign in to comment.