Sponsored Content
Top Forums Shell Programming and Scripting awk to split one field and print the last two fields within the split part. Post 302624773 by yifangt on Monday 16th of April 2012 09:41:57 PM
Old 04-16-2012
awk to split one field and print the last two fields within the split part.

Hello;
I have a file consists of 4 columns separated by tab. The problem is the third fields. Some of the them are very long but can be split by the vertical bar "|". Also some of them do not contain the string "UniProt", but I could ignore it at this moment, and sort the file afterwards. Here is the input file like (>100k rows!!!):
Code:
B_1019927    ES904790    Unknown    Unknown protein
B_1065532    JCVI_27855   JCVI_27855|gb|DW997385|gb|CD815403|gb|EV204160|gb|CD813268|gb|EX038736|UniProt|Q8RWM2    Putative uncharacterized protein At4g36940
B_1065533    JCVI_27856   JCVI_27856|gb|CD826551|gb|DY030639|gb|CX272767|gb|EV056662|gb|EX022554|gb|EX038271|UniProt|Q3E9T4    Protein At4g31430
B_1065534    JCVI_27857   JCVI_27857|gb|DY023794|gb|EE530065|gb|CX280035|UniProt|Q84WF5    Probable VAMP-like protein At1g33475
B_1065535    JCVI_27858   JCVI_27858|gb|ES941909|gb|ES940913|gb|EE514812|UniProt|Q8LBM4    Iron-sulfur assembly protein IscA-like 1, mitochondrial precursor

I want the output file with the truncated field-3, which is like following:
Code:
B_1019927    ES904790    Unknown    Unknown protein
B_1065532    JCVI_27855    UniProt    Q8RWM2     Putative uncharacterized protein At4g36940
B_1065533    JCVI_27856    UniProt    Q3E9T4     Protein At4g31430
B_1065534    JCVI_27857    UniProt    Q84WF5     Probable VAMP-like protein At1g33475
B_1065535    JCVI_27858    Uniprot    Q8LBM4     Iron-sulfur assembly protein IscA-like 1, mitochondrial precursor

What I did is by matching the "UniProt" string in the third field, split this field and only print the last two columns of the split fields, the original 3 fields ($1, $2, $4) need to be printed untouched.

Here is my code:
Code:
awk 'FS="\t" {for(i=1;i<=NF;i++) if(match($i,"UniProt")) split($i, a, "|"); {print $1, "\t", $2, "\t", $a[-2],"\t", $a[-1],"\t",$NF}}' infile.csv

Have struggled for a while, but did not work out. Appreciate your expertise to help me out. Thanks a lot! Yifang
 

10 More Discussions You Might Find Interesting

1. Shell Programming and Scripting

Split a field in awk script

Hi all, I have a field in the line, let's say argument $6, which is in the format 00.00 If i want to split the field to get rid of the "." in between of the amount, how can i do that i awk script? I have it like this split($6,a,".") but it will get rid of the last 2 digits after the... (4 Replies)
Discussion started by: CamTu
4 Replies

2. Shell Programming and Scripting

can Awk split my field on the .

I am trying to pull data out of a file to execute a series of repetitive commands. One of the fields is a {hex number}.{hex number} pattern. I need the first part of that structure to pass into the command. The file structure is like this: A 2007/10/20 09:50:00 li text ef6da.1 name ... (6 Replies)
Discussion started by: oly_r
6 Replies

3. Shell Programming and Scripting

How to split a field into two fields?

Hi, I have a comma delimited text file where character fields (as opposed to numeric and date fields) are always enclosed with double quotes. Records are separated by the newline character. In a shell script I would like to split a particular field into two separate fields (enclosed with double... (4 Replies)
Discussion started by: vbrown
4 Replies

4. Shell Programming and Scripting

Split field with awk

Hi there I have a text file with several fields separated by ";" as follow : 5466-43;5466;JOAN;LIGA;LIGA ESPANOLA;43;DEP LA CORUNA - OSASUNA;10/01/10 17:00 5466-44;5466;CARLES;LIGA;LIGA ESPANOLA;44;MALAGA - ATHL BILBAO;10/01/10 17:00 5466-45;5466;FAB;LIGA;LIGA ESPANOLA;45;REAL MADRID -... (4 Replies)
Discussion started by: capnino
4 Replies

5. Shell Programming and Scripting

AWK:Split fields separated by semicolon

Hi all, I have a .vcf file which contains 8 coulmns and the data under each column as shown below, CHROM POS ID REF ALT QUAL FILTER INFO 1 3000012 . A G 126 ... (6 Replies)
Discussion started by: mehar
6 Replies

6. Shell Programming and Scripting

awk split lines without knowing the number of fields a-priori

I want to use awk to split fields and put them into a file but I don't know the number of fields for example, in the following line Ports: 22/filtered/tcp//ssh///, 53/open/tcp//tcpwrapped///, 111/filtered/tcp//rpcbind///, 543/filtered/tcp//klogin///, 544/filtered/tcp//kshell///,... (3 Replies)
Discussion started by: esolvepolito
3 Replies

7. Shell Programming and Scripting

How to print 1st field and last 2 fields together and the rest of the fields after it using awk?

Hi experts, I need to print the first field first then last two fields should come next and then i need to print rest of the fields. Input : a1,abc,jsd,fhf,fkk,b1,b2 a2,acb,dfg,ghj,b3,c4 a3,djf,wdjg,fkg,dff,ggk,d4,d5 Expected output: a1,b1,b2,abc,jsd,fhf,fkk... (6 Replies)
Discussion started by: 100bees
6 Replies

8. Shell Programming and Scripting

awk split after second underscore in field

I am trying to split a tab-delimeted file using awk after the second _ in bold. The awk below is close but splits on the first _, and I am not sure how to use the second _. Thank you :). file chr1 92145889 92149424 NM_001195684_exon_0_10_chr1_92145900_r 0 - chr1 92161218 ... (4 Replies)
Discussion started by: cmccabe
4 Replies

9. Shell Programming and Scripting

awk to add plus or minus to fields and split another field

In the tab-delimited input below I am trying to use awk to -10 from $2 and +10 to $3. Something like awk -F'\t' -v OFS='\t' -v s=10 '{split($4,a,":"); print $1,$2-s,$3+s,a,$5,$6} | awk {split(a,b,"-"); print $1,$2-s,$3+s,b-s,b+s,$5,$6}' input should do that. I also need to -10 from $4... (2 Replies)
Discussion started by: cmccabe
2 Replies

10. UNIX for Beginners Questions & Answers

awk to split field twice using two deliminators

In the awk I am splitting on the : into array a, then splitting on the - into element b. I can not seem to duplicate b if there is no - after it. Lines 1,2,4 are examples. If there is a - after the number in b then the value to the right of it is $3 in the ouput. Thank you :). awk... (2 Replies)
Discussion started by: cmccabe
2 Replies
Bio::Tools::Analysis::Protein::Scansite(3pm)		User Contributed Perl Documentation	      Bio::Tools::Analysis::Protein::Scansite(3pm)

NAME
Bio::Tools::Analysis::Protein::Scansite - a wrapper around the Scansite server SYNOPSIS
use Bio::Tools::Analysis::Protein::Scansite; my $seq; # a Bio::PrimarySeqI object my $tool = Bio::Tools::Analysis::Protein::Scansite->new ( -seq => $seq->primary_seq ); # run Scansite prediction on a sequence $tool->run(); # alternatively you can say $tool->seq($seq->primary_seq)->run; die "Could not get a result" unless $tool->status =~ /^COMPLETED/; print $tool->result; # print raw prediction to STDOUT foreach my $feat ( $tool->result('Bio::SeqFeatureI') ) { # do something to SeqFeature # e.g. print as GFF print $feat->gff_string, " "; # or store within the sequence - if it is a Bio::RichSeqI $seq->add_SeqFeature($feat); } DESCRIPTION
This class is a wrapper around the Scansite 2.0 server which produces predictions for serine, threonine and tyrosine phosphorylation sites in eukaryotic proteins. At present this is a basic wrapper for the "Scan protein by input sequence" functionality, which takes a sequence and searches for motifs, with the option to select the search stringency. At present, searches for specific phosphorylation sites are not supported; all predicted sites are returned. Return formats The Scansite results can be obtained in several formats: 1. By calling my $res = $tool->result(''); $res holds a string of the predicted sites in tabular format. 2. By calling my $data_ref = $tool->result('value') $data_ref is a reference to an array of hashes. Each element in the array represents a predicted phosphorylation site. The hash keys are the names of the data fields,i.e., 'motif' => 'Casn_Kin1' # name of kinase 'percentile' => 0.155 # see Scansite docs 'position' => 9 # position in protein 'protein' => 'A1' # protein id 'score' => 0.3696 # see Scansite docs 'sequence' => 'ASYFDTASYFSADAT' # sequence surrounding site 'site' => 'S9' # phosphorylated residue 'zscore' => '-3.110' # see Scansite docs 3. By calling my @fts = $tool->Result('Bio::SeqFeatureI'); which returns an array of Bio::SeqFeatureI compliant objects with primary tag value 'Site' and tag names of 'motif', 'score', 'sequence', 'zscore' as above. See <http://scansite.mit.edu/>. This inherits Bio::SimpleAnalysisI which hopefully makes it easier to write wrappers on various services. This class uses a web resource and therefore inherits from Bio::WebAgent. SEE ALSO
Bio::SimpleAnalysisI, Bio::WebAgent 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/ AUTHORS
Richard Adams, Richard.Adams@ed.ac.uk, APPENDIX
The rest of the documentation details each of the object methods. Internal methods are usually preceded with a _ result Name : result Usage : $job->result (...) Returns : a result created by running an analysis Args : none (but an implementation may choose to add arguments for instructions how to process the raw result) The method returns a scalar representing a result of an executed job. If the job was terminated by an error, the result may contain an error message instead of the real data. This implementation returns differently processed data depending on argument: undef Returns the raw ASCII data stream but without HTML tags 'Bio::SeqFeatureI' The argument string defined the type of bioperl objects returned in an array. The objects are Bio::SeqFeature::Generic. 'parsed' Returns a reference to an array of hashes containing the data of one phosphorylation site prediction. Key values are: motif, percentile, position, protein, score, site, zscore, sequence. stringency Usage : $job->stringency(...) Returns : The significance stringency of a prediction Args : None (retrieves value) or 'High', 'Medium' or 'Low'. Purpose : Get/setter of the stringency to be sumitted for analysis. protein_id Usage : $job->protein_id(...) Returns : The sequence id of the protein or 'unnamed' if not set. Args : None Purpose : Getter of the seq_id. Returns the display_id of the sequence object. perl v5.14.2 2012-03-02 Bio::Tools::Analysis::Protein::Scansite(3pm)
All times are GMT -4. The time now is 03:51 AM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy