File tree Expand file tree Collapse file tree
Expand file tree Collapse file tree Original file line number Diff line number Diff line change @@ -163,14 +163,14 @@ LRSCs=$(readlink -f "$0")
163163LRSCf=$( dirname $LRSCs )
164164
165165echo " [" $( date) " ] Self-aligning the long reads (minimap2)"
166- # if [[ $minimapOptions == "PB" ]] ; then
167- # $LRSCf/minimap2/minimap2 --dual=yes -PD --no-long-join -w5 -g1000 -m30 -n1 -I100G -t"$nproc" "$reads" "$reads" > $tmpdir/"$alignments"
168- # else
169- # $LRSCf/minimap2/minimap2 -k15 -w5 -m100 -g10000 -r2000 --max-chain-skip 25 --dual=yes -PD --no-long-join -I100G -t"$nproc" "$reads" "$reads" > $tmpdir/"$alignments"
170- # fi
166+ if [[ $minimapOptions == " PB" ]] ; then
167+ $LRSCf /minimap2/minimap2 --dual=yes -PD --no-long-join -w5 -g1000 -m30 -n1 -I100G -t" $nproc " " $reads " " $reads " > $tmpdir /" $alignments "
168+ else
169+ $LRSCf /minimap2/minimap2 -k15 -w5 -m100 -g10000 -r2000 --max-chain-skip 25 --dual=yes -PD --no-long-join -I100G -t" $nproc " " $reads " " $reads " > $tmpdir /" $alignments "
170+ fi
171171
172172echo " [" $( date) " ] Correcting the long reads"
173- $LRSCf /bin/CONSENT -a Alignments_16868.paf -s " $minSupport " -S " $maxSupport " -l " $windowSize " -k " $merSize " -c " $commonKMers " -A " $minAnchors " -f " $solid " -m " $windowOverlap " -j " $nproc " -r " $reads " -M " $maxMSA " -p " $LRSCf " >> " $out "
173+ $LRSCf /bin/CONSENT -a $tmpdir / " $alignments " -s " $minSupport " -S " $maxSupport " -l " $windowSize " -k " $merSize " -c " $commonKMers " -A " $minAnchors " -f " $solid " -m " $windowOverlap " -j " $nproc " -r " $reads " -M " $maxMSA " -p " $LRSCf " >> " $out "
174174echo " [" $( date) " ] Removing the temporary files"
175175rm $tmpdir /" $alignments "
176176echo " [" $( date) " ] Exiting"
You can’t perform that action at this time.
0 commit comments