Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Frequency distribution of base quality scores from SAM file VeBeKay Bioinformatics 3 02-11-2014 12:48 AM
script to calculate A,T,G,C frequency for each position in an alingment pony2001mx Bioinformatics 8 12-02-2013 08:32 PM
How do i get the gap position from my .sam file? Petrichor Bioinformatics 1 09-23-2013 06:43 AM
Issue with Sam-Bam conversion samtools - how to remove last line of Sam file? TabeaK Bioinformatics 3 11-19-2012 10:05 AM
reads position in reference sequence(sam file) me91 Bioinformatics 4 08-27-2012 10:27 PM

Thread Tools
Old 02-02-2015, 03:50 PM   #1
Location: UK

Join Date: Jan 2014
Posts: 75
Default Frequency of A, C, T, or G at every position in SAM file

Dear all,

We would like to extract the read depth of each base type (A, C, T or G) for every position of our aligned reads (a SAM file created using Bowtie2).

Does anyone know of a tool that does this?

The closest link I found here was:

But I'm not sure that would be appropriate for use with aligned genome data.

Many thanks
M4TTN is offline   Reply With Quote
Old 02-02-2015, 11:49 PM   #2
Senior Member
Location: Cambridge, UK

Join Date: May 2010
Posts: 311

Hi- I think pysamstats can do what you need.
dariober is offline   Reply With Quote
Old 02-03-2015, 02:19 AM   #3
Location: UK

Join Date: Jan 2014
Posts: 75

Dear dariober,

We've got pysamstats up and running and it seems to be generating exactly the kind of output we require!

Many thanks for your help, Matt
M4TTN is offline   Reply With Quote

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 12:03 AM.

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