# HG changeset patch # User greg # Date 1448114247 18000 # Node ID 551630d1fae3f27b82548b06e2c8a0ff01924a68 # Parent a952b6740fb9f67845324a8054e8b3c0c4f75d83 Uploaded diff -r a952b6740fb9 -r 551630d1fae3 genetrack.py --- a/genetrack.py Sat Nov 21 08:57:17 2015 -0500 +++ b/genetrack.py Sat Nov 21 08:57:27 2015 -0500 @@ -1,13 +1,8 @@ """ genetrack.py -Input: either ssccidx or gff format of reads -.ssccidx format: tab-separated chromosome (chr##), index, + reads, - reads -.gff format: standard gff, score interpreted as number of reads - -Output: Called peaks in either gff or txt format -.txt format: tab-separated chromosome, strand, start, end, read count -.gff format: standard gff, score is read count +Input: either scidx or gff format of reads +Output: Called peaks in gff format """ import optparse import csv @@ -33,7 +28,7 @@ # Make sure the reads for each chromosome are sorted by index. input_path = genetrack_util.sort_chromosome_reads_by_index(dataset_path) else: - # We're processing ssccidx data. + # We're processing scidx data. input_path = dataset_path output_name = 's%se%su%sd%sF%s_on_data_%s' % (str(options.sigma), str(options.exclusion),