X-Git-Url: https://git.donarmstrong.com/?p=samtools.git;a=blobdiff_plain;f=bam_tview.c;h=1967b7c7c99f8e70400c4989a925a464dc4fa7e9;hp=20e6f466dfdf9ab0a2ba444e4181a20115b0632e;hb=0b184e9ada57345b3de076c688a18a913fb96551;hpb=0b3418cf166ce4a58cedf0d9a2df5ec3dd4cc5fa

diff --git a/bam_tview.c b/bam_tview.c
index 20e6f46..1967b7c 100644
--- a/bam_tview.c
+++ b/bam_tview.c
@@ -441,7 +441,8 @@ void error(const char *format, ...)
         fprintf(stderr, "\n");
         fprintf(stderr, "Usage: bamtk tview [options] <aln.bam> [ref.fasta]\n");
         fprintf(stderr, "Options:\n");
-        fprintf(stderr, "   -s STR      display only reads from this sample\n");
+        fprintf(stderr, "   -p chr:pos      go directly to this position\n");
+        fprintf(stderr, "   -s STR          display only reads from this sample\n");
         fprintf(stderr, "\n\n");
     }
     else
@@ -458,17 +459,24 @@ void error(const char *format, ...)
 int bam_tview_main(int argc, char *argv[])
 {
 	tview_t *tv;
-    char *samples=NULL;
+    char *samples=NULL, *position=NULL;
     int c;
-    while ((c = getopt(argc, argv, "s:")) >= 0) {
+    while ((c = getopt(argc, argv, "s:p:")) >= 0) {
         switch (c) {
             case 's': samples=optarg; break;
+            case 'p': position=optarg; break;
             default: error(NULL);
         }
     }
 	if (argc==optind) error(NULL);
 	tv = tv_init(argv[optind], (optind+1>=argc)? 0 : argv[optind+1], samples);
-	tv_draw_aln(tv, 0, 0);
+    if ( position )
+    {
+        int _tid = -1, _beg, _end;
+        bam_parse_region(tv->header, position, &_tid, &_beg, &_end);
+        if (_tid >= 0) { tv->curr_tid = _tid; tv->left_pos = _beg; }
+    }
+	tv_draw_aln(tv, tv->curr_tid, tv->left_pos);
 	tv_loop(tv);
 	tv_destroy(tv);
 	return 0;