bcl-convert and UMI's

Hi,

We're looking into using bcl-convert as an alternative for Picard wrt "UMI aware" demultiplexing.

Picard can demultiplex straight to bam with the UMI in the correct RX tag. bcl-convert appends the UMI sequence to the fastq header line when using the OverrideCycles option.

I'm looking for a way to get the UMI sequence from the fastq header into the sam RX tag (with or without alignment). I've already looked into fgbio, umi-tools, umis; but none of these tools seem to be able to handle the output format from bcl-convert.

Does anyone have a suggestion?

Thanks M

bcl-convert umi

• 2,160 views • link

updated 4.7 years ago by GenoMax

• written 4.7 years ago by matthias.desmet •

none of these tools seem to be able to handle the output format from bcl-convert

Output from bcl-convert will be fastq files. Indexes can be in separate files if needed. Are you using a specific library kit?

• 0 views • link

written 4.7 years ago by GenoMax

0 answers

No answers yet.

More posts like this

UMI workflow resulting in bams with empty reads
written 2.5 years ago by ww22runner •

Hello all, In my NGS workflow for UMI based reads, I first tried identifying and removing sequence adapters using bbmerge and cutcadapt: BBMERGE -Xmx1g -ignorejunk …
Removing umi split off as a separate fastq (RNA-seq)
written 3.6 years ago by YJ •

Hi, I used x-gen udi/umi adaptors from idt to generate my RNA-seq samples and ran my single end RNA-seq experiment. I received two fastq files …
How does bcl-convert handle UMIs in the 3' end of the insert
written 4.0 years ago by ingerslev •

How does bcl-convert handle reads with an UMI in the 3' end and adapter sequences? The bcl-convert [documentation](https://support-docs.illumina.com/SW/DRAGEN_v40/Content/SW/DRAGEN/SampleSheet.htm) is unclear on the order. If I …
How to combine UMIs split across reads 1 and 2 into a single UMI for deduplication?
written 5.2 years ago by SethG •

I am building small RNA libraries designed like this: [5' adapter]-[6N umi] --- [insert] --- [6N umi]-[3' adapter] I collect paired end data but map …
Transferring UMI from paired-end read 2 to header of read 1
written 5.7 years ago by robert.gnuegge •

I would like to deduplicate reads using UMIs. I have paired-end data where read 1 contains the mappable sequence, while read 2 only contains a …
How to cut some part of the header of reads in FASTQ file?
written 6.1 years ago by naeem40thju •

The head of reads contains long characters. Actually I wanted to use UMI-tools. If I extract UMI and add to the header ('_' separated), it …
UMI read consensus calling
written 6.5 years ago by SemiQuant

Hi I have amplicon sequence data where I included a UMI (unique molecular identifier) on my reads to allow me to correct sequencing errors. I …
Demultiplexing and deduplicating using barcodes and UMIs in the "mate" read
written 7.2 years ago by lech.kaczmarczyk •

Hi All, I'd like to demultiplex and deduplicate reads using in-line barcodes from read2 and UMIs from the read1 and read2. The format is like …
fgbio: Calling consensus reads using UMI
written 8.6 years ago by Nicolas Rosewick

Hi, I'm struggling a little bit with some data containing UMIs. It's enriched DNA-seq on specific genomic regions. To remove PCR duplicates and potential sequencing …
Demultiplexing bcl files into fastq files
written 9.0 years ago by Javad

Hello, I have bcl formatted files generated by illumina nextseq 500 machine and I want to demultiplex them into fastq files. I can not install …