Sponsored Content
Top Forums Shell Programming and Scripting Extracting DNA sequences from GenBank files using Perl Post 302328173 by akreibich07 on Tuesday 23rd of June 2009 03:28:18 PM
Old 06-23-2009
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 the genes and their positions in the species corresponding to NC_001666 (corn). Here is a sample of how the text file is formatted...

rbcL (56874..58304)
atpB -(54618..56114)
atpE -(54208..54621)
trnM (54020..54092)
trnV -(53158..53834)

For example, if in my command prompt I give input of the program name, the species number that I want, and the specific gene from that species whose DNA sequence I want:

perl nucleotide_bases.pl NC_001666 trnM

The program would go into NC_001666.txt, find trnM, see that it has a range from 54020 to 54092 and is on the positive strand(no negative sign). The program then goes into NC_001666.gb, goes to the long list of DNA bases at the bottom and starts at position 54020 and returns all base letters through 54092 (inclusively). So for this specific trnM, the output would be:

gcctacttaactcagtggttagagtattgctttcatacggcgggagtcattggttcaaatccaatagtaggta

If a gene has a negative next to the position range (meaning it's on the negative strand of DNA), the output should be reversed, starting from the higher position, going to the lower. Also, when a negative is there, in that output, all A's should be switched to T's, and all G's to C's and vice versa.

Also, if a gene appears more than once in a text file, give an error message that it appears more than once, and end the program.

If I could get a Perl script to return this information for any species (NC_number) I want, and any gene from that species that I want, it would be a great help in the research I am conducting. Thank you all for your time, and any help on how to write this script would be appreciated.

-akreibich07
 

10 More Discussions You Might Find Interesting

1. Shell Programming and Scripting

Perl - extracting data from .csv files

PROJECT: Extracting data from an employee timesheet. The timesheets are done in excel (for user ease) and then converted to .csv files that look like this (see color code key below): ,,,,,,,,,,,,,,,,,,, 9/14/2003,<-- Week Ending,,,,,,,,,,,,,,,,,, Craig Brennan,,,,,,,,,,,,,,,,,,,... (3 Replies)
Discussion started by: kregh99
3 Replies

2. Shell Programming and Scripting

GenBank Perl help...

Hey guys, I'm doing some Perl scripting for genomic data out of GenBank files...I have to extract the name of the plant, the file name, the number of bases, and all of the genes including their starting and ending positions...for example, with this GenBank file, LOCUS NC_010093 ... (7 Replies)
Discussion started by: akreibich07
7 Replies

3. Shell Programming and Scripting

To retrieve gene name & function thorugh Genbank id (gi|9910297)

Hi , I have list of genbank id's and ref number in this format. gi|9910297|ref|NM_019974.1| I want to retrive the gene name and fuction for each genbank list. I have around 1300 gi numbers in my excel sheet. So anybody can help me to retrive the information from NCBI through perl script... (0 Replies)
Discussion started by: shibujohn82
0 Replies

4. UNIX for Advanced & Expert Users

Extracting files with multiple links-perl

i want to write a perl script that gets/displays all those files having multiple links (in current directory) (4 Replies)
Discussion started by: guptesanket
4 Replies

5. Shell Programming and Scripting

Extracting column value from perl

Hello Kindly help me to find out the first column from first line of a flat file in perl I/P 9869912|20110830|00000000000013009|130|09|10/15/2010 12:36:22|W860944|N|00 9869912|20110830|00000000000013013|130|13|10/15/2010 12:36:22|W860944|N|00... (5 Replies)
Discussion started by: Pratik4891
5 Replies

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

7. Shell Programming and Scripting

Randomly selecting sequences and generating specific output files

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)
Discussion started by: Xterra
18 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. Shell Programming and Scripting

Extraction of sequences from files

hey!!! I have 2 files file1 is as ids.txt and is >gi|546473186|gb|AWWX01630222.1| >gi|546473233|gb|AWWX01630175.1| >gi|546473323|gb|AWWX01630097.1| >gi|546474044|gb|AWWX01629456.1| >gi|546474165|gb|AWWX01629352.1| file2 is sequences.fasta and is like >gi|546473233|gb|AWWX01630175.1|... (9 Replies)
Discussion started by: harpreetmanku04
9 Replies

10. Shell Programming and Scripting

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... (3 Replies)
Discussion started by: faizlo
3 Replies
Bio::SeqIO::Handler::GenericRichSeqHandler(3pm) 	User Contributed Perl Documentation	   Bio::SeqIO::Handler::GenericRichSeqHandler(3pm)

NAME
Bio::SeqIO::Handler::GenericRichSeqHandler - Bio::HandlerI-based data handler for GenBank/EMBL/UniProt (and other) sequence data SYNOPSIS
# MyHandler is a GenericRichSeqHandler object. # inside a parser (driver) constructor.... $self->seq_handler($handler || MyHandler->new(-format => 'genbank')); # in next_seq() in driver... $hobj = $self->seqhandler(); # roll data up into hashref chunks, pass off into Handler for processing... $hobj->data_handler($data); # or retrieve Handler methods and pass data directly to Handler methods... my $hmeth = $hobj->handler_methods; if ($hmeth->{ $data->{NAME} }) { my $mth = $hmeth->{ $data->{NAME} }; $hobj->$mth($data); } DESCRIPTION
This is an experimental implementation of a sequence-based HandlerBaseI parser and may change over time. It is possible (nay, likely) that the way handler methods are set up will change over development to allow more flexibility. Release pumpkins, please do not add this to a release until the API has settled. It is also likely that write_seq() will not work properly for some data. Standard Developer caveats: Do not use for production purposes. Not responsible for destroying (your data|computer|world). Do not stare directly at GenericRichSeqHandler. If GenericRichSeqHandler glows, back slowly away and call for help. Consider yourself warned! This class acts as a demonstration on how to handle similar data chunks derived from Bio::SeqIO::gbdriver, Bio::SeqIO::embldriver, and Bio::SeqIO::swissdriver using similar (or the same) handler methods. The modules currently pass all previous tests in t/genbank.t, t/embl.t, and t/swiss.t yet all use the same handler methods (the collected tests for handlers can be found in t/Handler.t). Some tweaking of the methods themselves is probably in order over the long run to ensure that data is consistently handled for each parser. Round-trip tests are probably in order here... Though a Bio::Seq::SeqBuilder is employed for building sequence objects no bypassing of data based on builder slots has been implemented (yet); this is planned in the near future. As a reminder: this is the current Annotation data chunk (via Data::Dumper): $VAR1 = { 'NAME' => 'REFERENCE', 'DATA' => '1 (bases 1 to 10001)' 'AUTHORS' => 'International Human Genome Sequencing Consortium.' 'TITLE' => 'The DNA sequence of Homo sapiens' 'JOURNAL' => 'Unpublished(2003)' }; ... This is the current SeqFeature data chunk (again via Data::Dumper): $VAR1 = { 'mol_type' => 'genomic DNA', 'LOCATION' => '<1..>10001', 'NAME' => 'FEATURES', 'FEATURE_KEY' => 'source', 'note' => 'Accession AL451081 sequenced by The Sanger Centre', 'db_xref' => 'taxon:9606', 'clone' => 'RP11-302I18', 'organism' => 'Homo sapiens' }; 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/ AUTHOR - Chris Fields Email cjfields at bioperl dot org APPENDIX
The rest of the documentation details each of the object methods. Internal methods are usually preceded with a _ new Title : new Usage : Function: Returns : Args : -format Sequence format to be mapped for handler methods -builder Bio::Seq::SeqBuilder object (normally defined in SequenceStreamI object implementation constructor) Throws : On undefined '-format' sequence format parameter Note : Still under heavy development Bio::HandlerBaseI implementing methods handler_methods Title : handler_methods Usage : $handler->handler_methods('GenBank') %handlers = $handler->handler_methods(); Function: Retrieve the handler methods used for the current format() in the handler. This assumes the handler methods are already described in the HandlerI-implementing class. Returns : a hash reference with the data type handled and the code ref associated with it. Args : [optional] String representing the sequence format. If set here this will also set sequence_format() Throws : On unimplemented sequence format in %HANDLERS data_handler Title : data_handler Usage : $handler->data_handler($data) Function: Centralized method which accepts all data chunks, then distributes to the appropriate methods for processing based on the chunk name from within the HandlerBaseI object. One can also use Returns : None Args : an hash ref containing a data chunk. reset_parameters Title : reset_parameters Usage : $handler->reset_parameters() Function: Resets the internal cache of data (normally object parameters for a builder or factory) Returns : None Args : None format Title : format Usage : $handler->format('GenBank') Function: Get/Set the format for the report/record being parsed. This can be used to set handlers in classes which are capable of processing similar data chunks from multiple driver modules. Returns : String with the sequence format Args : [optional] String with the sequence format Note : The format may be used to set the handlers (as in the current GenericRichSeqHandler implementation) get_params Title : get_params Usage : $handler->get_params('-species') Function: Convenience method used to retrieve the specified parameters from the internal parameter cache Returns : Hash ref containing parameters requested and data as key-value pairs. Note that some parameter values may be objects, arrays, etc. Args : List (array) representing the parameters requested set_params Title : set_params Usage : $handler->set_param({'-species') Function: Convenience method used to set specific parameters Returns : None Args : Hash ref containing the data to be passed as key-value pairs Methods unique to this implementation seqbuilder Title : seqbuilder Usage : Function: Returns : Args : Throws : Note : build_sequence Title : build_sequence Usage : Function: Returns : Args : Throws : Note : location_factory Title : location_factory Usage : Function: Returns : Args : Throws : Note : annotation_collection Title : annotation_collection Usage : Function: Returns : Args : Throws : Note : perl v5.14.2 2012-03-02 Bio::SeqIO::Handler::GenericRichSeqHandler(3pm)
All times are GMT -4. The time now is 05:48 AM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy