Login or Register to Ask a Question and Join Our Community

Sponsored Content
Full Discussion: Using = with sed to increase sequence count
Top Forums Shell Programming and Scripting Using = with sed to increase sequence count Post 303024203 by Xterra on Tuesday 2nd of October 2018 11:02:47 PM
Old 10-03-2018
Using = with sed to increase sequence count
I have a fasta file like this one:
Code:
>ID1
AAAAAA
>ID2
TTTTTT

And I am using this sed script to increase the count sequence
Code:
sed '/^>/s/.*//;/^$/=;/^$/d' text.txt | sed 's/[1-9].*/echo ">seq" \$(( ( & + 1 )\/2 ))/e'

I get the desired output:
Code:
>seq 1
AAAAAA
>seq 2
TTTTTT

However, this doesn't work and I do not understand why:
Code:
sed '/^>/s/.*//;/^$/=;/^$/d;s/[1-9].*/echo ">seq" \$(( ( & + 1 )\/2 ))/e' text.txt

I was hoping someone here would help me understand the issue. Moreover, I was hoping I could get a better, more elegant sed solution. While perl or awk might be more appropriate, I am actually looking for 100% sed approach.
Thanks in advance
 

10 More Discussions You Might Find Interesting

1. UNIX for Dummies Questions & Answers

IPv4 addresses: count/output and Awk/Sed

Hi forum. I am fairly new to scripting and use a simple script to process e-mails for my work. These e-mails contain a list of IPv4 IPs that I process and seperate into text files, which are then attached to a larger, 'digest' e-mail. I also put some of the output from the text files into the... (4 Replies)
Discussion started by: laebshade
4 Replies

2. UNIX for Dummies Questions & Answers

Record count problem using sed command

Hi, I have a script which removes 2 header records and 1 trailer record in a list of files. The commands doing the actions are sed '1,2d' $file > tempfile1.dat sed '$d' < tempfile1.dat > $output.txt Its working fine for all records except a file having size=1445509814 and number of... (2 Replies)
Discussion started by: ayanbiswas
2 Replies

3. UNIX for Dummies Questions & Answers

count number of fields not using SED or AWK

hi forums i need help with a little problem i am having. i need to count the number of fields that are in a saved variable so i can use that number to make a different function work properly. is there a way of doing this without using SED/AWK? anything would be greatly appreciated (4 Replies)
Discussion started by: strasner
4 Replies

4. Shell Programming and Scripting

Increase sed performance

I'm using sed to do find and replace. But since the file is huge and i have more than 1000 files to be searched, the script is taking a lot of time. Can somebody help me with a better sed command. Below is the details. Input: 1 1 2 3 3 4 5 5 Here I know the file is sorted. ... (4 Replies)
Discussion started by: gpaulose
4 Replies

5. UNIX for Dummies Questions & Answers

Grep char count & pipe to sed command

Hi I am having a 'grep' headache Here is the contents of my file: (PBZ,CP,(((ME,PBZ,BtM),ON),((ME,((PBZ,DG),(CW9,PG11))),CW9,TS2,RT1))) I would like to count out how many times 'PBZ' occurs and then place that number in the line above 3... (8 Replies)
Discussion started by: cavanac2
8 Replies

6. Shell Programming and Scripting

count and number instances of a character in sed or awk

I currently use LaTeX together with a sed script to set cloze test papers for my students. I currently pepend and equals sign to the front of the words I want to leave out in the finished test, =perpendicular, for example. I am able to number the blanks using a variable in LaTeX. I would like to... (8 Replies)
Discussion started by: maouinin
8 Replies

7. Shell Programming and Scripting

How to get count of replacements done by sed?

Hi , How can i get count of replacements done by sed in a file. I know grep -c is a method. But say if sed had made 10 replacement in a file, can i get number 10 some how? (8 Replies)
Discussion started by: abhitanshu
8 Replies

8. Shell Programming and Scripting

find common entries and match the number with long sequence and cut that sequence in output

Hi all, I have a file like this ID 3BP5L_HUMAN Reviewed; 393 AA. AC Q7L8J4; Q96FI5; Q9BQH8; Q9C0E3; DT 05-FEB-2008, integrated into UniProtKB/Swiss-Prot. DT 05-JUL-2004, sequence version 1. DT 05-SEP-2012, entry version 71. FT COILED 59 140 ... (1 Reply)
Discussion started by: manigrover
1 Replies

9. Shell Programming and Scripting

Ignore escape sequence in sed

Friends, In the file i am having more then 100 lines like, File1 had the values like this: #Example East.server_01=EAST.SERVER_01 East.server_01=EAST.SERVER_01 West.server_01=WEST.SERVER_01 File2 had the values like this: #Example EAST.SERVER_01=http://yahoo.com... (3 Replies)
Discussion started by: jothi basu
3 Replies

10. Shell Programming and Scripting

Count and search by sequence in multiple fasta file

Hello, I have 10 fasta files with sequenced reads information with read sizes from 15 - 35 . I have combined the reads and collapsed in to unique reads and filtered for sizes 18 - 26 bp long unique reads. Now i wanted to count each unique read appearance in all the fasta files and make a table... (5 Replies)
Discussion started by: empyrean
5 Replies

LEARN ABOUT DEBIAN

fixnt

fixnt(1)							      Debian								  fixnt(1)

NAME

       fixnt - Filter for the Windows NT postscript printer driver.

SYNOPSIS

       fixnt < BADFILE.ps > GOODFILE.ps

DESCRIPTION

       The  Windows  NT  postscript  driver has a tendency to make broken postscript files, that are incompatible with psutils.  fixnt is a filter
       that fixes these problems, allowing the use of psnup(1).

       The filter takes the broken postscript file on stdin, and outputs a fixed postscript file on stdout.  It has no other form  for	invocation
       and takes no options on the command-line.

OPTIONS

       fixnt  takes no options.

BUGS

       fixnt does not check for NTPSOct94.  For a workaround, use a sed(1) command to replace 'NTPSOct94' with 'NTPSOct95', like so:
	      sed 's/NTPSOct94/NTPSOct95/g'
       This is particularly important for Windows NT 3.5 users.

AUTHOR

       fixnt   was   written   by   Holger   Bauer   <Holger.Bauer@topmail.de>,  Michael  Rath	<rath@itsm.uni-stuttgart.de>,  and  Akim  Demaille
       <demaille@inf.enst.fr>.

REPORTING BUGS

       Report bugs to the Authors, but avoid sending large postscript files.

       Patches are always welcome; send to <bauer@itsm.uni-stuttgart.de>.

SEE ALSO

       psnup(1), sed(1)

a2ps								   February 2003							  fixnt(1)
All times are GMT -4. The time now is 05:28 PM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy