Dear Rudic,
Below script can rename the characters in gene1.fasta specified in list.txt.
However, it is not based on the position. Its purely based on the matching strings between the two files. But, here my problem is different, I tried workout like this
, but my list.txt is not having common strings, so I can not rename sequentially. That is why I seek your help.
Last edited by dineshkumarsrk; 11-13-2019 at 04:12 AM..
Hi ,
I have a typical situation. I have 4 files and with different headers (number of headers is varible ).
I need to make such a merged file which will have headers combined from all files (comman coluns should appear once only).
For example -
File 1
H1|H2|H3|H4
11|12|13|14
21|22|23|23... (1 Reply)
Hi Guys,
While I was writing one shell script , I just got struck at this point.
I need to extract words from a file at some specified position and do some comparison operation and need to replace the extracted word with another word.
Eg : I like Orange very much.
I need to replace... (19 Replies)
Hi,
I am new to unix. I want to delete 2 words placed at position say for example at 23rd and 45th position in a line. I used sed but couldnt achieve this.
Example: the file contains 2 lines
12345 98765 "12345" 876
12345 98765 "64578" 876
I want to delete " placed at position 13 and 19... (4 Replies)
I have a file with thousands of sequences that looks like this:
I need to replace the headers using a second file
Thus, I will end up having the following file:
I am looking for an AWK script that I can easily plug in my current pipeline.
Any help will be greatly appreciated! (6 Replies)
Hi, I have a file1 of many long sequences, each preceded by a unique header line. file2 is 3-columns list: headers name, start position, end position. I'd like to extract the sequence region of file1 specified in file2.
Based on a post elsewhere, I found the code:
awk... (2 Replies)
Hi,
I am unable to find the right option to extract the data in the fixed width file.
sample data
abcd1234xgyhsyshijfkfk
hujk9876 io xgla
loki8787eljuwoejroiweo
dkfj9098 dja
Search based on position 8-9="xg" and print the entire row
output
... (4 Replies)
OS : Linux 2.6x
Shell : Korn
In a single file , how can I identify all the Uniqe values at a specific character position and length of each record ,
and simultaneously SPLIT the records of the file based on each of these values and write them in seperate files .
Lets say :
a) I want to... (4 Replies)
I have two files. File1 is shown below.
>153L:B|PDBID|CHAIN|SEQUENCE
RTDCYGNVNRIDTTGASCKTAKPEGLSYCGVSASKKIAERDLQAMDRYKTIIKKVGEKLCVEPAVIAGIISRESHAGKVL
KNGWGDRGNGFGLMQVDKRSHKPQGTWNGEVHITQGTTILINFIKTIQKKFPSWTKDQQLKGGISAYNAGAGNVRSYARM
DIGTTHDDYANDVVARAQYYKQHGY
>16VP:A|PDBID|CHAIN|SEQUENCE... (7 Replies)
Hi,
I have a file with multiple lines(fixed width dat file). I want to search for '02' in the positions 45-46 and if available, in that lines, I need to replace value in position 359 with blank. As I am new to unix, I am not able to figure out how to do this. Can you please help me to achieve... (9 Replies)
Discussion started by: Pradhikshan
9 Replies
LEARN ABOUT DEBIAN
bp_search2tribe
BP_SEARCH2TRIBE(1p) User Contributed Perl Documentation BP_SEARCH2TRIBE(1p)NAME
search2tribe - Turn SearchIO parseable reports(s) into TRIBE matrix
SYNOPSIS
Usage:
search2tribe [-o outputfile] [-f reportformat] [-w/--weight] file1 file2 ..
DESCRIPTION
This script is probably too slow for most people's uses. It is better to use something like scripts/searchio/fastam9_to_table, -m 9 output
from BLAST, or the blast2table from the BLAST O'Reilly book to get a tabular output from these programs and then feed the table into MCL
with the mcxdeblast script and the --m9 option.
This script will turn a protein Search report (BLASTP, FASTP, SSEARCH) into a Markov Matrix for TribeMCL clustering.
The options are:
-o filename - the output filename [default STDOUT]
-f format - search result format (blast, fasta)
(ssearch is fasta format). default is blast.
-w or --weight VALUE - Change the default weight for E(0.0) hits
to VALUE (default=200 (i.e. 1e-200) )
-h - this help menu
Additionally specify the filenames you want to process on the command-line. If no files are specified then STDIN input is assumed. You
specify this by doing: search2tribe < file1 file2 file3
AUTHOR
Jason Stajich, jason-at-bioperl-dot-org
perl v5.14.2 2012-03-02 BP_SEARCH2TRIBE(1p)