Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
counting total bases in a fasta file morning latte Bioinformatics 6 09-23-2016 05:58 AM
Counting number of reads mapped to an allele jyu429 General 5 10-18-2014 10:43 AM
total number of counts with HTseq oliviera Bioinformatics 17 07-26-2013 07:33 AM
Counting number of short-read cover positions d_t_nguyen Bioinformatics 4 06-24-2013 12:55 AM
Question on counting number of reads per gene gen2prot Bioinformatics 3 06-25-2010 10:45 AM

Thread Tools
Old 08-11-2015, 05:47 AM   #1
Location: Europe

Join Date: Feb 2013
Posts: 58
Default Counting the total number of SNPs in a transcript


I have RNAseq data that has been aligned to reference transcriptome and the aligned bam file is sorted and indexed.

Using that bam file is it possible to estimate the total number SNPs in every read that has been aligned to reference transcriptome. I mean, for an example, if position 2 has base "A" in reference transcript and 20 reads aligned to that region region has "G" in mutated_transcript, then the total number of A_to_G for that position is 20. Likewise I would like to calculate total number of base "A"s is mutated to "G"s (and vice versa) in all reads at every possible position and sum it up so that I get count for total number of base A mutated to G for the entire transcript.

Kindly guide me

Thanks in advance
dena.dinesh is offline   Reply With Quote

bam file, count, mutation, rnaseq data, snps

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 09:38 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO