Hi,
I have a list with a specific mutations localised on one gene. I have 100 BAM RNA-seq samples, mapped to the reference genome. Let's assume, that I want to calculate number of transcript reads with a specific mutation (SNV), e.g. A/T that is localised at 120 position of gene (e.g. BRCA1). How can I do that? I do not want to do variant calling, just calculate based on mapped BAM file.
Thank you for your tips.
•
written
by
misterie
1 answer
convert the cDNA or proteic position to a genomic position: e.g: Getting genomic regions by local coordinates (protein/transcript) How to convert mutation position from gene coordinate to genome coordinate?
and then use something like Coverage In Bam File - Bases And Overall Count ; Finding bases at specific sites using pysam pileup ; ...
Log in to answer this question.
-
Germline variant calling from BAM without reference
written by sbstevenlee •
-
How to filter reads only with single SNV
written by sunhyunchang •
-
Finding mutational burden in TCGA dataset
written by Jamie Watson •
-
How to select a "representative" transcript in multiple transcript variants from annotation with AN…
written by svlachavas
-
Translate amino acid mutation to the nucleotide with position in reference genome
written by misterie
-
How to calculate Mutation Frequency for each Gene using TCGA SNV data?
written by mausamkumravat •
-
Amplicon based sequencing Mutation Calling
written by morovatunc
-
Plot read counts graph on transcript only using midpoint/start position of reads or the whole read
written by qyang •
-
Is It Possible For A Gene To Have A Mutation At The Same Position In Both Cosmic And Dbsnp?
written by Jordan
-
Population Rate Calculation Programme
written by jackuser1979
