Skip to content
This repository has been archived by the owner on Oct 17, 2022. It is now read-only.

Commit

Permalink
methcount section modified in methpipe manual
Browse files Browse the repository at this point in the history
  • Loading branch information
amalthomas111 committed Nov 12, 2015
1 parent 439233e commit 35bd2df
Show file tree
Hide file tree
Showing 2 changed files with 6 additions and 10 deletions.
Binary file modified docs/methpipe-manual.pdf
Binary file not shown.
16 changes: 6 additions & 10 deletions docs/methpipe-manual.tex
Original file line number Diff line number Diff line change
Expand Up @@ -419,22 +419,18 @@ \subsection{Computing single-site methylation levels}
of methylation is asymmetric since the cytosines on the complementary
strand do not necessarily have the same methylation status.

The input is in mapped read format, and the reads should be sorted
according to (chrom, end, start, strand). If your reads are not
The input is in mapped read format. The reads should be sorted and duplicate reads should be removed (same steps as in section 2.2). If your reads are not
sorted, run:

\begin{verbatim}
$ LC_ALL=C sort -k 1,1 -k 3,3n -k 2,2n -k 6,6 \
-o Human_ESC.mr.sorted_end_first Human_ESC.mr
$ LC_ALL=C sort -k 1,1 -k 2,2n -k 3,3n -k 6,6 \
-o Human_ESC.mr.sorted_start Human_ESC.mr
\end{verbatim}

Since \prog{methcounts} can only take one input file, if you have
multiple you can merge them using the \op{-m} option to the
\prog{sort} program:
If duplicate reads are not removed, run:

\begin{verbatim}
$ LC_ALL=C sort -m -k 1,1 -k 3,3n -k 2,2n -k 6,6 \
-o Human_ESC.mr.sorted_end_first Human_ESC.mr.1 Human_ESC.mr.2
$ duplicate-remover -S Human_ESC_dremove_stat.txt \
-o Human_ESC.mr.dremove Human_ESC.mr.sorted_start
\end{verbatim}

\paragraph{Running \prog{methcounts}:}
Expand Down

0 comments on commit 35bd2df

Please sign in to comment.