Mercurial > repos > yating-l > regtools_junctions_extract
diff validator.py @ 3:24fc8a8efe19 draft default tip
planemo upload commit ca69686dfafcabb815c93fd46d3c4dfe57459e39-dirty
author | yating-l |
---|---|
date | Tue, 20 Dec 2016 17:15:16 -0500 |
parents | cf258ca024ff |
children |
line wrap: on
line diff
--- a/validator.py Mon Dec 19 14:49:37 2016 -0500 +++ b/validator.py Tue Dec 20 17:15:16 2016 -0500 @@ -7,7 +7,6 @@ def checkAndFixBed(bedfile, revised_file): # Store the lines that have been removed removedLines = [] - scoreLines = [] # Remove the lines with invalid strand, create a score column to store the original scores and change scores in the original score column all to 1000 with open(revised_file, 'w') as tmp: with open(bedfile, 'r') as f: @@ -17,9 +16,8 @@ fields = line.split() strand = fields[5] score = fields[4] - if (int(fields[4]) > 1000): - scoreLines.append("line" + str(i) + ":" + line) - fields[4] = '1000' + fields[4] = '1000' + fields.append(score) if (strand == '+' or strand == '-'): tmp.write('\t'.join(map(str, fields))) tmp.write("\n") @@ -27,18 +25,15 @@ removedLines.append("line" + str(i) + ": " + line) i = i+1 - return removedLines, scoreLines + return removedLines def main(): inputfile = str(sys.argv[1]) outputfile = str(sys.argv[2]) - removed, changed = checkAndFixBed(inputfile, outputfile) + removed = checkAndFixBed(inputfile, outputfile) if (removed != []): print "\nRemoved invalid lines: \n" print "\n".join(removed) - if (changed != []): - print "\nThe following lines have scores > 1000, so they are changed to 1000:\n" - print "\n".join(changed) if __name__ == "__main__": main()