Sponsored Content
Top Forums UNIX for Dummies Questions & Answers How to change sequence name in along fasta file? Post 302663215 by baika on Wednesday 27th of June 2012 05:31:24 PM
Old 06-27-2012
How to change sequence name in along fasta file?

Hi
I have an alignment file (.fasta) with ~80 sequences. They look like this-

>JV101.contig00066(+):25302-42404|sequence_index=0|block_index=4|species=JV101|JV101_4_0

GAGGTTAATTATCGATAACGTTTAATTAAAGTGTTTAGGTGTCATAATTT
TAAATGACGATTTCTCATTACCATACACCTAAATTATCATCAATCTGAAT
TCAGATGTTTATTATAAAAATTAGATGAAAAATATGTTAATATACAAGTA

>JV501.contig00066(+):24356-42404|sequence_index=0|block_index=4|species=JV501|JV501_4_0
AATGACGATTTAGATGAAAAATAT...

The name of the sequences are too big and I want to just keep JV101, JV501 and delete rest of the words after the dot. I am new to unix, please suggest an easy unix command to do this.

Thanks

Baika
 

10 More Discussions You Might Find Interesting

1. Solaris

how to change the boot sequence?

Hi all, I have solaris and xp installed... Usually solaris occupies the first boot slot. So i wanted to know if it is possible to change the boot sequence to xp first and then solaris? (5 Replies)
Discussion started by: wrapster
5 Replies

2. Red Hat

change the shutdown sequence in red hat linux

Hi, I use red hat linux kernel 2.6 I want to add the application shutdown in shutdown sequence . I add the K script in /etc/rc.d/ all sub directory for all running level . But the auto shutdown application is not appear when I type "shutdown -r now" .. There is no indication the application... (5 Replies)
Discussion started by: chuikingman
5 Replies

3. Shell Programming and Scripting

How to change strt-up sequence of services on SLES10/11

Hi All, I want to change the start-up sequence of services on SLES10/11. I have my own start-up scripts for some services and I want them to start in a particular order(not in alphabetical order) Can anyone help me on this issue? (4 Replies)
Discussion started by: senrooy
4 Replies

4. Shell Programming and Scripting

How to change strt-up sequence of services on SLES10/11

Hi, Can anyone tell me what the following 2 lines are doing base=${0##*/} link=${base#*} I found this in a start up service script and I think it is giving the service link names which in turn will change the start up sequence of services. (3 Replies)
Discussion started by: senrooy
3 Replies

5. Shell Programming and Scripting

Parsing a fasta sequence with start and end coordinates

Hi.. I have a seperate chromosome sequences and i wanted to parse some regions of chromosome based on start site and end site.. how can i achieve this? For Example Chr 1 is in following format I need regions from 2 - 10 should give me AATTCCAAA and in a similar way 15- 25 should give... (8 Replies)
Discussion started by: empyrean
8 Replies

6. UNIX for Dummies Questions & Answers

Change sequence names in fasta file

I have fasta files with multiple sequences in each. I need to change the sequence name headers from: >accD:_59176-60699 ATGGAAAAGTGGAGGATTTATTCGTTTCAGAAGGAGTTCGAACGCA >atpA_(reverse_strand):_showing_revcomp_of_10525-12048 ATGGTAACCATTCAAGCCGACGAAATTAGTAATCTTATCCGGGAAC... (2 Replies)
Discussion started by: tyrianthinae
2 Replies

7. Shell Programming and Scripting

Extract sequence from fasta file

Hi, I want to match the sequence id (sub-string of line starting with '>' and extract the information upto next '>' line ). Please help . input > fefrwefrwef X900 AGAGGGAATTGG AGGGGCCTGGAG GGTTCTCTTC > fefrwefrwef X932 AGAGGGAATTGG AGGAGGTGGAG GGTTCTCTTC > fefrwefrwef X937... (2 Replies)
Discussion started by: ritakadm
2 Replies

8. 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

9. UNIX for Dummies Questions & Answers

Round up -FASTA file

I have the following script: awk 'FNR==NR{s+=$3;next;} { print $1 , $2, 100*$3/s }' and the following file: >P39PT-1224 Freq 900 cccctacgacggcattggtaatggctcagctgctccggatcccgcaagccatcttggatatgagggttcgtcggcctcttcagccaagg-cccccagcagaacatccagctgatcg >P39PT-784 Freq 2... (2 Replies)
Discussion started by: Xterra
2 Replies

10. UNIX for Beginners Questions & Answers

How to find a specific sequence pattern in a fasta file?

I have to mine the following sequence pattern from a large fasta file namely gene.fasta (contains multiple fasta sequences) along with the flanking sequences of 5 bases at starting position and ending position, AAGCZ-N16-AAGCZ Z represents A, C or G (Except T) N16 represents any of the four... (3 Replies)
Discussion started by: dineshkumarsrk
3 Replies
Bio::Tools::Run::Phylo::Semphy(3pm)			User Contributed Perl Documentation		       Bio::Tools::Run::Phylo::Semphy(3pm)

NAME
Bio::Tools::Run::Phylo::Semphy - Wrapper for Semphy SYNOPSIS
use Bio::Tools::Run::Phylo::Semphy; # Make a Semphy factory $factory = Bio::Tools::Run::Phylo::Semphy->new(); # Run Semphy with an alignment my $tree = $factory->run($alignfilename); # or with alignment object $tree = $factory->run($bio_simplalign); # you can supply an initial tree as well, which can be a newick tree file, # Bio::Tree::Tree object... $tree = $factory->run($bio_simplalign, $tree_obj); # ... or Bio::DB::Taxonomy object $tree = $factory->run($bio_simplalign, $bio_db_taxonomy); # (mixtures of all the above are possible) # $tree isa Bio::Tree::Tree DESCRIPTION
This is a wrapper for running the Semphy application by N. Friedman et a.. You can get details here: http://compbio.cs.huji.ac.il/semphy/. Semphy is used for phylogenetic reconstruction (making a tree with branch lengths from an aligned set of input sequences). You can try supplying normal Semphy command-line arguments to new(), eg. new(-hky => 1) or calling arg-named methods (excluding the initial hyphen(s), eg. $factory->hky(1) to set the --hky switch to true). Note that Semphy args are case-sensitive. To distinguish between Bioperl's -verbose and the Semphy's --verbose, you must set Semphy's verbosity with -semphy_verbose or the semphy_verbose() method. You will need to enable this Semphy wrapper to find the Semphy program. This can be done in (at least) three ways: 1. Make sure the Semphy executable is in your path. 2. Define an environmental variable SEMPHYDIR which is a directory which contains the Semphy application: In bash: export SEMPHYDIR=/home/username/semphy/ In csh/tcsh: setenv SEMPHYDIR /home/username/semphy 3. Include a definition of an environmental variable SEMPHYDIR in every script that will use this Semphy wrapper module, e.g.: BEGIN { $ENV{SEMPHYDIR} = '/home/username/semphy/' } use Bio::Tools::Run::Phylo::Semphy; FEEDBACK
Mailing Lists User feedback is an integral part of the evolution of this and other Bioperl modules. Send your comments and suggestions preferably to the Bioperl mailing list. Your participation is much appreciated. bioperl-l@bioperl.org - General discussion http://bioperl.org/wiki/Mailing_lists - About the mailing lists Support Please direct usage questions or support issues to the mailing list: bioperl-l@bioperl.org rather than to the module maintainer directly. Many experienced and reponsive experts will be able look at the problem and quickly address it. Please include a thorough description of the problem with code and data examples if at all possible. Reporting Bugs Report bugs to the Bioperl bug tracking system to help us keep track of the bugs and their resolution. Bug reports can be submitted via the web: http://redmine.open-bio.org/projects/bioperl/ AUTHOR - Sendu Bala Email bix@sendu.me.uk APPENDIX
The rest of the documentation details each of the object methods. Internal methods are usually preceded with a _ program_name Title : program_name Usage : $factory>program_name() Function: holds the program name Returns : string Args : None program_dir Title : program_dir Usage : $factory->program_dir(@params) Function: returns the program directory, obtained from ENV variable. Returns : string Args : None new Title : new Usage : $factory = Bio::Tools::Run::Phylo::Semphy->new() Function: creates a new Semphy factory Returns : Bio::Tools::Run::Phylo::Semphy Args : Most options understood by Semphy can be supplied as key => value pairs, with a true value for switches. These options can NOT be used with this wrapper (they are handled internally or don't make sense in this context): -h | --help | --fill-help -s | --sequence -t | --tree To distinguish between Bioperl's -verbose and the Semphy's --verbose, you must set Semphy's verbosity with -semphy_verbose run Title : run Usage : $result = $factory->run($fasta_align_file); -or- $result = $factory->run($align_object); -or- $result = $factory->run($fasta_align_file, $newick_tree_file); -or- $result = $factory->run($align_object, $tree_object); -or- $result = $factory->run($align_object, $db_taxonomy_object); Function: Runs Semphy on an alignment. Returns : Bio::Tree::Tree Args : The first argument represents an alignment, the second (optional) argument a species tree (to set an initial tree: normally the -t option to Semphy). The alignment can be provided as a multi-fasta format alignment filename, or a Bio::Align::AlignI compliant object (eg. a Bio::SimpleAlign). The species tree can be provided as a newick format tree filename or a Bio::Tree::TreeI compliant object. Alternatively a Bio::DB::Taxonomy object can be supplied, in which case the species tree will be generated by using the alignment sequence names as species names and looking for those in the supplied database. In all cases where an initial tree was supplied, the alignment sequence names must correspond to node ids in the species tree. _setparams Title : _setparams Usage : Internal function, not to be called directly Function: Creates a string of params to be used in the command string Returns : string of params Args : alignment and tree file names perl v5.12.3 2011-06-18 Bio::Tools::Run::Phylo::Semphy(3pm)
All times are GMT -4. The time now is 05:07 PM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy