Fastq: how can I check if they are from DNA or RNAseq data?
2
0
Entering edit mode
2.4 years ago
svp ▴ 680

I have (gave me) Illumina fastq files, which I want to use for variant calling, and I do not know if they are DNAseq or RNAseq data. How can I check this? I do not have any report or who to ask ?

WES Fastq • 1.1k views
ADD COMMENT
4
Entering edit mode
2.4 years ago
GenoMax 147k

Align them to the reference genome and check with a genome browser to see where the reads align. With good RNAseq data the reads should mainly align to exons, with DNAseq reads should more or less uniformly align to entire reference.

ADD COMMENT
1
Entering edit mode
2.4 years ago

I guess after mapping against a genome you will see that quite a few reads map to introns. Use a non splicing-aware program, such as bowtie2, BWA or the like

A clue will be given by analyzing the level of duplication, because in RNA is much higher. Some poorly sequenced DNA however, can also have overrepresented sequences

ADD COMMENT

Login before adding your answer.

Traffic: 1414 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6