Variant calling from RNASeq · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

Variant calling from RNASeq

I would like to know if it is possible to do variant calling using single end RNAseq.

rna-seq snp

• 1,416 views • link

written 5.9 years ago by ifudontmind_plzz

@Kevin, thank you. I was wondering using single end RNA sequencing, how efficient it will be handling splice junctions.

• 0 views • link

written 5.9 years ago by ifudontmind_plzz

I have no frame of reference for that. However, if your intention is actually to identify splice-junctions, then there is likely some program already produced for this that can take single-end reads. A quick search leads me to this, what seems to be a review of some sort: https://genomebiology.biomedcentral.com/articles/10.1186/s13059-015-0697-y

• 0 views • link

written 5.9 years ago by Kevin Blighe

1 answer

Yes, it is possible, but please consider the limitations: A: Inferring genotype based on RNA sequnces

Note also that Broad institute eventually removed the page to which I linked in my answer.

Kevin

• 0 views • link

written 5.9 years ago by Kevin Blighe

@Kevin, thank you. I was wondering using single end RNA sequencing, how efficient it will be handling splice junctions.

• 0 views • link

written 5.9 years ago by ifudontmind_plzz

Log in to answer this question.

More posts like this

Overlapping reads, how do we counter them while variant calling, haplotype calling, does it affect …
written 1.8 years ago by ayeraselvan •

I have Illumina paired end reads. My plan is on to do the mapping to the reference, sort them, remove the duplicates and use GATK …
Significance of "orientation" flag in FilterMutectCalls output
written 3.2 years ago by bt_cepo •

Hi! I have an annotated vcf file generated by mutect2 and flagged by FliterMutectCalls. Thing is, I have found my variant of interest but it …
variant calling using pseudoaligners
written 7.0 years ago by jmgoldstein7 •

Hello, I have recently been using a pseudoaligner (Salmon) to run RNA-Seq for differential gene expression, which has a lot of documentation and works very …
Normalizing ChIP-seq and ATAC-seq data?
written 7.5 years ago by star

I have some Single end and Paired end ChIP-seq and ATAC-seq data. I have done aligning using Bowtie2 and ignored duplicated reads using Picard, then …
SNP Calling of a single bacterial gene
written 8.3 years ago by Mirza.Jawad •

I have one sequenced bacterial gene (three copies) taken from three different countries. I want to see the SNPs and variant calling (if possible for …
How to compare the correlation between the RNAseq data from human and mouse tumor?
written 10.9 years ago by ifudontmind_plzz

Hello, I would like to compare the correlation between the gene expression pattern of mouse and human tumor type. We have generated a mouse tumor …
How to annotate gene names for genomic coordinates?
written 11.1 years ago by ifudontmind_plzz

I want to annotate gene names for a segmented file by finding overlap with it with a reference file. Would like to know how can …
gVCF from VCF and BAM
written 11.2 years ago by tejaswikoganti

Hello, Is it possible to get a gVCF from VCF file and a BAM file? I know I can get a gVCF from BAM file …
RNA-seq differential expression analysis, which aligner to choose between BWA/tophat/Bowtie?
written 11.5 years ago by EagleEye

I have RNAseq reads (single-end, 49-bp), I would like to do differential expression analysis. I know the standard pipeline is some rnaseq aligners ( e.g. …
Mapping To Different Genome Build
written 15.4 years ago by Mpiro •

<p>I have an illumina NGS dataset and would like to use bwa to align it to the human reference genome for variant calling. I know …