Bklein7 Week 3
From LMU BioDB 2015
Revision as of 01:10, 16 September 2015 by Bklein7 (Talk | contribs) (Added command sequences for negative reading frames and their outputs)
Complement of a Strand
Write a sequence of piped text processing commands that, when given a nucleotide sequence, returns its complementary strand. In other words, fill in the question marks:
cat sequence_file | sed "y/atcg/tagc/"
Reading Frames
Write 6 sets of text processing commands that, when given a nucleotide sequence, returns the resulting amino acid sequence, one for each possible reading frame for the nucleotide sequence.
Outputs generated using ~dondi/xmlpipedb/data/prokaryote.txt:
- +1 Reading Frame
cat sequence_file | sed "s/t/u/g" | sed "s/.../& /g" | sed -f genetic-code.sed | sed "s/[atcg]//g" | sed "s/ //g" Output: STIFQ-VRWPKKTILNLKRCLIPCSAYNPAASSAGGIL
- +2 Reading Frame
cat sequence_file | sed "s/t/u/g" | sed "s/^.//g" | sed "s/.../& /g" | sed -f genetic-code.sed | sed "s/[atcg]//g" | sed "s/ //g" Output: LLYFNRYDGQRRQY-T-NVA-YHVPRITQPPVPLAAF-
- +3 Reading Frame
cat sequence_file | sed "s/t/u/g" | sed "s/^..//g" | sed "s/.../& /g" | sed -f genetic-code.sed | sed "s/[atcg]//g" | sed "s/ //g" Output: YYISIGTMAKEDNIELETLPNTMFRV-PSRQFRWRHFN
- -1 Reading Frame
cat sequence_file | sed "y/atcg/tagc/" | sed "s/t/u/g" | rev | sed "s/.../& /g" | sed -f genetic-code.sed | sed "s/[atcg]//g" | sed "s/ //g" Output: VKMPPAELAAGLYAEHGIRQRFKENIVFFGHRTY-NIV
- -2 Reading Frame
cat sequence_file | sed "y/atcg/tagc/" | sed "s/t/u/g" | rev | sed "s/^.//g" | sed "s/.../& /g" | sed -f genetic-code.sed | sed "s/[atcg]//g" | sed "s/ //g" Output: LKCRQRNWRLGYTRNMVLGNVSSSILSSLAIVPIEI--
- -3 Reading Frame
cat sequence_file | sed "y/atcg/tagc/" | sed "s/t/u/g" | rev | sed "s/^..//g" | sed "s/.../& /g" | sed -f genetic-code.sed | sed "s/[atcg]//g" | sed "s/ //g" Output: -NAASGTGGWVIRGTWY-ATFQVQYCLLWPSYLLKYSR
Check Your Work
Fortunately, online tools are available for checking your work; we recommend the ExPASy Translate Tool, sponsored by the same people who run SwissProt. You’re free to use this tool to see if your text processing commands produce the same results.