#========== snap predictions/ training per species ======== set dp=dvir set species=virilis set dpid=${dp}_caf051209 set scd=$sc/${dp}3 set scaffoldgff=$scd/$dp-scaffolds.gff set chrfasta=$scd/$dpid.fa.gz (source snapmake.script) >& log.snapmake.$dpid & #---- set dp=dgri set species=grimshawi set dpid=${dp}_caf051209 set scd=$sc/${dp}3 set scaffoldgff=$scd/$dp-scaffolds.gff set chrfasta=$scd/$dpid.fa.gz (source snapmake.script) >& log.snapmake.$dpid & #---- set dp=dmoj set species=mojavensis set dpid=${dp}_caf051209 set scd=$sc/${dp}3 set scaffoldgff=$scd/$dp-scaffolds.gff set chrfasta=$scd/$dpid.fa.gz (source snapmake.script) >& log.snapmake.$dpid & #---- set dp=dana set species=ananassae set dpid=${dp}_caf051209 set scd=$sc/${dp}3 set scaffoldgff=$scd/$dp-scaffolds.gff set chrfasta=$scd/$dpid.fa.gz (source snapmake.script) >& log.snapmake.$dpid & #---- set dp=dere set species=erecta set dpid=${dp}_caf051209 set scd=$sc/${dp}3 set scaffoldgff=$scd/$dp-scaffolds.gff set chrfasta=$scd/$dpid.fa.gz (source snapmake.script) >& log.snapmake.$dpid & #---- set dp=dsec set species=sechellia set dpid=${dp}_br051028 set scd=$sc/${dp}1 set scaffoldgff=$scd/$dp-scaffolds.gff set chrfasta=$scd/$dpid.fa.gz (source snapmake.script) >& log.snapmake.$dpid & #---- set dp=dper set species=persimilis set dpid=${dp}_br051028 set scd=$sc/${dp}1 set scaffoldgff=$scd/$dp-scaffolds.gff set chrfasta=$scd/$dpid.fa.gz (source snapmake.script) >& log.snapmake.$dpid & #---- set dp=dyak set species=yakuba set dpid=${dp}_caf051213 set scd=$sc/${dp}3 set scaffoldgff=$scd/$dp-chromosomes.gff set chrfasta=$scd/$dpid.fa.gz (source snapmake.script) >& log.snapmake.$dpid & #---- set dp=dpse set species=pseudoobscura set dpid=${dp}_r201 set scd=$sc/${dp}2 set scaffoldgff=$scd/$dp-chromosomes.gff set chrfasta=$scd/$dp-chromosomes-r2.0.fa.gz set trainfasta=$scd/dpse_r201shred5.fa (source snapmake.script) >& log.snapmake.$dpid & #---- set dp=dsim set species=simulans set dpid=${dp}_wu050602 set scd=$sc/${dp}2 set scaffoldgff=$scd/$dp-chromosomes.gff set chrfasta=$scd/$dpid.fa.gz set trainfasta=$scd/xxx (source snapmake.script) >& log.snapmake.$dpid & #----