This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Simulate short-read RNA-seq data from long-read RNA-seq data

Hi,

I'm trying to use my self-developed tool to test it on long-read and short-read single-cell RNA-seq datasets. To ensure a fair comparison, I'm attempting to simulate short-length single-cell RNA-seq data from the long-read RNA-seq data. For example, for each sequence in the long-read FASTQ file, I will randomly select sequences that are 100 bp long to simulate a Smart-seq dataset. I will then use my tool for comparison. I'm wondering if there is a tool that can perform a similar task. I've noticed some people using Polyester to achieve this with transcript datasets, but I'm curious if I can directly use the FASTQ data to simulate short-read RNA-seq.

thank you!

Best, R

polyester single-cell simulation short-read long-read
• 1,726 views • link
written by rhonddaskl • rhonddaskl

1 answer

I think you could try doing this, but it wouldn't be appropriate. You'd be bringing long-read errors, which I believe are higher than short-read errors, into your "simulated" short reads. Separately, you wouldn't be creating any type of error that are hallmarks of short reads, so this would be a short-read sample only in name.

If I were to do this I'd assemble long reads first which will correct some of those errors, and simulate from the assembled material.

• 0 views • link
written by Mensur Dlakic Mensur Dlakic

Log in to answer this question.

More posts like this