Commit 4f044389 authored by Jeroen F.J. Laros's avatar Jeroen F.J. Laros

More cleaning.

parent 69ad5347
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