View Single Post
Old 01-13-2016, 09: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