{ "metadata": { "name": "" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "heading", "level": 3, "metadata": {}, "source": [ "Create a file formatted for visualization in the R package methylKit" ] }, { "cell_type": "heading", "level": 3, "metadata": {}, "source": [ "M1" ] }, { "cell_type": "code", "collapsed": false, "input": [ "!awk '{print $1,$2,$2+1,$3,$8,($7/$8),(1-($7/$8))}' /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_FilteredFiles/methylkit_out_M1.txt " ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 3 }, { "cell_type": "code", "collapsed": false, "input": [ "#Second line of code for creating a formatted file for methylKit\n", "!tr ' ' \"\\t\" /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_FilteredFiles/methylKit/methylkit_out_M1final.txt" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 4 }, { "cell_type": "code", "collapsed": false, "input": [ "!head /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_FilteredFiles/methylKit/methylkit_out_M1final.txt" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "C12768\t103\t104\t+\t6\t0.166667\t0.833333\r\n", "C12806\t76\t77\t+\t5\t0.2\t0.8\r\n", "C12806\t78\t79\t+\t5\t0.2\t0.8\r\n", "C12806\t142\t143\t+\t8\t0.375\t0.625\r\n", "C12924\t30\t31\t+\t5\t0\t1\r\n", "C12924\t38\t39\t+\t5\t0\t1\r\n", "C12924\t52\t53\t+\t6\t0\t1\r\n", "C12924\t60\t61\t+\t6\t0\t1\r\n", "C12924\t127\t128\t+\t6\t0\t1\r\n", "C12924\t136\t137\t+\t6\t0\t1\r\n" ] } ], "prompt_number": 5 }, { "cell_type": "heading", "level": 3, "metadata": {}, "source": [ "M3" ] }, { "cell_type": "code", "collapsed": false, "input": [ "!awk '{print $1,$2,$2+1,$3,$8,($7/$8),(1-($7/$8))}' /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_FilteredFiles/methylkit_out_M3.txt " ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 6 }, { "cell_type": "code", "collapsed": false, "input": [ "#Second line of code for creating a formatted file for methylKit\n", "!tr ' ' \"\\t\" /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_FilteredFiles/methylKit/methylkit_out_M3final.txt" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 7 }, { "cell_type": "code", "collapsed": false, "input": [ "!head /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_FilteredFiles/methylKit/methylkit_out_M3final.txt" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "C13128\t87\t88\t+\t5\t0\t1\r\n", "C13208\t83\t84\t+\t7\t0.428571\t0.571429\r\n", "C13220\t164\t165\t+\t6\t0\t1\r\n", "C13442\t65\t66\t+\t5\t0\t1\r\n", "C13442\t79\t80\t+\t5\t0\t1\r\n", "C13612\t89\t90\t+\t5\t0\t1\r\n", "C13992\t54\t55\t+\t5\t0\t1\r\n", "C14180\t158\t159\t+\t5\t0\t1\r\n", "C14220\t96\t97\t+\t12\t0\t1\r\n", "C14220\t143\t144\t+\t7\t0\t1\r\n" ] } ], "prompt_number": 8 }, { "cell_type": "heading", "level": 3, "metadata": {}, "source": [ "T1D3" ] }, { "cell_type": "code", "collapsed": false, "input": [ "!awk '{print $1,$2,$2+1,$3,$8,($7/$8),(1-($7/$8))}' /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_FilteredFiles/methylkit_out_T1D3.txt " ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 9 }, { "cell_type": "code", "collapsed": false, "input": [ "#Second line of code for creating a formatted file for methylKit\n", "!tr ' ' \"\\t\" /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_FilteredFiles/methylKit/methylkit_out_T1D3final.txt" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 10 }, { "cell_type": "code", "collapsed": false, "input": [ "!head /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_FilteredFiles/methylKit/methylkit_out_T1D3final.txt" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "C12960\t123\t124\t+\t8\t0\t1\r\n", "C13766\t145\t146\t+\t6\t0\t1\r\n", "C14468\t214\t215\t+\t6\t0.166667\t0.833333\r\n", "C14468\t216\t217\t+\t6\t0.166667\t0.833333\r\n", "C14796\t69\t70\t+\t5\t1\t0\r\n", "C14952\t148\t149\t+\t5\t0\t1\r\n", "C14966\t42\t43\t+\t5\t0.2\t0.8\r\n", "C14966\t54\t55\t+\t5\t0.2\t0.8\r\n", "C15814\t226\t227\t+\t5\t0\t1\r\n", "C15834\t91\t92\t+\t5\t0\t1\r\n" ] } ], "prompt_number": 11 }, { "cell_type": "heading", "level": 3, "metadata": {}, "source": [ "T1D5" ] }, { "cell_type": "code", "collapsed": false, "input": [ "!awk '{print $1,$2,$2+1,$3,$8,($7/$8),(1-($7/$8))}' /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_FilteredFiles/methylkit_out_T1D5.txt " ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 12 }, { "cell_type": "code", "collapsed": false, "input": [ "#Second line of code for creating a formatted file for methylKit\n", "!tr ' ' \"\\t\" /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_FilteredFiles/methylKit/methylkit_out_T1D5final.txt" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 13 }, { "cell_type": "code", "collapsed": false, "input": [ "!head /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_FilteredFiles/methylKit/methylkit_out_T1D5final.txt" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "C12960\t67\t68\t+\t5\t0\t1\r\n", "C12960\t123\t124\t+\t8\t0\t1\r\n", "C13766\t126\t127\t+\t9\t0\t1\r\n", "C13766\t145\t146\t+\t10\t0\t1\r\n", "C14220\t96\t97\t+\t5\t0\t1\r\n", "C14220\t143\t144\t+\t5\t0\t1\r\n", "C14796\t69\t70\t+\t8\t0.75\t0.25\r\n", "C14886\t55\t56\t+\t5\t0\t1\r\n", "C14944\t199\t200\t+\t7\t0\t1\r\n", "C14944\t214\t215\t+\t7\t0\t1\r\n" ] } ], "prompt_number": 14 }, { "cell_type": "heading", "level": 3, "metadata": {}, "source": [ "T3D3" ] }, { "cell_type": "code", "collapsed": false, "input": [ "!awk '{print $1,$2,$2+1,$3,$8,($7/$8),(1-($7/$8))}' /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_FilteredFiles/methylkit_out_T3D3.txt " ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 15 }, { "cell_type": "code", "collapsed": false, "input": [ "#Second line of code for creating a formatted file for methylKit\n", "!tr ' ' \"\\t\" /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_FilteredFiles/methylKit/methylkit_out_T3D3final.txt" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 16 }, { "cell_type": "code", "collapsed": false, "input": [ "!head /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_FilteredFiles/methylKit/methylkit_out_T3D3final.txt" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "C12960\t123\t124\t+\t7\t0\t1\r\n", "C13546\t41\t42\t+\t7\t0\t1\r\n", "C13546\t84\t85\t+\t7\t0\t1\r\n", "C13992\t120\t121\t+\t6\t0\t1\r\n", "C14180\t114\t115\t+\t5\t0\t1\r\n", "C14180\t125\t126\t+\t5\t0\t1\r\n", "C14180\t130\t131\t+\t5\t0\t1\r\n", "C14180\t158\t159\t+\t5\t0\t1\r\n", "C14220\t96\t97\t+\t7\t0\t1\r\n", "C14220\t143\t144\t+\t12\t0\t1\r\n" ] } ], "prompt_number": 17 }, { "cell_type": "heading", "level": 3, "metadata": {}, "source": [ "T3D5" ] }, { "cell_type": "code", "collapsed": false, "input": [ "!awk '{print $1,$2,$2+1,$3,$8,($7/$8),(1-($7/$8))}' /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_FilteredFiles/methylkit_out_T3D5.txt " ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 18 }, { "cell_type": "code", "collapsed": false, "input": [ "#Second line of code for creating a formatted file for methylKit\n", "!tr ' ' \"\\t\" /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_FilteredFiles/methylKit/methylkit_out_T3D5final.txt" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 19 }, { "cell_type": "code", "collapsed": false, "input": [ "!head /Volumes/web/Mollusk/bs_larvae_exp/AdapterTrimmed_FilteredFiles/methylKit/methylkit_out_T3D5final.txt" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "C13220\t164\t165\t+\t5\t0\t1\r\n", "C13490\t83\t84\t+\t5\t0\t1\r\n", "C13546\t84\t85\t+\t8\t0\t1\r\n", "C14220\t96\t97\t+\t5\t0\t1\r\n", "C14220\t143\t144\t+\t5\t0\t1\r\n", "C14450\t60\t61\t+\t6\t0\t1\r\n", "C14450\t94\t95\t+\t5\t0\t1\r\n", "C14450\t106\t107\t+\t6\t0\t1\r\n", "C14450\t118\t119\t+\t5\t0\t1\r\n", "C14450\t187\t188\t+\t6\t0.166667\t0.833333\r\n" ] } ], "prompt_number": 20 }, { "cell_type": "code", "collapsed": false, "input": [ "#NOTEBOOK NAME\n", "nb='Methylkit_Formatting_AdapterTrimmed'\n", "\n", "#NOTEBOOK LOCATION\n", "nbloc='/Users/claireolson/Dropbox/Lab\\ \\(1\\)/IPython_nbs/Claire/'\n", "\n", "\n", "\n", "#TIMESTAMP\n", "date=!date +%m%d_%H%M\n", "#LOCATION OF POST DIRECTORY\n", "post='/Volumes/web/Mollusk/iPython_Posts/'" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 21 }, { "cell_type": "code", "collapsed": false, "input": [ "cd {post}" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "/Volumes/web/Mollusk/iPython_Posts\n" ] } ], "prompt_number": 22 }, { "cell_type": "code", "collapsed": false, "input": [ "!ipython nbconvert {nbloc}{nb} --to html --template full --output {nb}{date}" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "[NbConvertApp] Using existing profile dir: u'/Users/claireolson/.ipython/profile_default'\r\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "[NbConvertApp] Converting notebook /Users/claireolson/Dropbox/Lab (1)/IPython_nbs/Claire/Methylkit_Formatting_AdapterTrimmed.ipynb to html\r\n", "[NbConvertApp] Support files will be in Methylkit_Formatting_AdapterTrimmed[0408_1324]_files/\r\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "[NbConvertApp] Loaded template html_full.tpl\r\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "[NbConvertApp] Writing 208694 bytes to Methylkit_Formatting_AdapterTrimmed[0408_1324].html\r\n" ] } ], "prompt_number": 23 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] } ], "metadata": {} } ] }