SAM file directionality

Understanding direction flags and how to extract them

• http://biobeat.wordpress.com/2013/03/14/separating-bam-files-to-forward-and-reverse-strand-files/
• http://biobeat.wordpress.com/2013/04/29/directional-rna-seq-part-1-extract-strand-information-from-sam-file/

 Paired-end specific discussion:

• http://ppotato.wordpress.com/2010/08/25/samtool-bitwise-flag-paired-reads/
• http://www.novocraft.com/wiki/tiki-index.php?page=Extracting+unmapped+reads+from+a+BAM+file+produced+by+novoalign

Useful samflag utility
http://picard.sourceforge.net/explain-flags.html

Trick : seqMINER

seqMINER is a Chip-seq density clustering and plotting software. A trick to keep clustering from previous computations and plot the same ordering on a new set of bam/bed files: Load your bam file, right click and choose view heatmap (bypassing the "cluster" button). Voila!

MISO sashimi plot (RNA-seq alternative splicing)

This is a modified python script from the MISO package for plotting alternative splicing events.

Modifications from the source:
1) Normalized junction counts (assumming that the junction reads are proportional to total read coverage)
2) Proper plotting of isoforms (lines do not extend beyond final exon/utr if one isoform is shorter than the other)