Sponsored Content
Top Forums Shell Programming and Scripting Convert a DNA sequence into Amino Acid Post 302957292 by faizlo on Thursday 8th of October 2015 05:12:58 PM
Old 10-08-2015
Convert a DNA sequence into Amino Acid

I am trying to write a bash script that would be able to read DNA sequences (each line in the file is a sequence) from a file, where sequences are separated by an empty line. I am then to find the amino acid that these DNA sequences encode per codon (each group of three literals.) For example, if I have a file with the sequence:

Code:
   GCATGCTGCGATAACTTTGGCTGAACTTTGGCTGAAGCATGCTGCGAAACTTTGGCTGAACTTTGGCTG

then starting from GCA (first three literals,) I want to decode the DNA into amino acids based on the following table:

Code:
    Codon(s)                  Amino-acid
    TTT,TTC                    Phe
    TTA,TTG,CTT,CTC,CTA,CTG    Leu
    ATT,ATC,ATA                Ile
    ATG                       Met
    GTT,GTC,GTA,GTG            Val
    TCT,TCC,TCA,TCG            Ser
    CCT,CCC,CCA,CCG            Pro
    ACT,ACC,ACA,ACG            Thr
    GCT,GCC,GCA,GCG            Ala
    TAT,TAC                    Tyr
    TAA,TAG                    Stop
    CAT,CAC                    His
    CAA,CAG                    Gln
    AAT,AAC                    Asn
    AAA,AAG                    Lys
    GAT,GAC                   Asp
    GAA,GAG                   Glu
    TGT,TGC                    Cys
    TGA                        Stop
    TGG                        Trp
    CGT,CGC,CGA,CGG            Arg
    AGT,AGC                    Ser
    AGA,AGG                    Arg
    GGT,GGC,GGA,GGG            Gly

that is, I need to get:
Code:
    AlaCysCysAspAsnPheGlyStopThrLeuAlaGluAlaCysCysGluThrLeuAlaGluLeuTrpLeu

Then I need to print the name of each Amino Acid and how many times it was used. For example:

Code:
    Ala: 4
    Cys: 4

and so on. I have 100s of files with DNA sequences in them, but I am not that good at bash. I tried awk and tr but I did not know how to code the table into a bash script.

Last edited by faizlo; 10-08-2015 at 09:08 PM.. Reason: Add CODE tags.
 

9 More Discussions You Might Find Interesting

1. Shell Programming and Scripting

escape sequence for $

Hi all, I have a requirement where the variable name starts with $, like $Amd=/home/student/test/ How to work wit it? can some one help me, am in gr8 confusion:confused: (5 Replies)
Discussion started by: shreekrishnagd
5 Replies

2. Shell Programming and Scripting

How to remove those sequence with same amino acid?What command line I should type?

My input is listed as: giNumber RefAminoAcid VarAminoAcid 10190711 P P 10190711 D D 109255248 I A 110349771 A ... (4 Replies)
Discussion started by: patrick chia
4 Replies

3. Shell Programming and Scripting

Extracting DNA sequences from GenBank files using Perl

Hi all, Using Perl, I need to extract DNA bases from a GenBank file for a given plant species. A sample GenBank file is here... Nucleotide This is saved on my computer as NC_001666.gb. I also have a file that is saved on my computer as NC_001666.txt. This text file has a list of all... (5 Replies)
Discussion started by: akreibich07
5 Replies

4. Shell Programming and Scripting

Tricky task with DNA sequences.

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)
Discussion started by: Xterra
8 Replies

5. Shell Programming and Scripting

find common entries and match the number with long sequence and cut that sequence in output

Hi all, I have a file like this ID 3BP5L_HUMAN Reviewed; 393 AA. AC Q7L8J4; Q96FI5; Q9BQH8; Q9C0E3; DT 05-FEB-2008, integrated into UniProtKB/Swiss-Prot. DT 05-JUL-2004, sequence version 1. DT 05-SEP-2012, entry version 71. FT COILED 59 140 ... (1 Reply)
Discussion started by: manigrover
1 Replies

6. Shell Programming and Scripting

How to convert multiple number ranges into sequence?

Looking for a simple way to convert ranges to a numerical sequence that would assign the original value of the range to the individual numbers that are on the range. Thank you given data 13196-13199 0 13200 4 13201 10 13202-13207 3 13208-13210 7 desired... (3 Replies)
Discussion started by: jcue25
3 Replies

7. Shell Programming and Scripting

Sequence generator

Thanks Guys This really helped (5 Replies)
Discussion started by: robert89
5 Replies

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

9. Red Hat

Rm -rf * sequence

If I run rm -rf * command under one parent directory. /data > rm -rf * Is there anyway to know which files will be deleted first ? Start using code tags please, ty. (2 Replies)
Discussion started by: sameermohite
2 Replies
Bio::Tools::Analysis::DNA::ESEfinder(3pm)		User Contributed Perl Documentation		 Bio::Tools::Analysis::DNA::ESEfinder(3pm)

NAME
Bio::Tools::Analysis::DNA::ESEfinder - a wrapper around ESEfinder server SYNOPSIS
use Bio::Tools::Analysis::DNA::ESEfinder; use strict; my $seq; # a Bio::PrimarySeqI or Bio::SeqI object $seq = Bio::Seq->new( -primary_id => 'test', -seq=>'atgcatgctaggtgtgtgttttgtgggttgtactagctagtgat'. -alphabet=>'dna'); my $ese_finder = Bio::Tools::Analysis::DNA::ESEfinder-> new(-seq => $seq); # run ESEfinder prediction on a DNA sequence $ese_finder->run(); die "Could not get a result" unless $ese_finder->status =~ /^COMPLETED/; print $ese_finder->result; # print raw prediction to STDOUT foreach my $feat ( $ese_finder->result('Bio::SeqFeatureI') ) { # do something to SeqFeature # e.g. print as GFF print $feat->gff_string, " "; # or store within the sequence - if it is a Bio::SeqI $seq->add_SeqFeature($feat) } DESCRIPTION
This class is a wrapper around the ESEfinder web server which uses experimentally defined scoring matrices to identify possible exonic splicing enhancers in human transcripts. The results can be retrieved in 4 ways. 1. "$ese_finder->result('')" retrieves the raw text output of the program 2. "$ese_finder->result('all')" returns a Bio::Seq::Meta::Array object with prediction scores for all residues in the sequence 3. "$ese_finder->result('Bio::SeqFeatureI')" returns an array of Bio::SeqFeature objects for sequences with significant scores. Feature tags are score, motif, SR_protein and method 4. "$ese_finder->result('raw')" returns an array of significant matches with each element being a reference to [SR_protein, position, motif, score] See <http://rulai.cshl.edu/tools/ESE2/> This the second implentation of Bio::SimpleAnalysisI which hopefully will make it easier to write wrappers on various services. This class uses a web resource and therefore inherits from Bio::WebAgent. SEE ALSO
Bio::SimpleAnalysisI, Bio::WebAgent 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 one of the Bioperl mailing lists. 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 the bugs and their resolution. Bug reports can be submitted via the web: https://redmine.open-bio.org/projects/bioperl/ AUTHORS
Richard Adams, Richard.Adams@ed.ac.uk, Heikki Lehvaslaiho, heikki-at-bioperl-dot-org APPENDIX
The rest of the documentation details each of the object methods. Internal methods are usually preceded with a _ perl v5.14.2 2012-03-02 Bio::Tools::Analysis::DNA::ESEfinder(3pm)
All times are GMT -4. The time now is 04:19 PM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy