Sponsored Content
Top Forums Shell Programming and Scripting Match and retrieve information from file Post 302944275 by Wan Fahmi on Monday 18th of May 2015 11:12:35 AM
Old 05-18-2015
Match and retrieve information from file

Hello

I just want to ask how to get the match of information column 2 file 1 and retrieve information from column 2 file 2. The column exon in file 1 and column color code in file 2.

File 1
Code:
//NODECOLORCODE         "Exon 1"        "ENST00000595813"
//NODECOLORCODE         "Exon 1"        "ENST00000596051"
//NODECOLORCODE         "Exon 2"        "ENST00000243643"
//NODECOLORCODE         "Exon 2"        "ENST00000421033"
//NODECOLORCODE         "Exon 3"        "ENST00000243643"
//NODECOLORCODE         "Exon 3"        "ENST00000421033"
//NODECOLORCODE         "Exon 4"        "ENST00000243643"
//NODECOLORCODE         "Exon 4"        "ENST00000421033"
//NODECOLORCODE         "Exon 5"        "ENST00000421033"
//NODECOLORCODE         "Exon 5"        "ENST00000440291"

File 2
Code:
"Exon 1"        "#FF0000"
"Exon 2"        "#FFD400"
"Exon 3"        "#FFff00"
"Exon 4"        "#2Bff00"
"Exon 5"        "#00ff00"

Desired output

Code:
//NODECOLORCODE         "Exon 1"        "ENST00000595813"        "#FF0000"
//NODECOLORCODE         "Exon 1"        "ENST00000596051"        "#FF0000"
//NODECOLORCODE         "Exon 2"        "ENST00000243643"        "#FFD400"
//NODECOLORCODE         "Exon 2"        "ENST00000421033"        "#FFD400"
//NODECOLORCODE         "Exon 3"        "ENST00000243643"        "#FFff00"
//NODECOLORCODE         "Exon 3"        "ENST00000421033"        "#FFff00"
//NODECOLORCODE         "Exon 4"        "ENST00000243643"        "#2Bff00"
//NODECOLORCODE         "Exon 4"        "ENST00000421033"        "#2Bff00"
//NODECOLORCODE         "Exon 5"        "ENST00000421033"        "#00ff00"
//NODECOLORCODE         "Exon 5"        "ENST00000440291"        "#00ff00"

I try using awk and it doesn't come out with the desired output I want. Thanks

Code:
awk 'FNR==NR {C[$2]=$1;next} FNR==1 {print "NODECLASSCOLOR EXON TRANSCRIPT_ID COLORCODE"; next} $1 in C {print C[$1], $1, $4, $2, $3}'

 

10 More Discussions You Might Find Interesting

1. Shell Programming and Scripting

A script pls( To retrieve database information)

KSH - Unix -AIX - db2 ************** Get the input from the user say '123' (It varies) Then i want to connect to a database say "D1" then i want th extract from the database by giving "select * from tablename where input = '123' I also want to connect to another database "D2" then i... (3 Replies)
Discussion started by: rollthecoin
3 Replies

2. Solaris

command to retrieve user information

Hi, I want the command to retrieve the existing user information such as * authorization * Profile * role * exipre(expiration date of login) * inactive please tell me how to do that Thank you. (3 Replies)
Discussion started by: S_venkatesh
3 Replies

3. UNIX for Dummies Questions & Answers

retrieve lines that match a pattern

Hi, I would like to know how can I get lines from a text file that match no more than 2 '>'. Example: Input file: a >cr1 4 a>b b>c a >cr2 5 a>b Output file: a >cr2 5 a>b Thanks in advance (2 Replies)
Discussion started by: fadista
2 Replies

4. Shell Programming and Scripting

Execute a C program and retrieve information

Hi I have the following script: #!/bin/sh gcc -o program program.c ./program & PID=$! where i execute a C program and i get its pid. I want to retrieve information about this program (e.g memory consumption) using command top. So far i have: top -d 1.0 -p $PID But i dont know how to... (6 Replies)
Discussion started by: nteath
6 Replies

5. Shell Programming and Scripting

Retrieve lines that match any occurence in a list of patterns

I have two files. The first containing a header and six columns of data. Example file 1: Number SNP ID dbSNP RS ID Chromosome Result_Call Physical Position 787066 SNP_A-8575395 RS6650104 1 NOCALL 564477 786872 SNP_A-8575125 RS10458597 1 AA ... (13 Replies)
Discussion started by: Selftaught
13 Replies

6. Shell Programming and Scripting

Match pattern1 in file, match pattern2, substitute value1 in line

not getting anywhere with this an xml file contains multiple clients set up with same tags, different values. I need to parse the file for client foo, and change the value of tag "64bit" from false to true. cat clients.xml <Client type"FIX"> <ClientName>foo</ClientName>... (3 Replies)
Discussion started by: jack.bauer
3 Replies

7. Shell Programming and Scripting

Retrieve information Text/Word from HTML code using awk/sed

awk/sed newbie here. I have a HTML file and from that file and I would like to retrieve a text word. <font face=arial size=-1><li><a href=/value_for_clients/Tokyo/abc_process.txt>abc</a> NDK Version: 4.0 </li> <font face=arial size=-1><li><a... (6 Replies)
Discussion started by: sk2code
6 Replies

8. Shell Programming and Scripting

Display match or no match and write a text file to a directory

The below bash connects to a site, downloads a file, searches that file based of user input - could be multiple (all that seems to work). What I am not able to figure out is how to display on the screen match found or no match found" and write a file to a directory (C:\Users\cmccabe\Desktop\wget)... (4 Replies)
Discussion started by: cmccabe
4 Replies

9. Solaris

Solaris commands to retrieve chipset information

I need to know what are the commands in Solaris to retrieve the below information about the hardware platform. 1. Chipset information (information about various hardware controller cards on the mother boards, system BIOS versions, PCI firmware version etc..) 2. Serial number of the work... (2 Replies)
Discussion started by: rajujayanthy
2 Replies

10. Shell Programming and Scripting

awk to update file based on partial match in field1 and exact match in field2

I am trying to create a cronjob that will run on startup that will look at a list.txt file to see if there is a later version of a database using database.txt as the source. The matching lines are written to output. $1 in database.txt will be in list.txt as a partial match. $2 of database.txt... (2 Replies)
Discussion started by: cmccabe
2 Replies
Bio::Tools::Sim4::Exon(3pm)				User Contributed Perl Documentation			       Bio::Tools::Sim4::Exon(3pm)

NAME
Bio::Tools::Sim4::Exon - A single exon determined by an alignment SYNOPSIS
# See Bio::Tools::Sim4::Results for a description of the context. # an instance of this class is-a Bio::SeqFeature::SimilarityPair # coordinates of the exon (recommended way): print "exon from ", $exon->start(), " to ", $exon->end(), " "; # the same (feature1() inherited from Bio::SeqFeature::FeaturePair) print "exon from ", $exon->feature1()->start(), " to ", $exon->feature1()->end(), " "; # also the same (query() inherited from Bio::SeqFeature::SimilarityPair): print "exon from ", $exon->query()->start(), " to ", $exon->query()->end(), " "; # coordinates on the matching EST (recommended way): print "matches on EST from ", $exon->est_hit()->start(), " to ", $exon->est_hit()->end(), " "; # the same (feature2() inherited from Bio::SeqFeature::FeaturePair) print "matches on EST from ", $exon->feature2()->start(), " to ", $exon->feature2()->end(), " "; # also the same (subject() inherited from Bio::SeqFeature::SimilarityPair): print "exon from ", $exon->subject()->start(), " to ", $exon->subject()->end(), " "; DESCRIPTION
This class inherits from Bio::SeqFeature::SimilarityPair and represents an exon on a genomic sequence determined by similarity, that is, by aligning an EST sequence (using Sim4 in this case). Consequently, the notion of query and subject is always from the perspective of the genomic sequence: query refers to the genomic seq, subject to the aligned EST hit. Because of this, $exon->start(), $exon->end() etc will always return what you expect. To get the coordinates on the matching EST, refer to the properties of the feature returned by est_hit(). 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 - Ewan Birney, Hilmar Lapp Ewan Birney <birney-at-sanger.ac.uk> Hilmar Lapp <hlapp-at-gmx.net> or <hilmar.lapp-at-pharma.novartis.com>. APPENDIX
The rest of the documentation details each of the object methods. Internal methods are usually preceded with a _ percentage_id Title : percentage_id Usage : $obj->percentage_id($newval) Function: This is a synonym for 100 * $obj->est_hit()->frac_identical(). Returns : value of percentage_id Args : newvalue (optional) est_hit Title : est_hit Usage : $est_feature = $obj->est_hit(); Function: Returns the EST hit pointing to (i.e., aligned to by Sim4) this exon (i.e., genomic region). At present, merely a synonym for $obj->feature2(). Returns : An Bio::SeqFeatureI implementing object. Args : perl v5.14.2 2012-03-02 Bio::Tools::Sim4::Exon(3pm)
All times are GMT -4. The time now is 11:21 PM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy