Printing into two files under difference situation
I want to printing into two files under difference situation.
For example, file 1 name.txt
File 2 a.gff
I want to check each line of file a.gff; if the two element"Glyma01g07930" and "Glyma02g13330" are both existing in file 1, print this line into file a_yes.gff; otherwise print the line into a_no.gff.
However, I have bunch of file 2 to check, like a.gff,b.gff,c.gff.....
How can I do that?
Both linux command line and perl script would be appreciated. Thank you.
Last edited by Scott; 10-09-2013 at 01:13 AM..
Reason: Code tags
Well, I hope this way you will respond to my inquiries.
I have 4 unix servers,with static ips (though i dont think this is an issue)....i can telnet and rlogin from one to the other....if i FTP from on et othe other and try to execute :
cd /user return
/user : no such file or... (1 Reply)
Hi,
I'm looking for someone who can think in sed. Basically, I need the trailing characters on every line in a file to be deleted. These characters are all in capitals, and always follow a number, but they often vary in number
For instance, on the line:
2006_10_9_p20_TALK
I'd want to... (4 Replies)
Hi All,
I am writing some data's into a file from C++ program. The files which i am writing is of fixed length . say 232 in length per line.
I am writing as . my c code is as
... (0 Replies)
Hi
I had deleted important files from my company server :(
the server is HPUX
and i don't know how to undo rm command or how to restore the files ..
iam appreciate for any help
Thanx ... (5 Replies)
Gurus,
I have one file which is having multiple columns and also this file is not always contain the exact columns; sometimes it contains 5 columns or 12 columns. Now, I need to find the difference from that particular file. Here is the sample file:
param1 | 10 | 20 | 30 |
param2 | 10 |... (6 Replies)
Hi Experts,
My requirement is to compare the second field/column in two files, if the second column is same in both the files then compare the first field. If the first is not matching then print the first and second fields of both the files.
first file (a .txt)
< 1210018971FF0000,... (6 Replies)
Hello,
I am PhD student (Biomedical sciences) and very new to Linux. I need some help with the following task :
I have files in the following format for their names :
An_A1_nnn_R1.txt; An_A1_nnm_R1.txt; An_A1_nnoo_R1.txt
An_A2_nnn_R1.txt; An_A2_nnm_R1.txt; An_A2_nno_R1.txt
... (8 Replies)
Hi,
I want to compare two files and print out their differences
e.g:
t1.txt
a,b,c,d
t2.txt
a,b,c,d,e,f
Output
e,f
Currently I do this long about way
tr ',' '\n' <t1.txt >t1.tmp
tr ',' '\n' <t2.txt >t2.tmp
diff t1.tmp t2.tmp > t12.tmp
I have to this comparison for 100 files, so... (3 Replies)
I have two file as given below which shows the ACL permissions of each file. I need to compare the source file with target file and list down the difference as specified below in required output. Can someone help me on this ?
Source File
*************
# file: /local/test_1
# owner: own
#... (4 Replies)
Discussion started by: sarathy_a35
4 Replies
LEARN ABOUT DEBIAN
bp_process_wormbase
BP_PROCESS_WORMBASE(1p) User Contributed Perl Documentation BP_PROCESS_WORMBASE(1p)NAME
process_wormbase.pl - Massage WormBase GFF files into a version suitable for the Generic Genome Browser
SYNOPSIS
% process_wormbase.pl ./WS61 > wormbase.gff
DESCRIPTION
This script massages the Wormbase GFF files located at ftp://www.wormbase.org/pub/wormbase/GENE_DUMPS into a version of the GFF format
suitable for display by the generic genome browser. It mainly adds comments to the annotations and designates certain well-spaced genetic
loci as framework landmarks.
This script requires the AcePerl distribution, which is available on CPAN (look for the "Ace" module).
To use this script, get the WormBase GFF files from the FTP site listed above and place them in a directory. It might be a good idea to
name the directory after the current release, such as WS61. You do not need to uncompress the files.
Then give that directory as the argument to this script and capture the script's output to a file:
% process_wormbase.pl ./WS61 > wormbase.gff
It may take a while before you see output from this script, since it must first fetch gene and protein database from the remote AceDB
running at www.wormbase.org. The wormbase.gff file can then be loaded into a Bio::DB::GFF database using the following command:
% bulk_load_gff.pl -d <databasename> wormbase.gff
SEE ALSO
Bio::DB::GFF, bulk_load_gff.pl, load_gff.pl
AUTHOR
Lincoln Stein <lstein@cshl.org>
Copyright (c) 2002 Cold Spring Harbor Laboratory
This library is free software; you can redistribute it and/or modify it under the same terms as Perl itself. See DISCLAIMER.txt for
disclaimers of warranty.
perl v5.14.2 2012-03-02 BP_PROCESS_WORMBASE(1p)