Sam Lipworth
Sam Lipworth
OK great thanks for looking into it
+ 1 for @dlaehnemann's suggestion - if this could be native it would be great. @dlaehnemann your stack answer is great - anyway to add a legend to show what...
thanks those are both really helpful suggestions - Example 5 does indeed work.
@krassowski very neat, thanks
Thanks very much - I can provide a sample dataset for diagnostics if helpful.
The problem appears to happen on one side of the gene - so the script behaves as expected to the left of the gene but not the right: e.g. group_4691,orf05,group_3988,yjcD,group_4950,aac6IIc~~~aac61bcr~~~~~~aacA4,blaOXA1,cat_1,group_5079,group_3814,group_1997,group_907,group_4665,group_2500,group_6706,group_7962,merB,group_8452,group_9561,group_6290,group_2197,group_7483,group_1416,group_323,group_9396,group_8617,group_9956,group_4839,sul2,strA~~~aph3Ib,group_4838,group_5062,group_8245,group_4951,group_2686,group_7436,yghA_1,group_10754,group_262,group_1432,group_2308,yghA~~~yghA_2~~~yghA_3,tnpA_2...
Great! thanks again for another speedy and helpful response!
Had this same issue - worked fine when compiled from source
e.g. group,icd10,dcsi CV_IHD,I2(4|0|5(?!2)),1
looks like this is a wider problem - R regex doesn't do lookaheads - could add an option to turn on perl? Workaround for now is to filter out any...