Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Newbie Question: How to count genome region coverage and CpG sites coverage? illuminaGA Bioinformatics 0 09-25-2015 08:25 AM
In the SAM format, how the POS field affected by insertion and deletion blaboon Bioinformatics 7 01-09-2015 05:23 PM
Getting insertion counts from DNA sequence data anjulka Bioinformatics 0 02-27-2013 02:14 PM
Per base sequence coverage from sam/bam file? ewilbanks Bioinformatics 7 06-06-2012 02:03 PM
PubMed: Tracking insertion mutants within libraries by deep sequencing and a genome-w Newsbot! Literature Watch 0 10-07-2009 03:00 AM

Thread Tools
Old 01-13-2016, 10:59 PM   #1
Junior Member
Location: Italy

Join Date: Jan 2016
Posts: 2
Default Sequence Coverage and Genome Insertion lenght from SAM

Hi to everyone,

I'm quite new of this science, and I'm looking for solve a problem.

I've a SAM file generated by BWA and must to calculate two different things:

1) Genome insertion lenght ( for plotting mean and standard dev )
2) Sequence Coverage

I've calculate the insertion lenght with this operation abs(col[3]-col[7]) made for each row.

The main problem is that i can understand how to calculate sequence coverage from SAM file, i must implement the algorithm in python or perl, but the prolem is not "the language" is wath operation i need to do?

Anyone have some ideas?

thanks in advance to everyone
sircamp is offline   Reply With Quote
Old 01-14-2016, 04:20 AM   #2
Senior Member
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,015

genomecov ( or coveragebed ( tools from Bedtools package can do this operation. You can find graphical representation of what needs to be done (since you want to write your own code) at those links.

Last edited by GenoMax; 01-14-2016 at 04:22 AM.
GenoMax is online now   Reply With Quote

coverage, insert length, sam, sam file

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 06:18 AM.

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