yanailab / celseq2

Generate the UMI count matrix from CEL-Seq2 sequencing data
https://yanailab.github.io/celseq2/
BSD 3-Clause "New" or "Revised" License
20 stars 13 forks source link

Add strand='rev' option to count-umi #1

Closed Puriney closed 7 years ago

Puriney commented 7 years ago

strand="rev" is needed when counting UMIs. This is requested by internal users.

2 possible solutions:

  1. demultiplexing. Trans-rev the R2.fastq reads before writing to demultiplexed reads files. Advantage is that the rest pipeline is untouched at all.
  2. count-umi. Add strand="rev" when iterating each read. Context is here And probably adapting these 7 lines to our pipeline will solve the issue.

2nd is favored as it will have wider scope.

Solution Refs: