About command line IGV screenshot
0
0
Entering edit mode
6.9 years ago
venu 7.0k

Dear all,

I am trying to take the screenshots of the regions of my interest loading aligned BAM file into IGV through command line. I've found very useful threads on this forum (Automated Plotting Of Seq. Reads In A Bam File). But I am interested to get a screenshot which shows bases in reads. The answers provided in the link is useful for getting screenshot 'showing mismatched bases'. I've gone through IGV webpage and found this and I've included show all bases command in batch_file.txt but no use. Isn't it possible to get the screenshot which shows bases in reads through command line?

Thank you

igv-screenshot • 3.2k views
ADD COMMENT
0
Entering edit mode

Not that I know of.

ADD REPLY
0
Entering edit mode

I tried several ways and at last I gave up on this. But is it possible to take screenshot covering the entire depth of a position through command line with little bit high clarity?

ADD REPLY
2
Entering edit mode

I also could not get IGV to display bases, but you could have a look at the much more basic samtools tview. You can at least get ASCII or HTML...

samtools tview -d H -p c1:10000 my.bam > my-region-of-interest.html
ADD REPLY
0
Entering edit mode

Thank you. But the output is somewhat obscure for me. I've looked at the region chr12:46231462 and the html file is showing

46231471     46231481      46231491     46231501     46231511

If I understand it correct, the first column is chr12:46231462, is it?

ADD REPLY
0
Entering edit mode

Yes, I think -p specifies the first column.

ADD REPLY

Login before adding your answer.

Traffic: 2133 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6