X-Git-Url: http://woldlab.caltech.edu/gitweb/?a=blobdiff_plain;ds=sidebyside;f=bam_md.c;h=3ca730993cbc01065eb2c76bd66df1640808e032;hb=63f0c5147985e4b1c01942d4525c04e017dbe966;hp=8d074870d336677cee2e606ab05ee81a2b4a6f1e;hpb=d363084f0412f3bcdeb0304aeb0974c9a10c7649;p=samtools.git

diff --git a/bam_md.c b/bam_md.c
index 8d07487..3ca7309 100644
--- a/bam_md.c
+++ b/bam_md.c
@@ -132,8 +132,11 @@ int bam_fillmd(int argc, char *argv[])
 				free(ref);
 				ref = fai_fetch(fai, fp->header->target_name[b->core.tid], &len);
 				tid = b->core.tid;
+				if (ref == 0)
+					fprintf(stderr, "[bam_fillmd] fail to find sequence '%s' in the reference.\n",
+							fp->header->target_name[tid]);
 			}
-			bam_fillmd1(b, ref, is_equal);
+			if (ref) bam_fillmd1(b, ref, is_equal);
 		}
 		samwrite(fpout, b);
 	}