03-25-2013
Both of these work, thank you.
The problem with my original code was obscured by my choice of an example file with unrealistically short sequences. In reality my sequences are longer, and the problem with my code is that it captures headers along with short sequences.
10 More Discussions You Might Find Interesting
1. UNIX for Dummies Questions & Answers
Hi ,
I have a peculiar case, where my sed command is working on a file which contains lines of small length.
sed "s/XYZ:1/XYZ:3/g" abc.txt > xyz.txt
when abc.txt contains lines of small length(currently around 80 chars) , this sed command is working fine.
when abc.txt contains lines of... (3 Replies)
Discussion started by: thanuman
3 Replies
2. UNIX for Dummies Questions & Answers
Hi Guys,
Could anyone shed some light on the length of a Unix filepath max length pls ?
thanks !
Wilson (3 Replies)
Discussion started by: wilsontan
3 Replies
3. Shell Programming and Scripting
Hi all,
I have a flat file of 1000 rows. I want to check the length of the 5th column. The one having the longest length , I want to set it as DEFINED PARAMETER.
So later I can check others with that particular number only.
Any ideas ?? (2 Replies)
Discussion started by: ganesh123
2 Replies
4. UNIX for Advanced & Expert Users
Hi,
This is my first post to this site. So kindly forgive if I am writing in a wrong section.
My query is that...
I want to modify the max username length size. I guess it is 32/64 on CentOS. Now I want to change it to 128. Is there any way to do that?
Thanks in advance!! :) (4 Replies)
Discussion started by: ajay303
4 Replies
5. UNIX for Dummies Questions & Answers
Hey Any one...
Do u know any way I can modify the max username length in unix? I guess it is 32/64 characters by default. Suppose I want to increase it to 128.
i hav tried /etc/skel
but no use...
How can I do that? (2 Replies)
Discussion started by: MayureshRisbud
2 Replies
6. Shell Programming and Scripting
I need a way to extract data from
X 4T Solution 21 OCT 2011 37 .00
to account 12345678 User1 user2
X 4T Solution Solution Unlimited 11 Sep 2009 248 .00
to account 87654321 user3 user4
I need it to extract 'X' '37.00' and account number 12345678.
I have extracted above stuff... (3 Replies)
Discussion started by: chakrapani
3 Replies
7. Shell Programming and Scripting
I have a file that stores data in pairs of lines, following this format:
line 1: header (preceded by ">")
line 2: sequence
Example.txt:
>seq1 name
GATTGATGTTTGAGTTTTGGTTTTT
>seq2 name
TTTTCTTC
I want to filter out the sequences and corresponding headers for all sequences that are less... (2 Replies)
Discussion started by: pathunkathunk
2 Replies
8. UNIX for Dummies Questions & Answers
Hi all,
I'm struggling to filter my data frame. I need to print only those lines whose max value (the number of columns may vary) is above a cut-off value.
My data looks like this:
chr22 17565753 17565754 5 4 5 5 6 2 5 5 6 2
chr22 17565754 ... (2 Replies)
Discussion started by: lsantome
2 Replies
9. Shell Programming and Scripting
Hi
Can someone tell me which script will work best (in terms of speed and simplicity to write and run) for a large text file to filter all the lines with a minimum specified length of words ? A sample script with be definitely of great help !!! Thanks in advance. :) (4 Replies)
Discussion started by: my_Perl
4 Replies
10. Shell Programming and Scripting
Hello Team,
Need your expertise on following:
Here is the set of data:
C1|4|C1SP1|A1|C1BP1|T1
C1|4|C1SP2|A1|C1BP2|T2
C2|3|C2SP1|A2|C2BP1|T2
C3|3|C3SP1|A3|C3BP1|T2
C2|2|C2SP2|A2|C2BP2|T1
I need to filter above date base on following two steps:
1. Group them by column 1 and 4
2.... (12 Replies)
Discussion started by: angshuman
12 Replies
LEARN ABOUT DEBIAN
alifmt
ALIFMT(5) User Manuals ALIFMT(5)
NAME
alifmt - Aligned sequences formats
DESCRIPTION
This document illustrates some common formats used for aligned sequences representation.
CLUSTAL
CLUSTAL W (1.82) multiple sequence alignment
MALK_ECOLI MASVQLQNVTKAWGEVVVSKDINLDIHEGEFVVFVGPSGCGKSTL
MALK_SALTY MASVQLRNVTKAWGDVVVSKDINLDIHDGEFVVFVGPSGCGKSTL
MALK_ENTAE MASVQLRNVTKAWGDVVVSKDINLEIQDGEFVVFVGPSGCGKSTL
MALK_PHOLU MSSVTLRNVSKAYGETIISKNINLEIQEGEF--------------
*:** *:**:**:*:.::**:***:*::***
MALK_ECOLI LRMIAGLETITSGDLACRRLHKEPGV
MALK_SALTY LRMIAGLETITSGDLACRRLHQEPGV
MALK_ENTAE LRMIAGLETVTSGDL-----------
MALK_PHOLU LRM-----------------------
***
Warning: Names must not contain spaces or exceed 30 characters.
FASTA
>MALK_ECOLI
MASVQLQNVTKAWGEVVVSKDINLDIHEGEFVVFVGPSGCGKSTLLRMIA
GLETITSGDLACRRLHKEPGV
>MALK_SALTY
MASVQLRNVTKAWGDVVVSKDINLDIHDGEFVVFVGPSGCGKSTLLRMIA
GLETITSGDLACRRLHQEPGV
>MALK_ENTAE
MASVQLRNVTKAWGDVVVSKDINLEIQDGEFVVFVGPSGCGKSTLLRMIA
GLETVTSGDL-----------
>MALK_PHOLU
MSSVTLRNVSKAYGETIISKNINLEIQEGEF--------------LRM--
---------------------
MEGA
#mega
!Title Multiple Sequence Alignment;
#MALK_ECOLI
MASVQLQNVTKAWGEVVVSKDINLDIHEGEFVVFVGPSGCGKSTLLRMIA
GLETITSGDLACRRLHKEPGV
#MALK_SALTY
MASVQLRNVTKAWGDVVVSKDINLDIHDGEFVVFVGPSGCGKSTLLRMIA
GLETITSGDLACRRLHQEPGV
#MALK_ENTAE
MASVQLRNVTKAWGDVVVSKDINLEIQDGEFVVFVGPSGCGKSTLLRMIA
GLETVTSGDL-----------
#MALK_PHOLU
MSSVTLRNVSKAYGETIISKNINLEIQEGEF-------------------
---------------------
MSF
!!AA_MULTIPLE_ALIGNMENT 1.0
PileUp of: @pep.list
msf.seq MSF: 55 Type: P Nov 22, 2001 11:02 Check: 2529 ..
Name: m73237 Len: 655 Check: 7493 Weight: 1.00
Name: l28824 Len: 655 Check: 5456 Weight: 1.00
Name: u04379 Len: 655 Check: 9580 Weight: 1.00
//
1 50
m73237 ~~~~~MADSA NHLPFFFGQI TREEAEDYLV QGGMSDGLYL LRQSRNYLGG
l28824 MASSGMADSA NHLPFFFGNI TREEAEDYLV QGGMSDGLYL LRQSRNYLGG
u04379 ~~~~~MPDPA AHLPFFYGSI SRAEAEEHLK LAGMADGLFL LRQCLRSLGG
51
m73237 ~~~~~
l28824 ~~~~~
u04379 AACG*
Warning: This format cannot handle more than 500 sequences in a single alignment.
NEXUS
#NEXUS
begin data;
dimensions ntax=2 nchar=89;
format datatype=Protein interleave gap=- missing='.';
matrix
[ 1 50]
btdDm -------AQQQQHHLHMQQAQHH-----------LHLSH------QQAQQ
TSp1 --------------------AEH-----------PSLR--------GTPL
[ 51 80]
btdDm YACPICSKKFSRSDHLSKHKKTHF------
TSp1 FACPICNKRFMRSDHLAKHVKTHN------
;
endblock;
Warning: Text enclosed in brackets is considered as comment, and thus ignored.
PHYLIP
Sequential (PHYLIPS):
2 109
ATISA1 GSPNLYQ-GGRKPWHSINFICAHDGFTLADLVTYNNKNNLANGEENNDG
ENHNYSWNCGEEGDFASISVKRLRKRQMRNFFVSLMVSQGVPMIYMGDE
YGHTKGGN---
POTISA1 GSPNLYQKGGRKPWNSINFVCAHDGFTLADLVTYNNKHNLANGEDNKDG
ENHNNSWNCGEEGEFASIFVKKLRKRQMRNFFLCLMVSQGVPMIYMGDE
YGHTKGGN---
Interleaved (PHYLIPI):
2 109
ATISA1 GSPNLYQ-GGRKPWHSINFICAHDGFTLADLVTYNNKNNLANGEENND
POTISA1 GSPNLYQKGGRKPWNSINFVCAHDGFTLADLVTYNNKHNLANGEDNKD
GENHNYSWNCGEEGDFASISVKRLRKRQMRNFFVSLMVSQGVPMIYMG
GENHNNSWNCGEEGEFASIFVKKLRKRQMRNFFLCLMVSQGVPMIYMG
DEYGHTKGGN---
DEYGHTKGGN---
Warning: Species names may not contain characters `( ) : ; , [ ]' and exceed 10 characters.
STOCKHOLM
# STOCKHOLM 1.0
MALK_ECOLI MASVQLQNVTKAWGEVVVSKDINLDIHEGEFVVFVGPSGCGKSTLLRMIA
MALK_SALTY MASVQLRNVTKAWGDVVVSKDINLDIHDGEFVVFVGPSGCGKSTLLRMIA
MALK_ENTAE MASVQLRNVTKAWGDVVVSKDINLEIQDGEFVVFVGPSGCGKSTLLRMIA
MALK_PHOLU MSSVTLRNVSKAYGETIISKNINLEIQEGEF-------------------
MALK_ECOLI RCHLFREDGTACRRLHKEPGV
MALK_SALTY RCHLFREDGSACRRLHQEPGV
MALK_ENTAE ---------------------
MALK_PHOLU ---------------------
//
SEE ALSO
squizz(1), seqfmt(5)
AUTHOR
Nicolas Joly (njoly@pasteur.fr), Institut Pasteur.
Unix 2009-05-19 ALIFMT(5)