Population Rate Calculation Programme

I have RNA seq of illumina paired reads mapped to reference genome (bam file) and VCF file from variant calling software. I would like to calculate population mutation rate (theta). Is there any good available programme or scripts already available?

biology genetics

• 3,116 views • link

updated 3.5 years ago by Ram

• written 13.8 years ago by jackuser1979

1 answer

I would not suggest calculating the population mutation rate with RNA data. WGS would be your best bet. your estimate of theta will be much lower if you are just considering mRNA.

• 0 views • link

written 13.8 years ago by Zev.Kronenberg

More posts like this

Genome assembly
written 5.1 years ago by ph19042 •

I have Illumina sequenced paired end reads of an arthropod genome. I already mapped my reads to reference genome and generated a bam file and …
Number of transcript reads with given mutation
written 6.6 years ago by misterie

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 …
Snp Contig-Synonymous & Nonsynonymous
written 13.8 years ago by jackuser1979

<p>I am working on non-model organism. I have called SNPs by read mapping illumina reads to reference genome(assembled contigs/scaffolds). I have VCF file with SNP …
Calculating Ds/Dn Ratio Through Contig Snp Annotation
written 13.8 years ago by jackuser1979

<p>I am newbie to NGS</p> <p>I have mapped my illumina reads to non-model assembled reference contigs. By variant calling I have VCF format with SNPs …
Tajima'S D Test-Ngs Data
written 13.8 years ago by jackuser1979

I have RNA seq illumina paired-end reads mapped to reference genome using bowtie and have bam file from samtools. I am interested in calculating Tajima's …
Snp-Annotation And Effect Prediction
written 13.8 years ago by jackuser1979

I have mapped RNA-seq illumina paired-end reads to reference genome contigs(which is assembled contigs/scaffolds). I have VCF file called from variant caller, which consists of …
No. Of Reads & No. Of Bases Per Locus-Calculation
written 13.8 years ago by jackuser1979

I have RNAseq illumina paired reads mapped to reference genome using bowtie and samtools. I have both bam file and mpileup file. Are there any …
Population Mutation Rate And Snp
written 13.9 years ago by jackuser1979

In our lab we have species which inbreeded over 100 generations. We have identified SNP, zygosity (heterozygous & homozyous) using illumina paired-end reads of same …
Snp Cluster Analysis
written 13.9 years ago by jackuser1979

I have illumina paired-end reads mapped to reference genome using bowtie and created mpileup using samtools and from mpileup identified SNPs using variant caller (varscan). …
Population Mutation Rate And Zygosity
written 13.9 years ago by jackuser1979

I have identified frequency of SNPs, heterozygotes and homozygotes per 100kb for my eukaryotic species (which is maintained for 100 generations by inbreeding) by read-mapping …