Daniel:Notebook/GenomeMiner/2013-7-31

From ZhangLabWiki
Revision as of 21:52, 1 August 2013 by >Djacobse
Jump to navigation Jump to search

HL155

Back to Calendar

Started with the sam file of only mapped reads. Then printed the third column into text file.

Awk third column from SAM

awk '{print $3}' hrcp_onlymapped > test.txt

Next I used perl to count the reads using this script: