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

Commit 35bd2df

Browse files
committed
methcount section modified in methpipe manual
1 parent 439233e commit 35bd2df

File tree

2 files changed

+6
-10
lines changed

2 files changed

+6
-10
lines changed

docs/methpipe-manual.pdf

6 Bytes
Binary file not shown.

docs/methpipe-manual.tex

Lines changed: 6 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -419,22 +419,18 @@ \subsection{Computing single-site methylation levels}
419419
of methylation is asymmetric since the cytosines on the complementary
420420
strand do not necessarily have the same methylation status.
421421

422-
The input is in mapped read format, and the reads should be sorted
423-
according to (chrom, end, start, strand). If your reads are not
422+
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
424423
sorted, run:
425424

426425
\begin{verbatim}
427-
$ LC_ALL=C sort -k 1,1 -k 3,3n -k 2,2n -k 6,6 \
428-
-o Human_ESC.mr.sorted_end_first Human_ESC.mr
426+
$ LC_ALL=C sort -k 1,1 -k 2,2n -k 3,3n -k 6,6 \
427+
-o Human_ESC.mr.sorted_start Human_ESC.mr
429428
\end{verbatim}
430-
431-
Since \prog{methcounts} can only take one input file, if you have
432-
multiple you can merge them using the \op{-m} option to the
433-
\prog{sort} program:
429+
If duplicate reads are not removed, run:
434430

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

440436
\paragraph{Running \prog{methcounts}:}

0 commit comments

Comments
 (0)