perl/etc/justonce/tmp/exSubfaBlast.sh

   1 #!/bin/bash
   2
   3 samtools faidx Obimaculoides_280.fa $1 >$1.fa
   4 blastn -task blastn -query ir813.fa -subject $1.fa -word_size 4 -outfmt 3 >$1.b3.txt
   5 zcat annotation/Obimaculoides_280.gene.gff3.gz|egrep "^$1\s" >$1.gff
   6
   7 # cat matches.lst|while read a;do ./t.sh $a;done
   8 # rar a -m5 -rr1% -mct -mcc -mt1 res.rar ir813.fa ir813.e10w4.txt matches.lst Scaffold* ir813.e1000w4.txt
   9 # $ cat matches.lst
  10 # Scaffold17571
  11 # Scaffold73065
  12 # Scaffold54938
  13