How to remove those sequence with same amino acid?What command line I should type?


 
Thread Tools Search this Thread
Top Forums Shell Programming and Scripting How to remove those sequence with same amino acid?What command line I should type?
# 1  
Old 01-20-2009
How to remove those sequence with same amino acid?What command line I should type?

My input is listed as:
giNumber RefAminoAcid VarAminoAcid
10190711 P P
10190711 D D
109255248 I A
110349771 A D

My desired output is:
giNumber RefAminoAcid VarAminoAcid
109255248 I A
110349771 A D

*Those with same amino acid, I want delete it and just remain those different amino acid one at the end.
What command line I should type?
Thanks you and appreciate your advise.
# 2  
Old 01-20-2009
Try this:

Code:
awk 'NR==FNR{a[$1]++;next}a[$1]==1' file file

Regards
# 3  
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..
# 4  
Old 01-20-2009
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
# 5  
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
Login or Register to Ask a Question

Previous Thread | Next Thread

10 More Discussions You Might Find Interesting

1. 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

2. 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

3. 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

4. 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

5. 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

6. 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

7. 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

8. 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

9. 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

10. 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
Login or Register to Ask a Question