Login or Register to Ask a Question and Join Our Community

Sponsored Content
Full Discussion: Changing from FASTA to PHYLIP format
Top Forums Shell Programming and Scripting Changing from FASTA to PHYLIP format Post 302499225 by Xterra on Wednesday 23rd of February 2011 05:35:10 PM
Old 02-23-2011
I need to combine both codes
Code:
awk '$1=substr($1"       ",1,8)"  "' FS="\n" OFS= RS=\> file

Code:
awk '{s=length($2)}END{print NR-1, s}' file

So I can get the desired output
Code:
3 100 
SeID1_____AACCATGACAGAGGAGATGTGAACAGATAGAGGGATGACAGATGACAGATAGACCCAGACTGACAGGTTCAAAGGCTGCAGTGCAGTGACGTGACGATTT 
Sequence__AACCATGACAGAGGAGATGTGAACAGATAGAGGGATGACAGATGACAGATAGACCCAGACTGACAGGTTCAAAGGCTGCAGTGCAGTGACGTGACGATTT 
Seq-39____AACCATGACAGAGGAGATGTGAACAGATAGAGGGATGACAGATGACAGATAGACCCAGACTGACAGGTTCAAAGGCTGCAGTGCAGTGACGTGACGATTT

I have been trying but I just cannot get the code to do what I want.
Can anyone explain me how can I combine them?
Thanks!
Last edited by Xterra; 02-23-2011 at 06:59 PM..
 

10 More Discussions You Might Find Interesting

1. Shell Programming and Scripting

Changing date format

Hi, Is there any way to change one date format to another ?? I mean I have a file having dates in the format (Thu Sep 29 2005) ... and i wud like to change these to YYYYMMDD format .. is there any command which does so ?? Or anything like enum which we have in C ?? Thanks in advance, ... (5 Replies)
Discussion started by: Sabari Nath S
5 Replies

2. UNIX for Dummies Questions & Answers

Changing the format of date

Hi, There are lots of threads about how to manipulate the date using date +%m %....... But how can I change the default format of the commad date? $ date Mon Apr 10 10:57:15 BST 2006 This would be on fedora and SunOs. Cheers, Neil (4 Replies)
Discussion started by: nhatch
4 Replies

3. Shell Programming and Scripting

changing format

Dear Experts, Currently my script is gereating the output like this as mentioned below. 8718,8718,0,8777 7450,7450,0,7483 5063,5063,0,5091 3840,3840,0,3855 3129,3129,0,3142 2400,2400,0,2419 2597,2597,0,2604 3055,3055,0,3078 4249,4249,0,4266 4927,4927,0,4957 8920,8920,0,8978... (4 Replies)
Discussion started by: shary
4 Replies

4. Shell Programming and Scripting

Changing date format

Hi, I have a column in a table of Timestamp datatype. For Example : Var1 is the column 2008-06-26-10.10.30.2006. I have Given query as date(var1) and time (var1) I got the file as in the below format : File1: Col1 Col2 2008-06-02|12.36.06 2008-06-01|23.36.35 But the problem is... (7 Replies)
Discussion started by: manneni prakash
7 Replies

5. Shell Programming and Scripting

changing month in Mmm format to mm FORMAT

i have an variable mydate=2008Nov07 i want o/p as in variable mymonth=11 (i.e nov comes on 11 number month) i want some command to do this for any month without using any loop. plz help me (1 Reply)
Discussion started by: RahulJoshi
1 Replies

6. UNIX for Dummies Questions & Answers

fasta format?

Hi, I'm in need of creating a file in the fasta format: >1A6A.A HVIIQAEFYLNPDQSGEFMFDFDGDEIFHVDMAKKETVWRLEEFGRFASFEAQGALANIAVDKANLEIMTKRSNYTPITN VPPEVTVLTNSPVELREPNVLICFIDKFTPPVVNVTWLRNGKPVTTGVSETVFLPREDHLFRKFHYLPFLPSTEDVYDCR VEHWGLDEPLLKHWEF >1A6A.B ... (5 Replies)
Discussion started by: lost
5 Replies

7. UNIX for Dummies Questions & Answers

Changing from Excel date format to MySQL date format

I have a list of dates in the following format: mm/dd/yyyy and want to change these to the MySQL standard format: yyyy-mm-dd. The dates in the original file may or may not be zero padded, so April is sometimes "04" and other times simply "4". This is what I use to change the format: sed -i '' -e... (2 Replies)
Discussion started by: figaro
2 Replies

8. Shell Programming and Scripting

Changing the date format

Hi all, I have a file with below data af23b|11-FEB-12|acc7 ad23b|12-JAN-12|acc4 as23b|15-DEC-11|acc5 z123b|18-FEB-12|acc1 I need the output as below:-(date in yyyymmdd format) af23b|20120211|acc7 ad23b|20120112|acc4 as23b|20111215|acc5 z123b|20120218|acc1 Please help me on this.... (7 Replies)
Discussion started by: gani_85
7 Replies

9. Shell Programming and Scripting

Shell script for changing the accession number of DNA sequences in a FASTA file

Hi, I am having a file of dna sequences in fasta format which look like this: >admin_1_45 atatagcaga >admin_1_46 atatagcagaatatatat with many such thousands of sequences in a single file. I want to the replace the accession Id "admin_1_45" similarly in following sequences to... (5 Replies)
Discussion started by: margarita
5 Replies

10. UNIX for Dummies Questions & Answers

Changing the file name format

Hello all, I am tryign to change the format of files (which are many in numbers). They at present are named like this: SomeProcess_M-130_100_1_3BR.root SomeProcess_M-130_101_2_3BX.root SomeProcess_M-130_103_3_3RY.root SomeProcess_M-130_105_1_3GH.root SomeProcess_M-130_99_1_3LF.root... (7 Replies)
Discussion started by: emily
7 Replies

LEARN ABOUT DEBIAN

glam2-purge

GLAM2-PURGE(1)							   glam2 Manual 						    GLAM2-PURGE(1)

NAME

       glam2-purge - Removes redundant sequences from a FASTA file

SYNOPSIS

       glam2-purge file score [options]

DESCRIPTION

       glam2-purge is a modified version of Andrew Neuwald's purge program that removes redundant sequences from a FASTA file. This is recommended
       in order to prevent highly similar sequences distorting the search for motifs. Purge works with either DNA or protein sequences and creates
       an output file such that no two sequences have a (gapless) local alignment score greater than a threshold specified by the user. The output
       file is named <file>.<score>. The alignment score is based on the BLOSUM62 matrix for proteins, and on a +5/-1 scoring scheme for DNA.
       Purge can also be used to mask tandem repeats. It uses the XNU program for this purpose.

OPTIONS

       -n
	   Sequences are DNA (default: protein).

       -b
	   Use blast heuristic method (default for protein).

       -e
	   Use an exhaustive method (default for DNA).

       -q
	   Keep first sequence in the set.

       -x
	   Use xnu to mask protein tandem repeats.

SEE ALSO

       glam2(1), glam2format(1), glam2mask(1), glam2scan(1), xnu(1)

       The full Hypertext documentation of GLAM2 is available online at http://bioinformatics.org.au/glam2/ or on this computer in
       /usr/share/doc/glam2/.

REFERENCES

       Purge was written by Andy Neuwald and is described in more detail in Neuwald et al., "Gibbs motif sampling: detection of bacterial outer
       membrane protein repeats", Protein Science, 4:1618-1632, 1995. Please cite it if you use Purge.

       If you use GLAM2, please cite: MC Frith, NFW Saunders, B Kobe, TL Bailey (2008) Discovering sequence motifs with arbitrary insertions and
       deletions, PLoS Computational Biology (in press).

AUTHORS

       Andrew Neuwald
	   Author of purge, renamed glam2-purge in Debian.

       Martin Frith
	   Modified purge to be ANSI standard C and improved the user interface.

       Timothy Bailey
	   Modified purge to be ANSI standard C and improved the user interface.

       Charles Plessy <plessy@debian.org>
	   Formatted this manpage in DocBook XML for the Debian distribution.

COPYRIGHT

       The source code and the documentation of Purge and GLAM2 are released in the public domain.

GLAM2 1056							    05/19/2008							    GLAM2-PURGE(1)
All times are GMT -4. The time now is 03:47 AM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy