Commit 974c7742 authored by Jeroen F.J. Laros's avatar Jeroen F.J. Laros

Added notes.

parent 21f5616e
AWK wikibook: http://en.wikibooks.org/wiki/An_Awk_Primer
Regular expressions: http://www.youtube.com/watch?v=EppQdkv4G2w
awk 'BEGIN {FS=""} /^>/ {next} {for (i=1;i<=NF;i++) A[$i]++; n++} END {for (i in A) print i,A[i]/n}' human.fasta
#!/bin/bash
cat $1 | grep '^>' | sed 's/_/ /' | cut -c 5- | sort | grep -v -E 'random|^17' | tr '_' ' ' | awk 'BEGIN {OFS="\t"} {print "chr", $1, $2, $2+$13, $3 "_" $4, $5 "_" $6}' | awk '{EXONS[$5]=EXONS[$5] " " $6} END {for (gene in EXONS) print gene ":" EXONS[gene]}'
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment