Sponsored Content
Top Forums Shell Programming and Scripting How to remove those sequence with same amino acid?What command line I should type? Post 302278644 by patrick chia on Tuesday 20th of January 2009 08:50:18 PM
Old 01-20-2009
Quote:
Originally Posted by Perderabo
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
Code:
$
$ cat file
giNumber RefAminoAcid VarAminoAcid
10190711 P P
10190711 D D
109255248 I A
110349771 A D
$
$
$  awk '$2 != $3' file
giNumber RefAminoAcid VarAminoAcid
109255248 I A
110349771 A D
$
$

The Single Letter Amino Acid Code
Perderabo, thanks a lot for your help. Your method is faster and effective. Really help me a lot. Smilie
 

10 More Discussions You Might Find Interesting

1. Solaris

How to remove user from a group using command line

Hi, Wanted to remove a user from a group , but no GUI , must use command line (2 Replies)
Discussion started by: civic2005
2 Replies

2. Programming

Need help in storing command line argument argv[2] to a variable of int type

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)
Discussion started by: frozensmilz
3 Replies

3. Shell Programming and Scripting

How can I calculate the total of nucleotide in Unix?What command line I should type?

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)
Discussion started by: patrick chia
2 Replies

4. Shell Programming and Scripting

How can I remove those duplicate sequence in UNIX?What command line I should type?

The input is: >HWI-EAS382_30FC7AAXX:4:1:1580:1465 AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA >HWI-EAS382_30FC7AAXX:4:1:1062:1640 AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA >HWI-EAS382_30FC7AAXX:4:1:272:629 AAAAAAAAGCTATAGTCTCGTCACACATACTCACAA >HWI-EAS382_30FC7AAXX:4:1:1033:1135... (4 Replies)
Discussion started by: patrick chia
4 Replies

5. UNIX for Dummies Questions & Answers

Command to remove First and Last line from a File

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)
Discussion started by: ragz_82
1 Replies

6. Shell Programming and Scripting

Remove Command-Line Option from String

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)
Discussion started by: mattmiller
7 Replies

7. Shell Programming and Scripting

How to remove unwanted elements from the command line?

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)
Discussion started by: jediwannabe
7 Replies

8. AIX

Finding the type of AIX server using command line

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)
Discussion started by: omkar.jadhav
9 Replies

9. Shell Programming and Scripting

Convert a DNA sequence into Amino Acid

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)
Discussion started by: faizlo
3 Replies

10. UNIX for Beginners Questions & Answers

Need a command to remove the last word in the first line of a file

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
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)
All times are GMT -4. The time now is 08:28 AM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy