Sponsored Content
Top Forums Shell Programming and Scripting replacement and "reading frames" Post 302433377 by jossojjos on Tuesday 29th of June 2010 11:02:03 AM
Old 06-29-2010
Thanks folks !

Actually, counting is just one of the things I'd like to do frequently ; replacing is another.

I can see quite easily how I could do this using the "fold" command (thanks for bringing it to my attention!)

Code:
echo "ACGTACGTACGT" | fold -w3 | sed -e 's/TAC/NNN/g' | tr -d '\n'

Maybe the perl solutions are quicker ? I'll have to test that, but I'll have to know perl a bit better.

Thanks again !
jos
 

9 More Discussions You Might Find Interesting

1. UNIX for Dummies Questions & Answers

Explain the line "mn_code=`env|grep "..mn"|awk -F"=" '{print $2}'`"

Hi Friends, Can any of you explain me about the below line of code? mn_code=`env|grep "..mn"|awk -F"=" '{print $2}'` Im not able to understand, what exactly it is doing :confused: Any help would be useful for me. Lokesha (4 Replies)
Discussion started by: Lokesha
4 Replies

2. UNIX for Dummies Questions & Answers

Vi - "The replacement pattern is too long"

Hi, I am trying to replace a value in a script with another value. I am performing a vi command from another script. vi - ${conf_path}/CANCEL_CD_PART2.txt<<! :%s/RANGE/${btch_range}/g :wq ! 'RANGE' is the current value that the parm in the other script has (PARM1=RANGE), along with... (3 Replies)
Discussion started by: hern14
3 Replies

3. Shell Programming and Scripting

cat $como_file | awk /^~/'{print $1","$2","$3","$4}' | sed -e 's/~//g'

hi All, cat file_name | awk /^~/'{print $1","$2","$3","$4}' | sed -e 's/~//g' Can this be done by using sed or awk alone (4 Replies)
Discussion started by: harshakusam
4 Replies

4. Shell Programming and Scripting

awk command to replace ";" with "|" and ""|" at diferent places in line of file

Hi, I have line in input file as below: 3G_CENTRAL;INDONESIA_(M)_TELKOMSEL;SPECIAL_WORLD_GRP_7_FA_2_TELKOMSEL My expected output for line in the file must be : "1-Radon1-cMOC_deg"|"LDIndex"|"3G_CENTRAL|INDONESIA_(M)_TELKOMSEL"|LAST|"SPECIAL_WORLD_GRP_7_FA_2_TELKOMSEL" Can someone... (7 Replies)
Discussion started by: shis100
7 Replies

5. AIX

"Frames" and "Words" in fcstat output

What are "Frames" and "Words" in the fcstat output? vio1:/home/padmin:# fcstat fcs0 <snip> Transmit Statistics Receive Statistics ------------------- ------------------ Frames: 122844229 363445456 Words: 50940091904 171210861568 <snip> The... (1 Reply)
Discussion started by: kah00na
1 Replies

6. Shell Programming and Scripting

"Encrypting" using number replacement

Hi, I have a task of "encrypting" a file by replacing the numbers with another set of numbers. Sort of swapping the values, say: 1 = 5 2 = 8 3 = 7 4 = 1 5 = 9 and so on.. so if i have 12345, my output should be 58719. problem is i get 98719 since after swapping 1 to 5, my sed... (7 Replies)
Discussion started by: agentgrecko
7 Replies

7. UNIX for Dummies Questions & Answers

Using "mailx" command to read "to" and "cc" email addreses from input file

How to use "mailx" command to do e-mail reading the input file containing email address, where column 1 has name and column 2 containing “To” e-mail address and column 3 contains “cc” e-mail address to include with same email. Sample input file, email.txt Below is an sample code where... (2 Replies)
Discussion started by: asjaiswal
2 Replies

8. Shell Programming and Scripting

Bash script - Print an ascii file using specific font "Latin Modern Mono 12" "regular" "9"

Hello. System : opensuse leap 42.3 I have a bash script that build a text file. I would like the last command doing : print_cmd -o page-left=43 -o page-right=22 -o page-top=28 -o page-bottom=43 -o font=LatinModernMono12:regular:9 some_file.txt where : print_cmd ::= some printing... (1 Reply)
Discussion started by: jcdole
1 Replies

9. AIX

Apache 2.4 directory cannot display "Last modified" "Size" "Description"

Hi 2 all, i have had AIX 7.2 :/# /usr/IBMAHS/bin/apachectl -v Server version: Apache/2.4.12 (Unix) Server built: May 25 2015 04:58:27 :/#:/# /usr/IBMAHS/bin/apachectl -M Loaded Modules: core_module (static) so_module (static) http_module (static) mpm_worker_module (static) ... (3 Replies)
Discussion started by: penchev
3 Replies
BLASTCLUST(1)						     NCBI Tools User's Manual						     BLASTCLUST(1)

NAME
blastclust - BLAST score-based single-linkage clustering SYNOPSIS
blastclust [-] [-C] [-L X] [-S X] [-W N] [-a N] [-b F] [-c filename] [-d filename] [-e F] [-i filename] [-l filename] [-o filename] [-p F] [-r filename] [-s filename] [-v [filename]] DESCRIPTION
blastclust automatically and systematically clusters protein or DNA sequences based on pairwise matches found using the BLAST algorithm in case of proteins or Mega BLAST algorithm for DNA. In the latter case a single Mega BLAST search is performed for all the sequences combined against a database created from the same sequences. blastclust finds pairs of sequences that have statistically significant matches and clusters them using single-linkage clustering. OPTIONS
A summary of options is included below. - Print usage message -C Complete unfinished clustering -L X Length coverage threshold (default = 0.9) -S X Score coverage threshold (bit score / length if < 3.0, percentage of identities otherwise; default = 1.75) -W N Use words of size N (length of best perfect match; zero invokes default behavior: 3 for proteins, 32 for nucleotides) -a N Number of CPU's to use (default = 1) -b F Do not require coverage on both neighbours -c filename Read advanced options from configuration file filename -d filename Input as a database -e F Disable id parsing in database formatting -i filename FASTA input file (program will format the database and remove files in the end; default = stdin) -l filename Restrict reclustering to id list in filename -o filename Output file for list of clusters (default = stdout) -p F Input is nucleotides, not proteins. -r filename Restore neighbors for reclustering from filename -s filename Save all neighbours to filename -v [filename] Print verbose progress messages (to filename) AUTHOR
The National Center for Biotechnology Information. SEE ALSO
blast(1), formatdb(1), /usr/share/doc/blast2/blastclust.html, <http://www.ncbi.nlm.nih.gov/BLAST/> NCBI
2004-06-25 BLASTCLUST(1)
All times are GMT -4. The time now is 02:42 PM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy