20 millions reads per sample for SNP

Question: 20 millions reads per sample for SNP

3.2 years ago by

archie • 100

India

archie • 100 wrote:

Dear all,

I have RNASEQ sequencing data (tumour) in which each sample has approximately 20 million reads. My target is to identify SNPs and predict the isoform expression variability using RNAseq data. Before doing all analysis , I wanted to check coverage of sequence data. Just to check coverage , I used formaula N*L/G where N no of reads per sample, L read length (150 as paired end data 75+75), and human genome size. For sample, with 26130832 reads count, I got coverage 1.3. Does this 1.3 indicates sequence coverage is 1X ? Am I right ? My target is to identity SNPs. Can I proceed with this much amount of data per sample?

Thank you in advance

snp rna-seq • 1.1k views

ADD COMMENT • link •

modified 3.2 years ago by swbarnes2 ♦ 9.6k • written 3.2 years ago by archie • 100

3.2 years ago by

swbarnes2 ♦ 9.6k

United States

swbarnes2 ♦ 9.6k wrote:

For what it's worth, in my lab, we shoot for 30M reads just to assess expression differences; you'd want more for variants. For variant calling, the average coverage is not helpful. You need a finer breakdown...you need to know what percentage of the transcriptome is adequately covered (say, 25x). You will need something like BEDTools to assess that. My guess is, you might have enough depth to examine the most highly expressed genes, but not much else.

ADD COMMENT • link written 3.2 years ago by swbarnes2 ♦ 9.6k

3.2 years ago by

Titus • 910

Titus • 910 wrote:

Hi ,

You should use G = EXOM , you count non coding part in your formula :)

Best

ADD COMMENT • link modified 3.2 years ago • written 3.2 years ago by Titus • 910

Dear Titus,

Yes. You are right. G=EXOM part will represent the actual targeted genome size . Thanks :)

Based on this, I think my data is good enough to proceed for analysis.

ADD REPLY • link written 3.2 years ago by archie • 100

RNA-seq has a very uneven coverage pattern. Some highly abundant genes will have high coverage (>1000 reads) while you'll also have genes with some expression (~40-80 reads), low abundant genes with just a few reads (<5) and many genes without expression in your tissue of interest (0 reads). As such, simply taking an average of the read count is pointless and meaningless.

I'm not saying that you can't, but RNA-seq is far from optimal for variant calling.

ADD REPLY • link written 3.2 years ago by WouterDeCoster ♦ 45k

Please log in to add an answer.

Similar posts • Search »