How to use samtools tview goto with long reference sequence name

Question

[quokka]

I'm sorry if the answer to this is blatantly obvious...

I'm trying to jump to a sequence region in samtools tview with the (g)oto command. The problem is I can't enter the full command due to size restrictions of the input section.

I want to enter:

Goto: d62ff96e-317e-46d3-92b5-f5a4b329694f:3598

but I can only enter:

Goto: d62ff96e-317e-46d3-92b5-f5a4b329694f:359

Any tips?

Answer 1

[Pierre Lindenbaum]

edit the file https://github.com/samtools/samtools/blob/50a0bafe616a822df2f98dc4635840e430cbea45/bam_tview.h

change https://github.com/samtools/samtools/blob/50a0bafe616a822df2f98dc4635840e430cbea45/bam_tview.h#L80

#define TV_MAX_GOTO 40

to a larger value like

#define TV_MAX_GOTO 50

and recompile samtools

[quokka]

Ah. Excellent. Not so obvious then. Thanks!

[Medhat]

Hi, If this was the answer to your question please accept it as a right answer.

Answer 2

[dariober]

A shameless plug... If you use samtools tview you may be interested in ASCIIGenome

Answer 3

[John Marshall]

Sadly this text entry box is too narrow for unfortunate contigs named via UUIDs. See PR samtools/samtools#968.

In the meantime, there is a workaround to jump to the position you're after using a vanilla samtools tview build: the goto command accepts =1234 to jump to a position on the current chromosome. So you can jump to where you want in two hops:

1. Goto d62ff96e-317e-46d3-92b5-f5a4b329694f which takes you to the start of that contig;
2. Goto =3598 to jump to the desired position.