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
