How to get DEPTH, as well as variant calling of transcriptome data
2
0
Entering edit mode
6.6 years ago

Hey for RNA-seq data analysis. how to get accurate depth, variant calling results.

any accurate practical guidelines for RNA-SEQ

THANKS

RNA-Seq • 1.5k views
ADD COMMENT
2
Entering edit mode
6.6 years ago

Read GATK best practice workflow for variant calling using RNA-seq. It is pretty much the same workflow as variant calling using genomic data but it includes splitting of the reads with 'N' operator in the CIGAR string into individual sub-reads followed by trimming of the overhanging portions if any. The VCF file will have information about the depth along with the variant calls. See: https://www.broadinstitute.org/gatk/guide/best-practices?bpm=RNAseq

ADD COMMENT
0
Entering edit mode

Hey please help me with this questions

How to calculate size of Human Transcriptome data?

ADD REPLY
0
Entering edit mode

What you mean by size here? I would strongly suggest you to look into Picard (http://broadinstitute.github.io/picard/) if your overall goal is to collect some metrics about RNA-seq data. If you need to call variants follow my answer above. 

ADD REPLY
0
Entering edit mode
6.6 years ago
Lynxoid ▴ 220

You could use samtools for depth, but you have to account for the fact that genome length > transcriptome length

ADD COMMENT

Login before adding your answer.

Traffic: 2545 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