@@ -10,37 +10,37 @@ recipe:
10
10
- |
11
11
mkdir -p srnaseq
12
12
cd srnaseq
13
- wget --random-wait --retry-connrefused -nv -N - c -O hsahg19.gff3 ftp://mirbase.org/pub/mirbase/20/genomes/hsa.gff3
13
+ wget --random-wait --retry-connrefused -nv -c -O hsahg19.gff3 ftp://mirbase.org/pub/mirbase/20/genomes/hsa.gff3
14
14
awk '$3=="miRNA"' hsahg19.gff3 | sed 's/=/ /g' > srna-transcripts.gtf
15
- wget --random-wait --retry-connrefused -nv -N - c -O wgRna.txt.gz http://hgdownload.soe.ucsc.edu/goldenPath/hg19/database/wgRna.txt.gz
15
+ wget --random-wait --retry-connrefused -nv -c -O wgRna.txt.gz http://hgdownload.soe.ucsc.edu/goldenPath/hg19/database/wgRna.txt.gz
16
16
zgrep -v 'hsa-' wgRna.txt.gz | awk '{print $2"\t.\tncrna\t"$3"\t"$4"\t.\t"$7"\t.\tname "$5";"}' >> srna-transcripts.gtf
17
- wget --random-wait --retry-connrefused -nv -N - c -O tRNAs.txt.gz http://hgdownload.soe.ucsc.edu/goldenPath/hg19/database/tRNAs.txt.gz
17
+ wget --random-wait --retry-connrefused -nv -c -O tRNAs.txt.gz http://hgdownload.soe.ucsc.edu/goldenPath/hg19/database/tRNAs.txt.gz
18
18
zcat tRNAs.txt.gz | awk '{print $2"\t.\ttRNA\t"$3"\t"$4"\t.\t"$7"\t.\tname "$5";"}' >> srna-transcripts.gtf
19
- wget --random-wait --retry-connrefused -nv -N - c -O rmsk.txt.gz http://hgdownload.soe.ucsc.edu/goldenPath/hg19/database/rmsk.txt.gz
19
+ wget --random-wait --retry-connrefused -nv -c -O rmsk.txt.gz http://hgdownload.soe.ucsc.edu/goldenPath/hg19/database/rmsk.txt.gz
20
20
zcat rmsk.txt.gz | awk '{print $6"\t.\trepeat\t"$7+1"\t"$8+1"\t.\t"$10"\t.\tname "$12";"}' >> srna-transcripts.gtf
21
- wget --random-wait --retry-connrefused -nv -N - c -O refGene.txt.gz http://hgdownload.soe.ucsc.edu/goldenPath/hg19/database/refGene.txt.gz
21
+ wget --random-wait --retry-connrefused -nv -c -O refGene.txt.gz http://hgdownload.soe.ucsc.edu/goldenPath/hg19/database/refGene.txt.gz
22
22
zcat refGene.txt.gz | awk '{print $3"\t.\tgene\t"$5"\t"$6"\t.\t"$4"\t.\tname "$13";"}' >> srna-transcripts.gtf
23
- # wget --random-wait --retry-connrefused -nv -N - c -O piR_hg19_v1.0.bed.gz http://www.regulatoryrna.org/database/piRNA/download/archive/v1.0/bed/piR_hg19_v1.0.bed.gz
23
+ # wget --random-wait --retry-connrefused -nv -c -O piR_hg19_v1.0.bed.gz http://www.regulatoryrna.org/database/piRNA/download/archive/v1.0/bed/piR_hg19_v1.0.bed.gz
24
24
# zcat piR_hg19_v1.0.bed.gz | awk '{print $1"\t.\tpiRNA\t"$2"\t"$3"\t.\t"$6"\t.\tname "$4";"}' >> srna-transcripts.gtf
25
25
# mirbase
26
- wget --random-wait --retry-connrefused -nv -N - c -O hairpin.fa.gz ftp://mirbase.org/pub/mirbase/21/hairpin.fa.gz
26
+ wget --random-wait --retry-connrefused -nv -c -O hairpin.fa.gz ftp://mirbase.org/pub/mirbase/21/hairpin.fa.gz
27
27
zcat hairpin.fa.gz | awk '{if ($0~/>hsa/){name=$0; print name} else if ($0~/^>/){name=0};if (name!=0 && $0!~/^>/){print $0;}}' | sed 's/U/T/g' > hairpin.fa
28
- wget --random-wait --retry-connrefused -nv -N - c -O mature.fa.gz ftp://mirbase.org/pub/mirbase/21/mature.fa.gz
28
+ wget --random-wait --retry-connrefused -nv -c -O mature.fa.gz ftp://mirbase.org/pub/mirbase/21/mature.fa.gz
29
29
zcat mature.fa.gz | awk '{if ($0~/>hsa/){name=$0; print name} else if ($0~/^>/){name=0};if (name!=0 && $0!~/^>/){print $0;}}' | sed 's/U/T/g' > mature.fa
30
- wget --random-wait --retry-connrefused -nv -N - c -O miRNA.str.gz ftp://mirbase.org/pub/mirbase/21/miRNA.str.gz
30
+ wget --random-wait --retry-connrefused -nv -c -O miRNA.str.gz ftp://mirbase.org/pub/mirbase/21/miRNA.str.gz
31
31
zcat miRNA.str.gz | awk '{if ($0~/hsa/)print $0}' > miRNA.str
32
- wget --random-wait --retry-connrefused -nv -N - c -O mirbase.gff3 ftp://mirbase.org/pub/mirbase/21/genomes/hsa.gff3
32
+ wget --random-wait --retry-connrefused -nv -c -O mirbase.gff3 ftp://mirbase.org/pub/mirbase/21/genomes/hsa.gff3
33
33
# tdrmapper
34
- wget --random-wait --retry-connrefused --no-check-certificate -nv -N - c -O trna_mature_pre.fa https://github.com/sararselitsky/tDRmapper/raw/master/hg19_mature_and_pre.fa
34
+ wget --random-wait --retry-connrefused --no-check-certificate -nv -c -O trna_mature_pre.fa https://github.com/sararselitsky/tDRmapper/raw/master/hg19_mature_and_pre.fa
35
35
# mintmap
36
36
wget --random-wait --retry-connrefused --no-check-certificate -nv -N -c https://github.com/TJU-CMC-Org/MINTmap/raw/master/LookupTable.tRFs.MINTmap_v1.txt
37
37
wget --random-wait --retry-connrefused --no-check-certificate -nv -N -c https://github.com/TJU-CMC-Org/MINTmap/raw/master/tRNAspace.Spliced.Sequences.MINTmap_v1.fa
38
38
wget --random-wait --retry-connrefused --no-check-certificate -nv -N -c https://github.com/TJU-CMC-Org/MINTmap/raw/master/OtherAnnotations.MINTmap_v1.txt
39
39
# mirdeep2
40
- wget --random-wait --retry-connrefused --no-check-certificate -nv -N - c -O Rfam_for_miRDeep.fa.gz https://github.com/lpantano/mirdeep2_core/raw/data/Rfam_for_miRDeep.fa.gz && gunzip Rfam_for_miRDeep.fa.gz
40
+ wget --random-wait --retry-connrefused --no-check-certificate -nv -c -O Rfam_for_miRDeep.fa.gz https://github.com/lpantano/mirdeep2_core/raw/data/Rfam_for_miRDeep.fa.gz && gunzip Rfam_for_miRDeep.fa.gz
41
41
# targetscan analysis
42
- wget --random-wait --retry-connrefused --no-check-certificate -nv -N - c -O Summary_Counts.txt.zip http://www.targetscan.org/vert_71/vert_71_data_download/Summary_Counts.all_predictions.txt.zip && unzip Summary_Counts.txt.zip
43
- wget --random-wait --retry-connrefused --no-check-certificate -nv -N - c -O miR_Family_Info.txt.zip http://www.targetscan.org/vert_71/vert_71_data_download/miR_Family_Info.txt.zip && unzip miR_Family_Info.txt.zip
42
+ wget --random-wait --retry-connrefused --no-check-certificate -nv -c -O Summary_Counts.txt.zip http://www.targetscan.org/vert_71/vert_71_data_download/Summary_Counts.all_predictions.txt.zip && unzip Summary_Counts.txt.zip
43
+ wget --random-wait --retry-connrefused --no-check-certificate -nv -c -O miR_Family_Info.txt.zip http://www.targetscan.org/vert_71/vert_71_data_download/miR_Family_Info.txt.zip && unzip miR_Family_Info.txt.zip
44
44
wget --random-wait --retry-connrefused --no-check-certificate -nv -N -c ftp://mirbase.org/pub/mirbase/21/database_files/mirna_mature.txt.gz
45
45
recipe_outfiles :
46
46
- srnaseq/srna-transcripts.gtf
0 commit comments