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
0
Entering edit mode

Not that I know of.

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?

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

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?

0
Entering edit mode

Yes, I think -p specifies the first column.