Allele specific expression (ASE) using unfiltered reads from RNA-seq data?
Hi, all!
I'm trying to measure the allele specific expression (ASE) levels of selected genes from RNA-seq data.
My question is if it's OK to use unfiltered read counts to calculate the ASE.
As I understand, GATK does various statistical(?) tests, which I can't understand, to filter out bad reads and determine genotypes.
Thus, to calculate the ASE, I think something similar steps should be included.
I found number of tools for ASE, but I'm curious how bad the results would be if I use the unfiltered reads.
Sorry for such a vague question. I just want to hear what people think.
Thank you!
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
RNA-seq scaling the data for heatmap
written by maplewj •
-
Missing SNPs on an array using PLINK
written by Vanna Hay •
-
Panther /GO analysis
written by omicsnstuff •
-
Normalizing NGS samples expected to have different total reads
written by mbk0asis
-
How can I access the consensus sequences of all human ERVs from Dfam?
written by mbk0asis
-
Feasible Annotation of 4 million unique SNPs
written by serpalma.v •
-
Phylogenetic diversity and Phylogenetic tree
written by ashakoor •
-
How to build a custom reference genome (from a BED file)
written by umn_bist •
-
Simple shell command for extracting indexed reads from fastq?
written by mbk0asis
-
Filtering ASE variants with those in high LD using plink
written by drautuna •
Can you edit the message commenting on what you mean with filtered reads in this context? Trimmed adapters, deduplicated, quality filtering, all of them?