Still doesn't produce what I need. I have 725 records, in each record contein the match ORGANISM with 4 lines on each. I would like to have at the end 725 lines with this four lines in one
example:
if I use
I got :
and what I need is this:
---------- Post updated at 01:28 PM ---------- Previous update was at 01:13 PM ----------
Thanks Alister! It worked
Last edited by Scott; 09-18-2010 at 05:43 PM..
Reason: Code tags, please...
Hi folks,
Lets say I have the following text file:
name, lastname, 1234, name.lastname@test.com
name1, lastname1, name2.lastname2@test.com, 2345
name, 3456, lastname, name3.lastname3@test.com
4567, name, lastname, name4.lastname4@test.com
I now need the following output:
1234... (5 Replies)
greetings citizens of Unix.com
I am perplexed with an issue.
The issue is trying to print the last 5 characters of a string in PERL.
Below are demonstrated my daft attempts at performing the forementioned task.
$row =~ m/^.*(.....)\s$/;
$row =~ m/\w{5}\s*$/i;$row =~... (3 Replies)
Hi all,
This should be very easy but I can't figure it out...
I have a file that looks like this:
@SRR057408.1 FW8Y5CK02R652T length=34
AGCAGTGGTATCAACGCAGAGTAAGCAGTGGTAT
+SRR057408.1 FW8Y5CK02R652T length=34
FIIHFF6666?=:88@@@BBD:::?@ABBAAA>8
@SRR057408.2 FW8Y5CK02TBMHV length=52... (1 Reply)
Hi again,
I'm looking for some help with nawk, I can print a line which has a regex match in it from a file using /pattern/ but I'm looking for a way to only print the $tring which contains the pattern, rather than the whole line.
This $tring may be of variable length, may occur at any point... (1 Reply)
I have a file of protein sequences with headers (my source file). Based on a list of IDs (which are included in some of the headers), I'd like to print out only the specified sequences, with only the ID as header.
In other words, I'd like to search source.txt for the terms in IDs.txt, and print... (3 Replies)
Hello,
I need an awk command to print only the lines that match regex on xth field from file.
For example if I use this command
awk -F"|" ' $22 == "20130117090000.*" 'It wont work, I think, because single quotes wont allow the usage of the metacharacter star * . On the other hand I dont know... (2 Replies)
Hi,
I am using the following code to fetch lines that are generated in last 1 hr . Hence, I am using date function to calculate -last 1 hr & the current hr and then somehow use awk (or sed-if someone could guide me better)
with some regex pattern.
dt_1=`date +%h" "%d", "%Y\ %l -d "1 hour... (10 Replies)
Hello.
I have been looking high and low for the solution for this. I seems there should be a simple answer, but alas.
I have a big xml file, and I need to extract certain information from specific items. The information I need can be found between a specific set of tags. let's call them... (2 Replies)
I have a directory of files, I can show the number of lines in each file and order them from lowest to highest with:
wc -l *|sort
15263 Image.txt
16401 reference.txt
40459 richtexteditor.txt
How can I also print the number of unique lines in each file?
15263 1401 Image.txt
16401... (15 Replies)
Discussion started by: spacegoose
15 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)