This is a test version of Biostars. For the public version, visit https://www.biostars.org.
igvtools count - how do I get it to show insertion and deletion count

Hello,

I am running the command:

igvtools count -w 1 --bases --query chr1:1326-1338 --minMapQuality 30 input.bam output.wig ref.fa

enter image description here

I get the correct number of A,C,G,T,N base counts as seen in IGV. However, I also need the number of insertions and deletions. They are shown in the GUI version, but I am unsure after looking at the documentation for igvtools count how to have it return these two values. I'm sure it's something simple I'm missing.

Any help you can provide would be appreciated.

Thanks!

igvtools
• 6,096 views • link
updated by Titus Titus • written by kwgsmith • kwgsmith

2 answers

I don't know about igvtools but have a look at pysamstats. See also these threads Calculate The Frequency Of Nucleotides At Each Position In An Mpileup File and Nucletotide distribution, at each position, in a .sam/.bam file ?

• 0 views • link
written by dariober dariober

I've tried pysamstats and it gives the incorrect (according to IGV) number of insertions and deletions. Oddly though it gives the same insertions and deletions that my perl script, which I wrote after I tried pileup2indel and pysamstats and before I ran across igvtools, calculates.

I tried to talk to the creator of pysamstats, but unfortunately I couldn't share my data with him or come up with something similar enough that he could help.

I've also looked at mpileup, but it will not do everything I want it to do, or at least that is my understanding after talking with someone much more knowledgeable in my lab.

Thanks for the reply though!

• 0 views • link
written by kwgsmith • kwgsmith

Hi , I m using igvtools count and if i had this options --strands read i have INDEL and the information number of count in foward and reverse ! It looks like the information you want to :)

And just a few thinks when i use the command i got all the information right (number of A,C,G,T,DEL,INS) except for INS foward information witch is always equal to 0 . Is it a bug ?

Best regards

Tristan

• 0 views • link
written by Titus Titus

Log in to answer this question.

More posts like this