File: create.sh

package info (click to toggle)
snpeff 5.4.b%2Bdfsg-2
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 757,500 kB
  • sloc: java: 62,572; perl: 2,279; sh: 1,184; python: 744; xml: 507; makefile: 48
file content (26 lines) | stat: -rwxr-xr-x 626 bytes parent folder | download | duplicates (3)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
#!/bin/sh

# wget ftp://ftp.solgenomics.net/genomes/Solanum_lycopersicum/wgs/assembly/build_2.40/S_lycopersicum_chromosomes.2.40.fa.gz

# GFF has to be downloaded manually from 
#
# http://solgenomics.net/itag/release/2.3/list_files#
#          File ITAG2.3_gene_models.gff3
#
#

# Create genes.gff
echo Decompressing the file
cp ITAG2.3_gene_models.gff3.gz genes.gff.gz
rm -f genes.gff
gunzip genes.gff.gz 

echo Fixing start-end problem
cat genes.gff | ./fixStartEnd.pl > g
mv g genes.gff

# Append FASTA Sequences
echo "###" >> genes.gff
echo "##FASTA" >> genes.gff
zcat S_lycopersicum_chromosomes.2.40.fa.gz >> genes.gff