Dear Grand Father, You are genius and your code is perfect, it is filling up the gaps (taking care of even-odd numbers). However, in the range plus/minus 8, for example for PIP=53-337 in the case of Contig100 if I want to have all "1"s what can I do? Lets say 45 to 345 all take 1 and before 45 all 0s. Thanks again, Pedro

.
.
.
.
          
        CLS_S3_Contig100        40                      0
        CLS_S3_Contig100        41                      0
        CLS_S3_Contig100        42                      0
        CLS_S3_Contig100        43                      0
        CLS_S3_Contig100        44                      0
        CLS_S3_Contig100        45                      0
        CLS_S3_Contig100        46                      0
        CLS_S3_Contig100        47                      0
        CLS_S3_Contig100        48                      0
        CLS_S3_Contig100        49                      0
        CLS_S3_Contig100        50                      0
        CLS_S3_Contig100        51                      0
        CLS_S3_Contig100        52                      0
CLS_S3_Contig100_st     CLS_S3_Contig100        53      10      0.3717
+  1
        CLS_S3_Contig100        54                      0
CLS_S3_Contig100_at     CLS_S3_Contig100        55      11      0.4321
+  1
        CLS_S3_Contig100        56                      0
CLS_S3_Contig100_st     CLS_S3_Contig100        57      10      0.3223
+  1
        CLS_S3_Contig100        58                      0
CLS_S3_Contig100_at     CLS_S3_Contig100        59      11      0.4055
+  1
        CLS_S3_Contig100        60                      0
CLS_S3_Contig100_st     CLS_S3_Contig100        61      11      0.4511
+  1
        CLS_S3_Contig100        62                      0
CLS_S3_Contig100_at     CLS_S3_Contig100        63      11      0.474 
+  1
        CLS_S3_Contig100        64                      0
.
.
.
.data flow...
[download]