Actually that first field is not an amino acid. Field 2 and field 3 must be different from each other. If this understanding is right, then
The Single Letter Amino Acid Code
Perderabo, thanks a lot for your help. Your method is faster and effective. Really help me a lot.
The following program takes two command line arguments.
I want the second argument (fileCount) to be stored/printed as a int value. I tried my best to typecast the char to int (check the printf statement at last) but is not working...the output is some junk value.
This program is in its... (3 Replies)
For example, if I have the file whose content are:
>HWI-EAS382_30FC7AAXX:7:1:927:1368
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
>HWI-EAS382_30FC7AAXX:7:1:924:1373
ACGAACTTTAAAGCACCTCTTGGCTCGTATGCCGTC
I want my output calculate the total of nucleotide. So my output should look like this:... (2 Replies)
I have a file from which the Header and the Trailer lines need to be removed. They are confirmed to be the first and the last lines in the file.
I have tried a few commands, but not successful yet. It needs to be implemented urgently, hence any help is greatly appreciated.
Raghu
----------... (1 Reply)
I want to add a "-r <remote_host>" option to my ksh script, causing the script to run a script of the same name on the specified remote host. The remote invocation should itself include all the command-line options of the original invocation, less the -r option.
For example, this invocation:
... (7 Replies)
CREATE TYPE "SUPERB"."OE_L2_FOOTERDATA_REC"
OID '1232323232' IS OBJECT
(
AWLSTrxId NUMBER,
AdjPurchaseAmount NUMBER,
Network NUMBER,
ConfigVersion NUMBER
);
/
I wanted my output to be as follow:
CREATE TYPE OE_L2_FOOTERDATA_REC
OID '1232323232'... (7 Replies)
Hi All,
I am trying to find out type of AIX servers , for eg : server A is eServer p5 520 , server B is pSeries 630-6C4 etc.
Can some one suggest me the best way to retrive this information through command line.
is there any direct command or script by which i can retrieve this information... (9 Replies)
I am trying to write a bash script that would be able to read DNA sequences (each line in the file is a sequence) from a file, where sequences are separated by an empty line. I am then to find the amino acid that these DNA sequences encode per codon (each group of three literals.) For example, if I... (3 Replies)
I have a eg file op.txt
This is a cat
This is a fat cat
This is a fat black cat
I want to remove only the word cat from the first alone .can somebody help.
please do wrap your samples in CODE TAGS
As per forum rules. (7 Replies)
Discussion started by: Sharks
7 Replies
LEARN ABOUT DEBIAN
bio::variation::aareversemutate
Bio::Variation::AAReverseMutate(3pm) User Contributed Perl Documentation Bio::Variation::AAReverseMutate(3pm)NAME
Bio::Variation::AAReverseMutate - point mutation and codon
information from single amino acid changes
SYNOPSIS
$aamut = Bio::Variation::AAReverseMutate->new
(-aa_ori => 'F',
-aa_mut => 'S',
-codon_ori => 'ttc', # optional
-codon_table => '3' # defaults to 1
);
@points = $aamut->each_Variant;
if (scalar @points > 0 ) {
foreach $rnachange ( @points ) {
# $rnachange is a Bio::Variation::RNAChange object
print " ", $rnachange->allele_ori->seq, ">",
$rnachange->allele_mut->seq, " in ",
$rnachange->codon_ori, ">", $rnachange->codon_mut,
" at position ", $rnachange->codon_pos, "
";
}
} else {
print "No point mutations possible
",
}
DESCRIPTION
Bio::Variation::AAReverseMutate objects take in reference and mutated amino acid information and deduces potential point mutations at RNA
level leading to this change. The choice can be further limited by letting the object know what is the the codon in the reference sequence.
The results are returned as Bio::Variation::RNAChange objects.
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 lists 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 the bugs and their resolution. Bug reports can be submitted via the
web:
https://redmine.open-bio.org/projects/bioperl/
AUTHOR - Heikki Lehvaslaiho
Email: heikki-at-bioperl-dot-org
APPENDIX
The rest of the documentation details each of the object methods. Internal methods are usually preceded with a _
aa_ori
Title : aa_ori
Usage : $obj->aa_ori();
Function:
Sets and returns original aa sequence. If value is not
set, returns false.
Amino acid sequences are stored in upper case characters,
others in lower case.
Example :
Returns : string
Args : single character amino acid code
aa_mut
Title : aa_mut
Usage : $obj->aa_mut();
Function:
Sets and returns the mutated allele sequence. If value is not
set, returns false.
Example :
Returns : string
Args : single character amino acid code
codon_ori
Title : codon_ori
Usage : $obj->codon_ori();
Function:
Sets and returns codon_ori triplet. If value is not set,
returns false. The string has to be three characters
long. The chracter content is not checked.
Example :
Returns : string
Args : string
codon_table
Title : codon_table
Usage : $obj->codon_table();
Function:
Sets and returns the codon table id of the RNA
If value is not set, returns 1, 'universal' code, as the default.
Example :
Returns : integer
Args : none if get, the new value if set
each_Variant
Title : each_Variant
Usage : $obj->each_Variant();
Function:
Returns a list of Variants.
Example :
Returns : list of Variants
Args : none
perl v5.14.2 2012-03-02 Bio::Variation::AAReverseMutate(3pm)