BP_MUTATE(1p) User Contributed Perl Documentation BP_MUTATE(1p)NAME
mutate.pl - randomly mutagenize a single protein or DNA sequence
SYNOPSIS
./mutate.pl -p 25 -i test.fa -n 5 -f swiss -o muts.swiss
#or
./mutate.pl --percent=25 --input=test.fa --number=5 -output=x.fa
DESCRIPTION
Randomly mutagenize a single protein or DNA sequence one or more times. Specify percentage mutated and number of resulting mutant
sequences. Print mutagenized sequences to STDOUT or write to an output file.
-h|--help Help
-p|--percent Percent mutagenized
-n|--number Number of mutant sequences created
-o|--output Output file (optional)
-f|--format Output format (default: fasta)
-i|--input Input file
FEEDBACK
User feedback is an integral part of the evolution of this and other Bioperl scripts. Send your comments and suggestions 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
Brian Osborne, bosborne at alum.mit.edu
perl v5.14.2 2012-03-02 BP_MUTATE(1p)
Check Out this Related Man Page
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)
Hi,
I have string like this:
$str=' DNA OR ("rna AND binding AND protein")';
I just wanted to substitute AND with a blank. How can i do that?
I want the output like this:
$string= DNA OR ("rna binding protein") (3 Replies)
Hi,
I'm working with DNA sequences, which as you might know are composed of "codons", which are the "words" if you like, and all codons are 3 letters long. So something I frequently need in my scripts, is pattern searching that respects this typical word size.
E.g. : searching "TAC" in the... (6 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)
I am trying to reverse and complement my DNA sequences. The file format is FASTA, something like this:
Now, to reverse the sequence, I should start reading from right to left. At the same should be complemented. Thus, "A" should be read as "T"; "C" should be read as "G"; "T" should be converted... (8 Replies)
Hi everyone,
I have a large text file containing DNA sequences in fasta format as follows:
>someseq
GAACTTGAGATCCGGGGAGCAGTGGATCTC
CACCAGCGGCCAGAACTGGTGCACCTCCAG
GCCAGCCTCGTCCTGCGTGTC
>another seq
GGCATTTTTGTGTAATTTTTGGCTGGATGAGGT
GACATTTTCATTACTACCATTTTGGAGTACA
>seq3450... (4 Replies)
Hi All, here's a question from newbie
I have a data like this, which set of small DNA sequences separated by new line
GAATCCGGAAACAGCAACTTCAAANCA
GTNATTCGGGCCAAACTGTCGAA
TTNGGCAACTGTTAGAGCTCATGCGACA
CCTGCTAAACGAGTTCGAGTTGAANGA
TTNCGGAAGTGGTCGCTGGCACGG
ACNTGCATGTACGGAGTGACGAAACCI... (6 Replies)
I have two files containing hundreds of different sequences with the same Identifiers (ID-001, ID-002, etc.,), something like this:
Infile1:
ID-001 ATGGGAGCGGGGGCGTCTGCCTTGAGGGGAGAGAAGCTAGATACA
ID-002 ATGGGAGCGGGGGCGTCTGTTTTGAGGGGAGAGAAGCTAGATACA
ID-003... (18 Replies)
Hi,
I have two files: a.doc and b.txt
I wish to search the strings from file b.txt in a.doc and want to highlight them in a.doc with different colours using Perl or bash./awk/sed?
Please guide me. :)
Thanks!!!!! (10 Replies)
Hi Tech Guru,
I have a test file as below , which needs some more fields to be populated randomly :
dks3243;12;20130823;1420;25m;0;syt
dks3243;rocy;10
dks3243;kiop;18
sde21p4;77;20151210;8479;7py;9;vfr
sde21p4;temp;67
sfq6i01;12;20120123;3412;4rd;7;jui
sfq6i01;uymk;90
sfq6i01;kiop;51
... (8 Replies)
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)
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,
Attached is my very simple C++ code to remove any substrings (DNA sequence) of each other, i.e. any redundant sequence is removed to get unique sequences. Similar to sort | uniq command except there is reverse-complementary for DNA sequence. The program runs well with small dataset, but... (11 Replies)
I want to reverse some DNA sequences and complement them at the same time. Thus, A changes to T; C to G; T to A and G to C.
example:
infile
>GHL8OVD01CMQVT SHORT1
TTGATGT
>GHL8OVD01CMQVT SHORT2
TTGATGT
outfile:
>GHL8OVD01CMQVT SHORT1
ACATCAA
>GHL8OVD01CMQVT SHORT2
ACATCAA
The... (4 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)