Re: [Samtools-devel] Wgsim patch

[Nils H. <nil...@gm...>]

I found another spot too.

Index: wgsim.c
===================================================================
--- wgsim.c    (revision 543)
+++ wgsim.c    (working copy)
@@ -250,6 +250,7 @@
                     int n = (c[1]&mutmsk) >> 12, ins = c[1] >> 4;
                     while(n > 0) {
                         putchar("ACGTN"[ins & 0x3]);
+                        ins >>= 2;
                         n--;
                     }
                     printf("\t-\n");
@@ -266,6 +267,7 @@
                     int n = (c[1]&mutmsk) >> 12, ins = c[1] >> 4;
                     while (n > 0) {
                         putchar("ACGTN"[ins & 0x3]);
+                        ins >>= 2;
                         n--;
                     }
                     printf("\t+\n");



On 3/25/10 5:48 PM, "Nils Homer" <nil...@gm...> wrote:

> The simulated insertion was not being reported correctly if it is a homozygous
> insertion.  Below is the patch.
> 
> Sincerely,
> 
> Nils Homer
> 
>  svn diff:
> Index: wgsim.c
> ===================================================================
> --- wgsim.c    (revision 543)
> +++ wgsim.c    (working copy)
> @@ -250,6 +250,7 @@
>                      int n = (c[1]&mutmsk) >> 12, ins = c[1] >> 4;
>                      while(n > 0) {
>                          putchar("ACGTN"[ins & 0x3]);
> +                        ins >>= 2;
>                          n--;
>                      }
>                      printf("\t-\n");