You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
getAssemblyInfo(){
#Have user estimate/enter assembly parameters if unenteredif [ -z"$CUTOFF" ];thenforiin {2..20};doecho$i>> pfile
done
cat pfile | parallel -j $NUMProc --no-notice "echo -n {}xxx && mawk -v x={} '\$1 >= x' uniq.seqs | wc -l"| mawk '{gsub("xxx","\t",$0); print;}'|$sort -g > uniqseq.data
rm pfile
#Plot graph of above data
gnuplot << \EOF set terminal dumb size 120, 30 set autoscale set xrange [2:20] unset label set title "Number of Unique Sequences with More than X Coverage (Counted within individuals)" set xlabel "Coverage" set ylabel "Number of Unique Sequences" set lmargin 10 plot 'uniqseq.data' with lines notitle pause -1EOFecho -en "\007"echo -en "\007"echo -en "\007"echo -e "Please choose data cutoff. In essence, you are picking a minimum (within individual) coverage level for a read (allele) to be used in the reference assembly"read CUTOFF
fi
The text was updated successfully, but these errors were encountered:
cbird808
changed the title
getAssemblyInfo is missing a }getAssemblyInfo is missing a closing }Sep 22, 2019
The text was updated successfully, but these errors were encountered: