Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #Vraag 11
- cat "genes_CDS.fa" | tr '\n' ' ' | sed 's/ >/\n>/g' | awk '
- BEGIN {
- print "Het aantal codons per gen dat codeert voor cysteine:"
- }
- {
- tel=0
- for (i=1;i<=length($2);i+=3) {
- if (substr($2,i,3) == "TGT" || substr($2,i,3) "TGC") {
- tel++
- }
- }
- print substr($1,2)":",tel
- }
- '
- #Vraag 12
- cat "Zebrafish_proteins_single_line.txt" | awk '
- {
- tel=0
- for (i=1;i<=length($2);i++) {
- if (substr($2,i,1) == "M") {
- tel++;
- }
- }
- percentage=tel/length($2)*100;
- print substr($1,2)":",tel, "("percentage") %"
- }
- ' | sort -nrk2 | head -5
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement