|
1 | 1 | ''' |
2 | | -============================================================================== |
| 2 | +==================================================================== |
3 | 3 | prepare_for_em - make the output from dedup compatible with EM tools |
4 | | -=============================================================================== |
5 | | -prepare_for_rsem - make output from dedup or group compatible with RSEM or |
| 4 | +==================================================================== |
| 5 | +
|
| 6 | +prepare_for_em - make output from dedup or group compatible with RSEM or |
6 | 7 | Salmon |
7 | 8 |
|
8 | | -Usage: umi_tools prepare_for_rsem [OPTIONS] [--stdin=IN_BAM] [--stdout=OUT_BAM] |
| 9 | +Usage: umi_tools prepare_for_em [OPTIONS] [--stdin=IN_BAM] [--stdout=OUT_BAM] |
9 | 10 |
|
10 | 11 | note: If --stdout is ommited, standard out is output. To |
11 | 12 | generate a valid BAM/SAM/CRAM file on standard out, please |
12 | 13 | redirect log with --log=LOGFILE or --log2stderr |
13 | 14 |
|
14 | | -
|
15 | | -
|
16 | 15 | The SAM format specification states that the mnext and mpos fields should point |
17 | 16 | to the primary alignment of a read's mate. However, not all aligners adhere to |
18 | 17 | this standard. |
|
42 | 41 |
|
43 | 42 | prepare-for-em specific options |
44 | 43 | ------------------------------- |
45 | | -
|
46 | | -==================== |
47 | | ---tags=TAG[,TAG....] |
48 | | -==================== |
| 44 | +""""""""""""""""""""""""" |
| 45 | +``--tags =TAG[,TAG....]`` |
| 46 | +""""""""""""""""""""""""" |
49 | 47 | List of SAM tags that are transfered from read1 to read2. The default |
50 | 48 | is UG and BX, which is the numeric UMI group, and the infered true UMI |
51 | 49 | respectively. |
|
54 | 52 | ''' |
55 | 53 |
|
56 | 54 | from umi_tools import Utilities as U |
| 55 | +from umi_tools import Documentation |
57 | 56 | from collections import defaultdict, Counter |
58 | 57 | import pysam |
59 | 58 | import sys |
60 | | -from umi_tools import Documentation |
61 | 59 |
|
62 | 60 | __doc__ = __doc__ + Documentation.GENERIC_DOCSTRING_SBCRAM_INPUT + Documentation.GENERIC_DOCSTRING_SBCRAM_OUTPUT |
63 | 61 |
|
64 | 62 | usage = ''' |
65 | | -prepare_for_rsem - make output from dedup or group compatible with RSEM |
| 63 | +prepare_for_em - make output from dedup or group compatible with EM tools |
66 | 64 |
|
67 | | -Usage: umi_tools prepare_for_rsem [OPTIONS] [--stdin=IN_BAM] [--stdout=OUT_BAM] |
| 65 | +Usage: umi_tools prepare_for_em [OPTIONS] [--stdin=IN_BAM] [--stdout=OUT_BAM] |
68 | 66 |
|
69 | 67 | note: If --stdout is ommited, standard out is output. To |
70 | 68 | generate a valid BAM file on standard out, please |
|
0 commit comments