10 likes | 168 Views
Create a program that reads a DNA sequence from the keyboard, and reports the sequence length and the G/C content of the sequence (as a fraction) Adapt the G/C script so multiple sequences in FASTA format are read from a file
E N D
Create a program that reads a DNA sequence from the keyboard, and reports the sequence length and the G/C content of the sequence (as a fraction) • Adapt the G/C script so multiple sequences in FASTA format are read from a file • Modify the script to process a file containing any number of sequences in EMBL format • Now let the program generate the reverse complement of the sequence(s), and report sequence length and G/C content • Using the rand function of Perl to shuffle the nucleotides of the input sequence, while maintaining sequence composition; again report sequence length and G/C content • Create a script to find the DNA fragments you get after cutting the sequence in the example1.fasta file with AluI and with AvaI • Find the open reading frames in the example1.fasta sequence • Translate the open reading frames to protein, using the standard genetic code from the Geneticcode database (http://srs.bioinformatics.nl)