Mercurial > repos > portiahollyoak > temp
comparison scripts/cmd.total.sh @ 0:28d1a6f8143f draft
planemo upload for repository https://github.com/portiahollyoak/Tools commit 132bb96bba8e7aed66a102ed93b7744f36d10d37-dirty
| author | portiahollyoak |
|---|---|
| date | Mon, 25 Apr 2016 13:08:56 -0400 |
| parents | |
| children |
comparison
equal
deleted
inserted
replaced
| -1:000000000000 | 0:28d1a6f8143f |
|---|---|
| 1 rm tmp | |
| 2 #for i in flamBGFM flamKGFM flamEmbryo flamHets flamTranshets | |
| 3 #for i in harwich wXh24 wXh14 wXh21 whXw14 whXw21 | |
| 4 #for i in armiTranshetsOvary armiTranshetsSoma armiHetsOvary armiHetsSoma rhinoTranshetsOvary rhinoTranshetsSoma rhinoHetsOvary rhinoHetsSoma qinTranshetsOvary qinHetsOvary w1118Ovary w1118Soma orerOvary orerSoma orerEmbryo | |
| 5 #for i in w1118Ovary w1118Soma qinTranshetsOvary qinHetsOvary qintestTranshetsOvary qintestHetsOvary | |
| 6 for i in harwich.ovary harwichG20.ovary W1.ovary W1G20.ovary wXh1g.ovary whXh3g.ovary whXh5g.ovary whXh7g.ovary whXw3g.ovary whXw5g.ovary whXw7g.ovary | |
| 7 #for i in W1.ovary wXh1g.ovary whXh3g.ovary whXh5g.ovary whXw3g.ovary whXw5g.ovary armiHets.ovary armiHets.carcass armiTranshets.ovary armiTranshets.carcass flamBGFM.ovary flam.embryo flamHets.ovary flamKGFM.ovary flamTranshets.ovary harwich.ovary introgression2.ovary introgression2X3.ovary introgression3.ovary orer.embryo orer.ovary orer.carcass qinDf.ovary qinHets.ovary qinTMB.ovary qinTranshets.ovary rhinoHets.ovary rhinoHets.carcass rhinoTranshets.ovary rhinoTranshets.carcass w1.ovary w1.carcass whXw14d.ovary whXw21d.ovary wXh14d.ovary wXh21d.ovary wXh2_4d.ovary | |
| 8 | |
| 9 do | |
| 10 | |
| 11 awk -F "\t" -v sample=$i '{OFS="\t"; print $1,$2,$3,sample,$5,$6}' $i.downsample.bam.unpair.uniq.transposons.filtered.bed >> tmp | |
| 12 | |
| 13 ## Filter BS A{36} | |
| 14 grep FBgn0000224_BS tmp | egrep "\+51|\-51" > tmp.BS | |
| 15 | |
| 16 ## Merge Stalker | |
| 17 ediff tmp diff tmp.BS > tmp2 | |
| 18 | |
| 19 done | |
| 20 | |
| 21 perl /home/wangj2/jpp_findTransposonJumping/mergeTagsWithoutGap.pl tmp2 > dysgenic.uniq.transposons.filtered.woGap.bed | |
| 22 perl /home/wangj2/jpp_findTransposonJumping/mergeTagsWithGap.pl dysgenic.uniq.transposons.filtered.woGap.bed 500 > dysgenic.uniq.transposons.filtered.wGap.bed | |
| 23 | |
| 24 rm tmp2 tmp.BS tmp | |
| 25 | |
| 26 perl get_class.pl dysgenic.uniq.transposons.filtered.wGap.bed > dysgenic.uniq.transposons.filtered.wGap.class.bed |