Question: Finding SNV showing allele specific bias
gravatar for ishmahe16
4 months ago by
ishmahe160 wrote:

I have a fastq.gz file. I want to know the SNV's in this file showing an allele specific bias ( Bias as SNV allele frequency >0.7) .How can I do that?

rna-seq gene • 257 views
ADD COMMENTlink written 4 months ago by ishmahe160

You're going to have better luck getting useful answers if you a) show evidence that you've tried to research this problem yourself b) tell us what you've tried, and c) provide more details

ADD REPLYlink written 4 months ago by Chris Miller20k

So i had 6 files 3 control and 3 primed files and had to find the most differentially expressed gene. Usong Galaxy I did FastQc then trimmed the seauences. Once that was done I did RNA STAR to aloign them followed by htseq Count . Once I had the bam files I did deseq2 and got the genes that I am looking for.

For the SNV thing I tried converting the bam file to bcf using samtools and then to vcf file . Once I had my vcf file I had done grep -c AF1 but I am stuck at that point.

ADD REPLYlink written 4 months ago by ishmahe160
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1776 users visited in the last hour