Sponsored Content
Top Forums UNIX for Dummies Questions & Answers Search for a particular word and replace the first character Post 302904592 by Fahmida on Thursday 5th of June 2014 04:12:34 AM
Old 06-05-2014
Search for a particular word and replace the first character

Hi Unix gurus,

I've a dna sequence in a file format known as fasta format (sequence header starts with > and ignored), an example shown below:
Code:
>sequence_1
CGTATTCTCCGAATACC
ATACG
>sequence_2
CAGATTTTCAAATACCCCC

In a file like this I want to do the following three search and replace. The original file becomes three different files based on each search condition.
1) search two character pattern CG and replace that C by the word NULL
2) search three character pattern CXG and replace that C by the word NULL (where X=A,C or T but not G)
3) search three character pattern CXX and replace that C by the word NULL (where X=A,C or T but not G)

Example output based on above file:
File_CG_replace
Code:
>sequence_1
NULLGTATTCTCNULLGAATACC
ATANULLG
>sequence_2
CAGATTTTCAAATACCCCC

File_CXG_replace
Code:
>sequence_1
CGTATTCTNULLCGAATACC
ATACG
>sequence_2
CAGATTTTCAAATACCCCC

File_CXX_replace
Code:
>sequence_1
CGTATTNULLTCNULLGAATANULLC
ATACG
>sequence_2
CAGATTTTNULLAAATANULLCCCC

Thanks for your help.
 

10 More Discussions You Might Find Interesting

1. UNIX for Dummies Questions & Answers

to search for a particular character in a word

Hi I would like to accept in a string from user like username/pwd@dbname suppose the user does not input @ then i should throw an error that @ symbol missing . How to achieve this Thanks in advance Suresh (6 Replies)
Discussion started by: ssuresh1999
6 Replies

2. Shell Programming and Scripting

sed search and replace word assistance...

Hi, I am trying to write a shell script designed to take input line by line by line from a file with a word on each line for editing with sed. Example file: 1.ejverything 2.bllown 3.maikling 4.manegement 5.existjing 6.systems My design currently takes input from the user, and... (2 Replies)
Discussion started by: mkfitzwilliams
2 Replies

3. Shell Programming and Scripting

Need to search and replace based on character count

Hi, I wanted to add a newline character after every 100 characters in a file using a awk or shell without reading each line of the file. I want to run a command on the complete file. This does based on a string but i want to add a new line after every 100 characters ir-respective of the... (3 Replies)
Discussion started by: vijaykrc
3 Replies

4. UNIX for Dummies Questions & Answers

vi Search for text, Replace with <CR> or control character.

Greetings, Using vi, how can I change the following text: -I/myviews/nexus_7400rel/vobs/nexus/platforms/97400/include -I/myviews/nexus_7400rel/vobs/nexus/modules/i2c/7400/include -I/myviews/nexus_7400rel/vobs/nexus/modules/surface/7400/include Into this:... (4 Replies)
Discussion started by: omega949
4 Replies

5. Shell Programming and Scripting

Interesting question - Search and replace the word after sign "="

Hi Guys, Req your help in searching and replacing the word that comes after equals(=) symbol I would like to replace the sting in bold with a string in variable. d=ABCDF8C44C22 # grep -i NIM_MASTERID ${_NIMINFO} export NIM_MASTERID=00CDF8C44C00 I'm looking to replace any word that... (4 Replies)
Discussion started by: ajilesh
4 Replies

6. Shell Programming and Scripting

Search for word in a xml file and replace it with something else

Hello Unix Users, I am very new to Unix so I am not sure how do I do the following. I need a script such that when I type the following in the command prompt > . scriptName.sh wordToBeReplaced DirectoryLocation will find the word someword located in a somefile.xml in DirectoryLocation... (8 Replies)
Discussion started by: 5211171
8 Replies

7. Shell Programming and Scripting

Regex:search/replace but not for escaped character

Hi Input: - -- --- ---- aa-bb-cc aa--bb--cc aa---bb---cc aa----bb----cc Output: . - -. -- aa.bb.cc (7 Replies)
Discussion started by: chitech
7 Replies

8. Shell Programming and Scripting

Search a character and replace it with multiple lines

This is for AIX 6.1, I've a flat file and the format is like this DECLARE some statements; BEGIN some statements; END; I've to search BEGIN and replace it with the following 4 lines BEGIN For x in 1..1 LOOP BEGIN Similarly I've to search END and replace it with the... (7 Replies)
Discussion started by: Mukul Sharma
7 Replies

9. Shell Programming and Scripting

Search and replace the string with new word using xml tags

Hi All i need to replace the url1 inside <remote> tag in below xml in first instance and in the second instance with url2. any help appreciated <locations> <hudson.scm.SubversionSCM_-ModuleLocation> <remote>https://svn2015.com/svn/repos/internalshard</remote> ... (4 Replies)
Discussion started by: madankumar.t@hp
4 Replies

10. Post Here to Contact Site Administrators and Moderators

Search for a pattern and replace a space at specific position with a Character in File

In file, we have millions of records each of 1000 in length. And at specific position say 800 there is a space, we need to replace it with Character X if the ID in that row starts with 123. So far i have used the below which is replacing space at that position to X but its not checking for... (3 Replies)
Discussion started by: Jagmeet Singh
3 Replies
ABACAS(1)							   User Commands							 ABACAS(1)

NAME
abacas - Algorithm Based Automatic Contiguation of Assembled Sequences SYNOPSIS
abacas -r ref -q qs -p prog [OPTIONS] OR abacas -r ref -q psf -e ref reference sequence in a single fasta file qs contigs in multi-fasta format rog MUMmer program to use: 'nucmer' or 'promer' psf pseudomolecule/ordered sequence file in fasta format OPTIONS -h print usage -d use default nucmer/promer parameters -s int minimum length of exact matching word (nucmer default = 12, promer default = 4) -m print ordered contigs to file in multifasta format -b print contigs in bin to file -N print a pseudomolecule without "N"s -i int mimimum percent identity [default 40] -v int mimimum contig coverage [default 40] -V int minimum contig coverage difference [default 1] -l int minimum contig length [default 1] -t run tblastx on contigs that are not mapped -g string (file name) print uncovered regions (gaps) on reference to file name -a append contigs in bin to the pseudomolecule -o prefix output files will have this prefix -P pick primer sets to close gaps -f int number of flanking bases on either side of a gap for primer design (default 350) -R int Run mummer [default 1, use -R 0 to avoid running mummer] -e Escape contig ordering i.e. go to primer design -c Reference sequence is circular DESCRIPTION
ABACAS is intended to rapidly contiguate (align, order, orientate), visualize and design primers to close gaps on shotgun assembled contigs based on a reference sequence. ABACAS uses MUMmer to find alignment positions and identify syntenies of assembled contigs against the reference. The output is then pro- cessed to generate a pseudomolecule taking overlapping contigs and gaps in to account. ABACAS generates a comparision file that can be used to visualize ordered and oriented contigs in ACT. Synteny is represented by red bars where colour intensity decreases with lower values of percent identity between comparable blocks. Information on contigs such as the orientation, percent identity, coverage and overlap with other contigs can also be visualized by loading the outputted feature file on ACT. AUTHOR
ABACAS IS Copyright (C) 2008-10 The Wellcome Trust Sanger Institute, Cambridge, UK. This manual page was written by Andreas Tille <tille@debian.org>, for the Debian project (and may be used by others). 1.3.1 2011-02-11 ABACAS(1)
All times are GMT -4. The time now is 02:15 PM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy