This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Extracting UMI sequences from paired-end reads

Hello,

I have a paired end fastq file and my experiment is designed in a way that each PAIRED READ has ONE barcode in it. This barcode might be on the forward read or on the reverse read (not both of them) and the barcode has a specific sequence before and after it that helps to identify the barcode. So, some of the forward reads have the barcode and some of them do not. This is also true for reverse reads.

The sequence of barcode is located at the beginning of reads (of course after pre processing trimming) and it is like this: GTC NNN NNN G

Does any one know a reliable tool for extraction of UMI sequences in this experimental design and quantifying the number of unique UMIs aligned to each gene?

Thanks in advance

rna-seq next-gen
• 7,564 views • link
written by Javad Javad

seems like I was a minute late if you deem fit please retire my answer. Thanks

• 0 views • link
written by ivivek_ngs ivivek_ngs

1 answer

I reckon these are scRNA-Seq data, you should probably also mention the nature of the sequencing chemistry if its 10x Chromium, Fluidigm etc since some of them do have such issues.

However coming to your answer, @ian Sudbury and Tom Smith at CGAT Oxford has provided a pretty handy python script to extract UMI information from both SE and PE data along with an optional way of extracting cell barcodes. Check below

Extract UMI from Fastq

• 0 views • link
written by ivivek_ngs ivivek_ngs

Thanks for your answer. I have already tried their tool but the problem is when it comes to paired end reads, UMI-tools needs barcodes on both reverse and forward reads which is not the case in my experiment. Unfortunately, Their tool is not appropriate for this type of experiment design.

• 0 views • link
written by Javad Javad

Thanks. The post is updated.

• 0 views • link
written by Javad Javad

Log in to answer this question.

More posts like this