10 More Discussions You Might Find Interesting
1. Solaris
Good morning everyone,
I'm looking for some help to retrieve data in a scenario where I might have made a big mistake. I'm hoping to understand what I did wrong.
My system is made of two Solaris 11 Express servers (old free version for evaluation). The first if for data and the second is... (7 Replies)
Discussion started by: rnd
7 Replies
2. UNIX for Advanced & Expert Users
I have the following data set about the snps ID txt file
POS ID
78599583 rs987435
33395779 rs345783
189807684 rs955894
33907909 rs6088791
75664046 rs11180435
218890658 rs17571465
127630276 rs17011450
90919465 rs6919430
and a gene... (7 Replies)
Discussion started by: marwah
7 Replies
3. Shell Programming and Scripting
I am trying to modify the awk below to include the gene name ($5) for each target and can not seem to do so. Also, I'm not sure the calculation is right (average of all targets that are the same is $4 using the values in $7)? Thank you :).
awk '{if((NR>1)&&($4!=last)){printf("%s\t%f\t%s\n",... (1 Reply)
Discussion started by: cmccabe
1 Replies
4. UNIX for Dummies Questions & Answers
Hello, I'm currently implementing the && function in a shell using C. For example, if we input cmd1 && cmd2, then cmd2 executes only when cmd1 exits successfully. I'm thinking about:
int main() {
int i;
char **args;
while(1) {
printf("yongfeng's shell:~$ ");
args =... (5 Replies)
Discussion started by: Yongfeng
5 Replies
5. Shell Programming and Scripting
Hi ,
I am trying to write a script which does the wincvs login and does some checkout operations.
When i do it from command prompt with command
wincvs -d :pserver:username@server:/cvs/repository
it pops up a cvs window and i can do some operations there. I wanted to automate this through... (1 Reply)
Discussion started by: namishtiwari
1 Replies
6. Shell Programming and Scripting
Hi all,
Using Perl, I need to extract DNA bases from a GenBank file for a given plant species. A sample GenBank file is here...
Nucleotide
This is saved on my computer as NC_001666.gb. I also have a file that is saved on my computer as NC_001666.txt. This text file has a list of all... (5 Replies)
Discussion started by: akreibich07
5 Replies
7. Shell Programming and Scripting
Hey guys,
I'm doing some Perl scripting for genomic data out of GenBank files...I have to extract the name of the plant, the file name, the number of bases, and all of the genes including their starting and ending positions...for example, with this GenBank file,
LOCUS NC_010093 ... (7 Replies)
Discussion started by: akreibich07
7 Replies
8. UNIX for Dummies Questions & Answers
HI GUrus,
I am very mcuh new for this forum ....i need you all guys help for resolving my queries.
As i have DNS server , platform is sunsolaris----8. We are using that serevr for resolving all network nodes.
I am going to explain the process what i am following ....
I could have... (2 Replies)
Discussion started by: leenagoyal
2 Replies
9. Shell Programming and Scripting
hi
i have a function.
parse()
{
qname="lpdma111"
qmgr="lpdma111"
qalias="lpqqali"
}
Now the requirement is i want ot get the vaues form the function
start()
{
//retireve the values from the parse function
} (5 Replies)
Discussion started by: Satyak
5 Replies
10. UNIX for Dummies Questions & Answers
Hello all,
Can someone explain to me the advantage between using subshell over a function call in scripts? To me these are the same. Am I wrong to think this? (4 Replies)
Discussion started by: larry
4 Replies
Bio::SeqIO::entrezgene(3pm) User Contributed Perl Documentation Bio::SeqIO::entrezgene(3pm)
NAME
Bio::SeqIO::entrezgene - Entrez Gene ASN1 parser
SYNOPSIS
use Bio::SeqIO;
# don't instantiate directly - instead do
my $seqio = Bio::SeqIO->new(-format => 'entrezgene',
-file => $file);
my $gene = $seqio->next_seq;
DESCRIPTION
This is EntrezGene ASN bioperl parser. It is built on top of Bio::ASN1::EntrezGene, a low level ASN parser built by Mingyi Liu
(<http://sourceforge.net/projects/egparser>). The easiest way to use it is shown above.
You will get most of the Entrez Gene annotation such as gene symbol, gene name and description, accession numbers associated with the gene,
etc. Almost all of these are given as Bio::AnnotationI objects.
If you need all the data do:
my $seqio = Bio::SeqIO->new(-format => 'entrezgene',
-file => $file,
-debug => 'on' );
my ($gene,$genestructure,$uncaptured) = $seqio->next_seq;
The second variable returned, $genestructure, is a Bio::Cluster::SequenceFamily object. It contains all Refseqs and the genomic contigs
that are associated with the particular gene. The third variable, $uncaptured, is a reference to a plain array.
You can also modify the output to allow back compatibility with the old LocusLink parser:
my $seqio = Bio::SeqIO->new(-format => 'entrezgene',
-file => $file,
-locuslink => 'convert');
The "-debug" and "-locuslink" options slow down the parser.
Example code which looks for ontology terms:
my $eio = new Bio::SeqIO(-file => $file,
-format => 'entrezgene',
-service_record => 'yes');
while (my $seq = $eio->next_seq) {
my $gid = $seq->accession_number;
foreach my $ot ($ann->get_Annotations('OntologyTerm')) {
next if ($ot->term->authority eq 'STS marker'); # No STS markers
my $evid = $ot->comment;
$evid =~ s/evidence: //i;
my @ref = $ot->term->get_references;
my $id = $ot->identifier;
my $fid = 'GO:' . sprintf("%07u",$id);
print join(" ",$gid, $ot->ontology->name, $ot->name, $evid,
$fid, @ref?$ref[0]->medline:''), "
";
}
}
FEEDBACK
Mailing Lists
User feedback is an integral part of the evolution of this and other Bioperl modules. Send your comments and suggestions preferably to the
Bioperl mailing list. Your participation is much appreciated.
bioperl-l@bioperl.org - General discussion
http://bioperl.org/wiki/Mailing_lists - About the mailing lists
Support
Please direct usage questions or support issues to the mailing list:
bioperl-l@bioperl.org
rather than to the module maintainer directly. Many experienced and reponsive experts will be able look at the problem and quickly address
it. Please include a thorough description of the problem with code and data examples if at all possible.
Reporting Bugs
Report bugs to the Bioperl bug tracking system to help us keep track of the bugs and their resolution. Bug reports can be submitted via the
web:
https://redmine.open-bio.org/projects/bioperl/
AUTHOR - Stefan Kirov
Email skirov at utk.edu
CONTRIBUTORS
Hilmar Lapp, hlapp at gmx.net
APPENDIX
This parser is based on Bio::ASN1::EntrezGene module.
The rest of the documentation details each of the object methods. Internal methods are usually preceded with a _
_process_products_coordinates
To do:
_process_prop
To do: process GO
perl v5.14.2 2012-03-02 Bio::SeqIO::entrezgene(3pm)