How to analyze DNA data with UMI

Hi all!

I'm facing with analyze DNA exomes sequencing data with UMI. From the demultiplexing, I have three fastq files: one for reads R1, one for the UMI tags, and one for the R2. Now, I have to align data to the reference genome and call variants. Have you some suggestions for me?

Does anyone know some "best practices" for that?

Thanks
Denise

next-gen-sequencing umi

• 2,133 views • link

updated 3.2 years ago by Ram

• written 6.4 years ago by deniselavezzari •

Have you looked at documentation of standard software such as UMI-tools?

Also, see How to take in consideration UMI in WES pipeline?. By the way, Ian (the guy who answers this question) is the author of UMI-tools.

• 0 views • link

written 6.4 years ago by ATpoint

0 answers

No answers yet.

More posts like this

Bacterial contamination in human DNA sample
written 4.1 years ago by priya.bmg •

Hello I need some suggestions. We have been facing problems with sequencing old DNA samples. Poor read length and yield. We suspect bacterial contamination in …
Forward and reverse single cell data analysis
written 5.5 years ago by elb

Hi, I have to analyze published scRNA-Seq data. Data are Chromium 10X on NovaSeq6000 and: for each sample I have forward (R1) in lane1 and …
Best way to deal with chimeric reads?
written 5.7 years ago by deniselavezzari •

Hi everyone! I have some DNA data (custom panel on human gene) that have a lot of chimeras reads, due to the library preparation procotol. …
Paired end data with R1 demultiplexed but not R2
written 6.6 years ago by singlegenome •

Hello! I just got some paired-end sequencing data that has R1 demultiplexed and in separate sample-specific folders, but R2 all in a single big file. …
fastq files from pair-end sequencing, for 4c-seq analysis, How to analyze ?
written 6.8 years ago by htchd •

Hi, I have fastq files from pair-end sequencing, for 4c-seq/Chip-seq analysis. How to analyze it? The data contains R1 and R2 fastq files. Should I …
Trouble preparing paired-end (MiSeq 2x250), double-tagged and dual-indexed reads from metabarcoding…
written 8.7 years ago by jack1120 •

I have paired-end reads (file_name_R1.fastq & file_name_R2.fastq) from a metabarcoding study that were originally amplified with forward and reverse primers that each had tags (6-8 …
Trim read-specific adapter in paired-end reads
written 8.8 years ago by Nicolas Rosewick

Hi, I've a specific enriched DNA-seq library to analyze ( 2x76 bp sequenced on a NextSeq500). The library is defined as : R1 R2 ==============>-----------------<===========#####@@@@@ …
How to align non overlapping fastq files
written 9.3 years ago by Rongen •

Hi, We have sequenced 250 bp read length library at 100 X 2 cycles on Miseq by mistake, because of which the R1 and R2 …
MCF-7 Cell Line DNA Sequence Data
written 9.7 years ago by bioinformatics.cancer

Hi, We are planning to analyze some cancer cell lines for variants such as mutations, gene fusions and structural variants. I am looking for sequencing …
Reversing R1.fastq and R2.fastq
written 9.8 years ago by L. A. Liggett

Just as a check for the substitutions I was observing, I reversed R1.fastq and R2.fastq being the paired end reads from as sequencing experiment. Let's …