countdigi

#!/usr/bin/env Rscript --vanilla

suppressPackageStartupMessages({
  library(argparse)
  library(tidyverse)
})

# ----------------------------------------------------------------------------------------------------

parser<-ArgumentParser()

countdigi

idx     count   name    ptlp    opp     mism    mol     alt     chrom   patch   top     rsids
1       9139281 38.p12  True    False   False   genomic False   True    False   True    782,783,795
2       8315064 37.p13  False   False   False   genomic False   True    False   True    782,783,800
3       310452  38.p12  False   False   False   genomic True    False   False   True    1964,2719,3484
4       206704  37.p13  False   False   False   genomic False   False   True    True    7245,10248,10794
5       198592  37.p13  False   True    False   genomic False   True    False   True    10083,117222,131488
6       120216  37.p13  False   False   False   genomic False   False   False   True    366279,366777,372180
7       73063   38.p12  False   False   False   genomic False   False   False   False   504728,545765,656963
8       72307   38.p12  False   False   False   genomic False   False   True    True    10451,16947,131198
9       66116   37.p13  False   True    False   genomic False   False   False   True    7

countdigi

#!/usr/bin/env python3                                                                                
                                                                                                      
import gzip                                                                                           
import os                                                                                             
import sys                                                                                            
import collections                                                                                    
                                                                                                      
def load_rsmerge_records(fname):                                                                      
    fields = ["rshigh", "rslow", "build_id", "orien", "create_time", "last_updated_time", "rscurrent", "orien2current", "comment"]
    Record = collections.namedtuple("Record", fields)