Sponsored Content
Top Forums Shell Programming and Scripting How to remove those sequence with same amino acid?What command line I should type? Post 302278432 by koti_rama on Tuesday 20th of January 2009 08:07:06 AM
Old 01-20-2009
try like this:
s1-> this contains you source date
out-> out put file


while read line
do
acid=`echo $line | awk ' FS=" " {print $1}'`
count=`grep $acid s1 | wc -l`
if [ "$count" -gt "1" ]; then
echo " Acid $acid exit in out file"
continue;
else
echo "$line" >>out.log
fi
done<s1

Last edited by koti_rama; 01-20-2009 at 09:16 AM..
 

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::Phenotype::OMIM::OMIMentryAllelicVariant(3pm)	User Contributed Perl Documentation	Bio::Phenotype::OMIM::OMIMentryAllelicVariant(3pm)

NAME
Bio::Phenotype::OMIM::OMIMentryAllelicVariant - Representation of a allelic variant of the OMIM database SYNOPSIS
use Bio::Phenotype::OMIM::OMIMentryAllelicVariant; $av = Bio::Phenotype::OMIM::OMIMentryAllelicVariant->new( -number => ".0001", -title => "ALCOHOL INTOLERANCE", -symbol => "ALDH2*2", -description => "The ALDH2*2-encoded ...", -aa_ori => "GLU", -aa_mut => "LYS", -position => 487, -additional_mutations => "IVS4DS, G-A, +1" ); DESCRIPTION
This class models the allelic variant of the OMIM database. This class is intended to be used together with a OMIM entry class. 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 one of 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
Christian M. Zmasek Email: czmasek-at-burnham.org or cmzmasek@yahoo.com WWW: http://monochrome-effect.net/ Address: Genomics Institute of the Novartis Research Foundation 10675 John Jay Hopkins Drive San Diego, CA 92121 APPENDIX
The rest of the documentation details each of the object methods. new Title : new Usage : $av = Bio::Phenotype::OMIM::OMIMentryAllelicVariant->new( -number => ".0001", -title => "ALCOHOL INTOLERANCE", -symbol => "ALDH2*2", -description => "The ALDH2*2-encoded ...", -aa_ori => "GLU", -aa_mut => "LYS", -position => 487, -additional_mutations => "IVS4DS, G-A, +1" ); Function: Creates a new OMIMentryAllelicVariant object. Returns : A new OMIMentryAllelicVariant object. Args : -number => the OMIM allelic variant number -title => the title -symbol => a symbol -description => a description -aa_ori => the original amino acid -aa_mut => the mutated amino acid -position => the position of the mutation -additional_mutations => free form description of additional mutations init Title : init() Usage : $av->init(); Function: Initializes this OMIMentryAllelicVariant to all "". Returns : Args : number Title : number Usage : $av->number( ".0001" ); or print $av->number(); Function: Set/get for the OMIM allelic variant number of this OMIMentryAllelicVariant. Returns : The OMIM allelic variant number. Args : The OMIM allelic variant number (optional). title Title : title Usage : $av->title( "ALCOHOL INTOLERANCE" ); or print $av->title(); Function: Set/get for the title of this OMIMentryAllelicVariant. Returns : The title. Args : The title (optional). symbol Title : symbol Usage : $av->symbol( "ALDH2*2" ); or print $av->symbol(); Function: Set/get for the symbol of this OMIMentryAllelicVariant. Returns : A symbol. Args : A symbol (optional). description Title : description Usage : $av->description( "The ALDH2*2-encoded protein has a change ..." ); or print $av->description(); Function: Set/get for the description of this OMIMentryAllelicVariant. Returns : A description. Args : A description (optional). aa_ori Title : aa_ori Usage : $av->aa_ori( "GLU" ); or print $av->aa_ori(); Function: Set/get for the original amino acid(s). Returns : The original amino acid(s). Args : The original amino acid(s) (optional). aa_mut Title : aa_mut Usage : $av->aa_mut( "LYS" ); or print $av->aa_mut(); Function: Set/get for the mutated amino acid(s). Returns : The mutated amino acid(s). Args : The mutated amino acid(s) (optional). position Title : position Usage : $av->position( 487 ); or print $av->position(); Function: Set/get for the position of the mutation. Returns : The position of the mutation. Args : The position of the mutation (optional). additional_mutations Title : additional_mutations Usage : $av->additional_mutations( "1-BP DEL, 911T" ); or print $av->additional_mutations(); Function: Set/get for free form description of (additional) mutation(s). Returns : description of (additional) mutation(s). Args : description of (additional) mutation(s) (optional). to_string Title : to_string() Usage : print $av->to_string(); Function: To string method for OMIMentryAllelicVariant objects. Returns : A string representations of this OMIMentryAllelicVariant. Args : perl v5.14.2 2012-03-02 Bio::Phenotype::OMIM::OMIMentryAllelicVariant(3pm)
All times are GMT -4. The time now is 04:26 AM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy