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?