python qsub是啥意思_qsub|pasta|

cd /xxx/genome_stat/Annotation

ln -s /xxx/02.annotation/gff_v2/*.homolog.v2.gff /xxx/genome_stat/Annotation

ls *.gff | while read l;do echo "perl /xxx/02.annotation/build_pipeline/bin/stat_gff.pl $l" > $l.sh;done

ls *.sh | while read l;do qsub -cwd -l vf=1G,p=1 -P P18Z10200N0100 -q st.q $l;done

ls *.sh.o* | while read l;do cat $l;done > o.all

cd /xxx/genome_stat/Assembly/fa

ln -s /xxx/01.assembly/1.genomes/*.genomic.fa .

ls *.fa|while read l;do echo "sh /xxx/pipeline/Assemble/Evaluation/N50/n50.sh $l .fa";done > $l.sh

ls *.sh |while read l;do qsub -cwd -l vf=1G,p=1 -P P18Z10200N0100 -q st.q *.sh

ls *.stat | while read l; do cat $l ; done > stat.all

paste stat.all /xxx/genome_stat/Annotation/o.all > output

awk '{print $10,$1,$2,$3,$4,$5,$6,$7,$7-$8,$9,$11,$12,$13,$14,$15,$16,$17,$18}' output > excel

mv excel /xxx/genome_stat

cd /xxx/genome_stat

cat excel.txt | head -n 1 > id

cat id excel |awk '{gsub(/.homolog.v2.gff/,""); print $0}' > excel.out

perl check.pl excel.txt excel.out >err