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:
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
File_CXG_replace
File_CXX_replace
Thanks for your help.
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)
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)
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)
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)
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)
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)
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)
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
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
LEARN ABOUT DEBIAN
2ndscore
2NDSCORE(1) User Contributed Documentation 2NDSCORE(1)NAME
2ndscore - find the best hairpin anchored at each position.
SYNOPSIS
2ndscore in.fasta > out.hairpins
DESCRIPTION
For every position in the sequence this will output a line:
-0.6 52 .. 62 TTCCTAAAGGTTCCA GCG CAAAA TGC CATAAGCACCACATT
(score) (start .. end) (left context) (hairpin) (right contenxt)
For positions near the ends of the sequences, the context may be padded with 'x' characters. If no hairpin can be found, the score will be
'None'.
Multiple fasta files can be given and multiple sequences can be in each fasta file. The output for each sequence will be separated by a
line starting with '>' and containing the FASTA description of the sequence.
Because the hairpin scores of the plus-strand and minus-strand may differ (due to GU binding in RNA), by default 2ndscore outputs two sets
of hairpins for every sequence: the FORWARD hairpins and the REVERSE hairpins. All the forward hairpins are output first, and are
identified by having the word 'FORWARD' at the end of the '>' line preceding them. Similarly, the REVERSE hairpins are listed after a '>'
line ending with 'REVERSE'. If you want to search only one or the other strand, you can use:
--no-fwd Don't print the FORWARD hairpins
--no-rvs Don't print the REVERSE hairpins
You can set the energy function used, just as with transterm with the --gc, --au, --gu, --mm, --gap options. The --min-loop, --max-loop,
and --max-len options are also supported.
FORMAT OF THE .BAG FILES
The columns for the .bag files are, in order:
1. gene_name
2. terminator_start
3. terminator_end
4. hairpin_score
5. tail_score
6. terminator_sequence
7. terminator_confidence: a combination of the hairpin and tail score that
takes into account how likely such scores are in a random sequence. This
is the main "score" for the terminator and is computed as described in
the paper.
8. APPROXIMATE_distance_from_end_of_gene: The *approximate* number of base
pairs between the end of the gene and the start of the terminator. This
is approximate in several ways: First, (and most important) TransTermHP
doesn't always use the real gene ends. Depending on the options you give
it may trim some off the ends of genes to handle terminators that
partially overlap with genes. Second, where the terminator "begins"
isn't that well defined. This field is intended only for a sanity check
(terminators reported to be the best near the ends of genes shouldn't be
_too far_ from the end of the gene).
USING TRANSTERM WITHOUT GENOME ANNOTATIONS
TransTermHP uses known gene information for only 3 things: (1) tagging the putative terminators as either "inside genes" or "intergenic,"
(2) choosing the background GC-content percentage to compute the scores, because genes often have different GC content than the intergenic
regions, and (3) producing slightly more readable output. Items (1) and (3) are not really necessary, and (2) has no effect if your genes
have about the same GC-content as your intergenic regions.
Unfortunately, TransTermHP doesn't yet have a simple option to run without an annotation file (either .ptt or .coords), and requires at
least 2 genes to be present. The solution is to create fake, small genes that flank each chromosome. To do this, make a fake.coords file
that contains only these two lines:
fakegene1 1 2 chome_id
fakegene2 L-1 L chrom_id
where L is the length of the input sequence and L-1 is 1 less than the length of the input sequence. "chrom_id" should be the word directly
following the ">" in the .fasta file containing your sequence. (If, for example, your .fasta file began with ">seq1", then chrom_id =
seq1).
This creates a "fake" annotation with two 1-base-long genes flanking the sequence in a tail-to-tail arrangement: --> <--. TransTermHP can
then be run with:
transterm -p expterm.dat sequence.fasta fake.coords
If the G/C content of your intergenic regions is about the same as your genes, then this won't have too much of an effect on the scores
terminators receive. On the other hand, this use of TransTermHP hasn't been tested much at all, so it's hard to vouch for its accuracy.
SEE ALSO transterm(1)AUTHOR
Alex Mestiashvili <alex@biotec.tu-dresden.de>
2011-02-19 2NDSCORE(1)