BP_AACOMP(1p) User Contributed Perl Documentation BP_AACOMP(1p)NAME
aacomp - amino acid composition of protein sequences
SYNOPSIS
aacomp [-f/--format FORMAT] [-h/--help] filename
or
aacomp [-f/--format FORMAT] < filename
or
aacomp [-f/--format FORMAT] -i filename
DESCRIPTION
This scripts prints out the count of amino acids over all protein sequences from the input file.
OPTIONS
The default sequence format is fasta.
The sequence input can be provided using any of the three methods:
unnamed argument
aacomp filename
named argument
aacomp -i filename
standard input
aacomp < filename
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
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:
https://redmine.open-bio.org/projects/bioperl/
AUTHOR - Jason Stajich
Email jason@bioperl.org
HISTORY
Based on aacomp.c from an old version of EMBOSS
perl v5.14.2 2012-03-02 BP_AACOMP(1p)
Check Out this Related Man Page
BP_CHAOS_PLOT(1p) User Contributed Perl Documentation BP_CHAOS_PLOT(1p)NAME
chaos_plot - a chaos plot from DNA and RNA sequences
SYNOPSIS
chaos_plot.pl -i/--input=INPUTFILE -f/--format=SEQFORMAT
-o/--output=OUTPUTFILE -g/--graphics=GRAPHIC FORMAT
-w/--width=WIGHT -h/--height=HEIGHT
DESCRIPTION
This scripts generates image files using GD image library to visualize nucleotide sequences using chaos plot.
OPTIONS
Valid graphics formats are currently gd, gd2, png, wbmp, jpeg and gif.
The default size of the image file is 600x400.
The sequence input can be provided using any of the three methods:
unnamed argument
chaos_plot filename
named argument
chaos_plot -i filename
standard input
chaos_plot < filename
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
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:
https://redmine.open-bio.org/projects/bioperl/
AUTHOR - Jason Stajich
Email jason@bioperl.org
HISTORY
This code is based on EMBOSS C code for chaos.c by Ian Longden. Included are documentation from EMBOSS code:
Chaos produces a chaos plot. The original application is part of the ACEDB genome database package, written by ** Richard Durbin (MRC LMB,
UK) rd@mrc-lmba.cam.ac.uk, and Jean Thierry-Mieg (CRBM du CNRS, France) mieg@crbm1.cnusc.fr
perl v5.14.2 2012-03-02 BP_CHAOS_PLOT(1p)
HI,
I have sentences like this:
@str=("An ribonucleic acid (RNA)-binding protein, started its expression in the daughter cells","Elucidation of the mechanism of retinal degeneration of RNA-binding protein","Rna binding protein is an important protein","In the retinal degenerative process... (1 Reply)
KINDLY HELP ME FOR SHELL SCRIPTING FOR THIS TASK.
My input file consists of thousands of sequence in this format. The given input file consists of four sequences which are starting with ‘>’ symbol (each sequence shown in different colour for easy understanding). I have to use a command at $... (3 Replies)
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)
Hi all,
Looking for suggestions on a better way to sum numbers in a key value pair formated file. What I have works but seems really clunky to me. Any suggestions would be greatly appreciated.
cat test.txt | perl -ne 'm/(M=)(\d+\.?\d?\d?)/ && print "$2\n"' | awk '{ sum+=$1} END {printf... (7 Replies)
Hello,
A bioperl problem I thought could be done with awk: convert the fasta format (Note: the length of each row is not the same for each entry as they were combined from different files!) to tabular format.
input.fasta:
>YAL069W-1.334 Putative promoter sequence... (6 Replies)
I have a file that contains a number of fasta sequences, each one starts with '>'. I would like to count the number of sequences I have in the file. When I use the code below I get:
$ grep -c '>' total_Protein.txt
1
even though the file contains hundreds
If I use:
$ grep -c {>}... (3 Replies)
I have a file of protein sequences with headers (my source file). Based on a list of IDs (which are included in some of the headers), I'd like to print out only the specified sequences, with only the ID as header.
In other words, I'd like to search source.txt for the terms in IDs.txt, and print... (3 Replies)
I have a fasta file as follows
>sp|O15090|FABP4_HUMAN Fatty acid-binding protein, adipocyte OS=Homo sapiens GN=FABP4 PE=1 SV=3
MCDAFVGTWKLVSSENFDDYMKEVGVGFATRKVAGMAKPNMIISVNGDVITIKSESTFKN
TEISFILGQEFDEVTADDRKVKSTITLDGGVLVHVQKWDGKSTTIKRKREDDKLVVECVM
KGVTSTRVYERA
>sp|L18484|AP2A2_RAT AP-2... (3 Replies)
Hi
How can I extract sequences from a fasta file with respect a certain criteria? The beginning of my file (containing in total more than 1000 sequences) looks like this:
>H8V34IS02I59VP
SDACNDLTIALLQIAREVRVCNPTFSFRWHPQVKDEVMRECFDCIRQGLG
YPSMRNDPILIANCMNWHGHPLEEARQWVHQACMSPCPSTKHGFQPFRMA... (6 Replies)
Hello,
I am trying to create a file in windows and i want the filename to have timestamp as well but something is wrong and i can not understand waht. The code that i use is the following
($cwkday,$cmonth,$cday,$ctime,$cyear) = split(/\s+/, localtime);
$current_date =... (5 Replies)
Hi,
I have a file which has wrong time format and we want to correct it before we load it.
WRONG FORMAT : 93:0:00
CORRECT FORMAT :09:30:00
If you notice the 0 at the front is missing. Its the case always. (6 Replies)
I have a fasta file as follows
>sp|Q8WWQ8|STAB2_HUMAN Stabilin-2 OS=Homo sapiens OX=9606 GN=STAB2 PE=1 SV=3
MMLQHLVIFCLGLVVQNFCSPAETTGQARRCDRKSLLTIRTECRSCALNLGVKCPDGYTM
ITSGSVGVRDCRYTFEVRTYSLSLPGCRHICRKDYLQPRCCPGRWGPDCIECPGGAGSPC
NGRGSCAEGMEGNGTCSCQEGFGGTACETCADDNLFGPSCSSVCNCVHGVCNSGLDGDGT... (3 Replies)