Clinical-Genomics / demultiplexing

To keep scripts associated with execution of the Illumina demultiplexing pipeline
5 stars 0 forks source link

3.2 iSeq Sample sheet in the correct format #90

Closed emiliaol closed 4 years ago

emiliaol commented 5 years ago

As a bioinfo prod member i want the iSeq sample sheet to be in the correct format so that I don't have to do a manual intervention.

Problem: The current demux sheet fetch -a iseq -p --longest <FC ID> returns index2 in the incorrect format for direct use in demultiplexing. It needs to be reverse complement from how it is today. The reverse complement has to be done on the padded sequence.

Suggest solution:

  1. In the existing code add a reverse complement step after padding of index2 (only).
  2. update document 1161 demultiplexing

Test case on thalamus: /home/hiseq.clinical/iSeq/runs/20190528_FS10000534_7_BPC29611-3324 FC ID = BPC29611-3324 ACC5726A1 current index2 = ACCTTAGTGT, should be = ACACTAAGGT

DoD: