01-14-2008
try zcat , if you have it.
10 More Discussions You Might Find Interesting
1. UNIX for Dummies Questions & Answers
Give shell script....which takes two file names as input and compares the contents, is both are same delete second file's contents.....
I try with "diff"...... but confusion how to use "diff" with if ---else
Thanking you (5 Replies)
Discussion started by: krishnampkkm
5 Replies
2. Shell Programming and Scripting
Hi,
I am new to shell programming and need help. I have File1 with some ID numbers and File2 with ID number and some associated information.
I want to match the ID numbers from File1 to contents in File2 and output a third file which pulls out the ID numbers and the associated information with... (2 Replies)
Discussion started by: newpro
2 Replies
3. Shell Programming and Scripting
Hi,
I want to delete the contents of a file which are matching with contents of other file in shell scripting.
Ex.
file1
sheel,sumit,1,2,3,4,5,6,7,8
sumit,rana,2,3,4,5,6,7,8,9
grade,pass,2,3,4,5,6,232,1,1
name,sur,33,1,4,12,3,5,6,8
sheel,pass,2,3,4,5,6,232,1,1
File2... (3 Replies)
Discussion started by: ranasheel2000
3 Replies
4. Shell Programming and Scripting
Hi,
I am facing issue while reading data from a file in UNIX. my requirement is to compare two files and for the text pattern matching in the 1st file, replace the contents in second file by the contents of first file from start to the end and write the contents to thrid file.
i am able to... (2 Replies)
Discussion started by: seeki
2 Replies
5. Shell Programming and Scripting
Hi one of the output of the command is as below
# sed -n "/CCM-ResourceHealthCheck:/,/---------/{/CCM-ResourceHealthCheck:/d;/---------/d;p;}" Automation.OutputZ$zoneCounter | sed 's/$/<br>/'
Resource List : <br>
*************************** 1. row ***************************<br>
... (2 Replies)
Discussion started by: vivek d r
2 Replies
6. UNIX for Advanced & Expert Users
When I use this command:
vi /home/bob/.vimI expect to see.
" ============================================================================
" Netrw Directory Listing (netrw v149)
" /home/bob/.vim
" Sorted by name
" Sort sequence:... (4 Replies)
Discussion started by: cokedude
4 Replies
7. UNIX for Beginners Questions & Answers
I have run the following command : od -c Results_May18.fixrank | head
Here is the result. I wanted the results in tab delimited. Thanks
$ od -c Results_May18.fixrank | head
0000000 M 0 1 6 0 1 : 1 2 9 : 0 0 0 0 0
0000020 0 0 0 0 - A T T D Y ... (2 Replies)
Discussion started by: Benard
2 Replies
8. Shell Programming and Scripting
Shell script logic
Hi
I have 2 input files like with file 1 content as (file1)
"BRGTEST-242" a.txt "BRGTEST-240" a.txt "BRGTEST-219" e.txt
File 2 contents as fle(2)
"BRGTEST-244" a.txt "BRGTEST-244" b.txt "BRGTEST-231" c.txt "BRGTEST-231" d.txt "BRGTEST-221" e.txt
I want to get... (22 Replies)
Discussion started by: pottic
22 Replies
9. Shell Programming and Scripting
So I'm stumped.
First... APOLOGIES... my work is offline in an office that has zero internet connectivity, as required by our client. If need be, I could print out my script attempts and retype them here. But on the off chance... here goes.
I have a text file (file_source) of terms, each line... (3 Replies)
Discussion started by: Brusimm
3 Replies
10. UNIX for Beginners Questions & Answers
I am trying to consolidate multiple information files (<hostname>.Linux.nfslist) into one file so that I can import it into Excel. I can get the file contents with cat *Linux.nfslist >> nfslist.txt. I need each line prefaced with the hostname. I am unsure how to do this.
--- Post updated at... (5 Replies)
Discussion started by: Kentlee65
5 Replies
LEARN ABOUT DEBIAN
bp_classify_hits_kingdom
BP_CLASSIFY_HITS_KINGDOM(1p) User Contributed Perl Documentation BP_CLASSIFY_HITS_KINGDOM(1p)
NAME
classify_hits_kingdom - classify BLAST hits by taxonomic kingdom
USAGE
classify_hits_kingdom [-i tab_file] [-i second_BLAST_file] [-e evalue_cutoff]
[-t dir_where_TAXONOMY_files_are] [-g gi2taxid]
[-z PATH_TO_zcat] [-v]
DESCRIPTION
Will print out the taxonomic distribution (at the kingdom level) for a set of hits against the NR database. This script assumes you've
done a search against the protein database, you'll have to make minor changes in the gi_taxid part to point to the gi_taxid_nuc.dump file.
This expects BLAST files in tabbed -m9 or -m8 format. Output with -m 8 or use blast2table.pl to convert (or fastam9_to_table.PLS if using
FASTA).
Input values:
-t/--taxonomy directory where the taxonomy .dmp files are (from NCBI)
-g/--gi Location of gi_taxid_prot.dmp (or gi_taxid_nucl.dmp if
the search was against a NT db)
-i/--in The name of the tab delimited -m8/-m9 output files to
process.
-e/--evalue Provide an E-value cutoff for hits to be considered
-z/--zcat Path to the 'zcat' executable, can also be 'gunzip -c'
if no zcat on your system.
Flags
-v/--verbose To turn on verbose messages
-h/--help Display this helpful information
This is intended to be useful starting script, but users may want to customize the output and parameters. Note that I am summarizing the
kingdoms here and Eukaryota not falling into Metazoa, Viridiplantae, or Fungi gets grouped into the general superkingdom Eukaryota. for
simplicity. There are comments in the code directing you to where changes can be made if you wanted to display hits by phylum for example.
Note that you must wipe out the cache file 'gi2class' that is craeed in your directory after making these changes.
AUTHOR
Jason Stajich jason_at_bioperl_dot_org
perl v5.14.2 2012-03-02 BP_CLASSIFY_HITS_KINGDOM(1p)