BP_PROCESS_WORMBASE(1p) User Contributed Perl Documentation BP_PROCESS_WORMBASE(1p)NAME
process_wormbase.pl - Massage WormBase GFF files into a version suitable for the Generic Genome Browser
SYNOPSIS
% process_wormbase.pl ./WS61 > wormbase.gff
DESCRIPTION
This script massages the Wormbase GFF files located at ftp://www.wormbase.org/pub/wormbase/GENE_DUMPS into a version of the GFF format
suitable for display by the generic genome browser. It mainly adds comments to the annotations and designates certain well-spaced genetic
loci as framework landmarks.
This script requires the AcePerl distribution, which is available on CPAN (look for the "Ace" module).
To use this script, get the WormBase GFF files from the FTP site listed above and place them in a directory. It might be a good idea to
name the directory after the current release, such as WS61. You do not need to uncompress the files.
Then give that directory as the argument to this script and capture the script's output to a file:
% process_wormbase.pl ./WS61 > wormbase.gff
It may take a while before you see output from this script, since it must first fetch gene and protein database from the remote AceDB
running at www.wormbase.org. The wormbase.gff file can then be loaded into a Bio::DB::GFF database using the following command:
% bulk_load_gff.pl -d <databasename> wormbase.gff
SEE ALSO
Bio::DB::GFF, bulk_load_gff.pl, load_gff.pl
AUTHOR
Lincoln Stein <lstein@cshl.org>
Copyright (c) 2002 Cold Spring Harbor Laboratory
This library is free software; you can redistribute it and/or modify it under the same terms as Perl itself. See DISCLAIMER.txt for
disclaimers of warranty.
perl v5.14.2 2012-03-02 BP_PROCESS_WORMBASE(1p)
Check Out this Related Man Page
BP_GENERATE_HISTOGRAM(1p) User Contributed Perl Documentation BP_GENERATE_HISTOGRAM(1p)NAME
generate_histogram.pl -- Generate a histogram of Bio::DB::GFF features
SYNOPSIS
generate_histogram.pl -d gadfly variation gene:curated
DESCRIPTION
Use this utility to generate feature density histograms from Bio::DB::GFF databases. The result is a GFF data file that is suitable for
loading with load_gff.pl.
OPTIONS
The following options are recognized:
Option Description
-----------------
--dsn <dsn> Data source (default dbi:mysql:test)
--adaptor <adaptor> Schema adaptor (default dbi::mysqlopt)
--user <user> Username for mysql authentication
--pass <password> Password for mysql authentication
--aggregator <list> Comma-separated list of aggregators
BUGS
Please report them.
SEE ALSO
Bio::DB::GFF
AUTHOR
Lincoln Stein <lstein@cshl.org>
Copyright (c) 2001 Cold Spring Harbor Laboratory
This library is free software; you can redistribute it and/or modify it under the same terms as Perl itself. See DISCLAIMER.txt for
disclaimers of warranty.
perl v5.14.2 2012-03-02 BP_GENERATE_HISTOGRAM(1p)
Hello,
I need help writing a script that will use the touch -t in
ddmmyyyy HH:mm format.
Basically I have to go into a directory and give all of the files their original time stamp, hence the touch -t command. But where to start ?
I was planning on using this:
touch -t
pwd (as... (2 Replies)
hi,
I have a script which export a database evernight in a directory.
The script alos ftpies the files in the directory to another machine.
Problem is I am copying/ftping everything in the directory. This directory will grow overtime so I want to know I ftp last export in directory and... (5 Replies)
Hi UNIX gurus,
the scenario is that i have written a script which takes as input a directory structure.
And after that a tar is made ,then zipped and FTP it to a server.But how do i uncompress it from my script as FTP doesnot support any uncompress command during FTP session.
I have to automate... (2 Replies)
I'm new to scripting. I'm trying to write a script to download files from ftp site, the following is the script and the message i get after running the script. no files were downloaded :(
Thanks advance!
script:
#!/usr/bin/ksh
DAY=`date --date="-1 days" +%y%m%d`
ftp -v -n "ftp.address" <<... (5 Replies)
I have 20 or so files that are ftp'd each day to a particular directory. I want to create a script that will tell me if any of these files fails to arrive or if any additional files arrive. I'm thinking I'd have a list of file names that should arrive each day, and the script would check each days... (2 Replies)
Hi there. i'm new user at here
I need help for this. I need to write a script that maintains a database of cryptographic checksums. In other words, I need this script to check the files checksum whether the files has been modified or not. But i got no idea where to start. Hope anyone here can... (8 Replies)
Hi All,
Using this command wc -l *e* > create2.txt i'm getting the following output >>>create2.txt listed output files, my requirement is how to cut the first coloum in all the files mentioned in create2.txt in perl.
50 allignment.pl
3 create.txt
4 application.txt
... (8 Replies)
I would like to write a script that would create 1000 files like clm0001.txt to clm1000.txt. All the files would contain the same contents. How do I achieve this using a script?Please help. (2 Replies)
Hi, my problem is that I have two files. File no. 1 is a gff text file (say gi1) that has gene information like :
********************
gene 39389788..39395643
/gene="RPSA"
/note="Derived by automated computational analysis using
... (2 Replies)
Hi all.....
I have a file which contains large data...like
I want to print the rows starting from "pixel" till the file read the letter "TER" into a new output file....
can anyone plz help in doing this ?? (5 Replies)
Hello,
if this input:
gff art ex 2833 2966 . + . ID=A_172736
gff art ex 2976 3165 . + . ID=A_172736
gff art ex 3195 3941 . + . ID=A_713953
gff art ex 8615 8753 . + . ... (5 Replies)
Hi,
I have 1000 files names data1.txt through data1000.txt inside a folder. I want to write a script that will take each first line from the files and write them as output into a new file. How do I go about doing that? Thanks! (2 Replies)
I am writing a shell script to uncompress files in a directory, then call a Perl script to search the files for given terms, store those terms in a different output file , and compress the output. I get a syntax error with my use of foreach. Below is my script.
#!/bin/csh -fxv
if (!... (2 Replies)
Hello everyone,
I'm am a newbie to coding so I am reaching out in hopes that I can get some help from this forum.
I am trying to run the script below from a single directory, however the directory has many subfolders. In each of those subfolders is a file, uniquely named to that specific... (3 Replies)
All,
I am trying to create a report on the duration of an ETL load from the file arrival to the final dump in to a database for SLA's.
Does anyone have any guidance or ideas on how metadata can be extracted; information of a file: like file name, created timestamp, count of records and load... (1 Reply)