1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
|
.. _collapse:
***************
Collapse fastq(.gz) files
***************
**Definition**
Normally quality values are lost in small RNA-seq pipelines due to collapsing after adapter recognition. This option allow to collapse reads after adapter removal with ``cutadapt`` or any other tool. This way the mapping can use quality values, allowing to map using ``bwa`` for instance, or any other alignment tool that doesn't support FASTA files.
**Methods**
The new quality values are the average of each of the sequence collapse.
**Example**
::
seqcluster collapse -f sample_trimmed.fastq -o collapse
* ``-f`` is the fastq(.gz) file
* ``-o`` the folder where the outout will be created. A new FASTQ file, where the name stand for::
@seq_[0-9]_x[0-9]
The number right after ``_x`` means the abundance of this sequence in the sample
|
