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
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
Comment