예제 #1
0
def main():
    argn = len(sys.argv)
    if argn != 2:
        usage()
        sys.exit()
    seqfile = sys.argv[1]
    allSeqs = parseSequences(seqfile)
    for i in range(len(allSeqs)):
        print ">seq" + str(i), str(i)
        print allSeqs[i]
예제 #2
0
def main():
    argn = len(sys.argv)
    if argn != 2:
        usage()
        sys.exit()
    seqfile = sys.argv[1]
    allSeqs = parseSequences(seqfile)
    for i in range(len(allSeqs)):
        print ">seq"+str(i), str(i)
        print allSeqs[i]
예제 #3
0
def main():
    argn = len(sys.argv)
    if argn != 2:
        usage()
        sys.exit()
    seqfile = sys.argv[1]
    allSeqs = parseSequences(seqfile)
    print "Input Sequence File:", seqfile
    for seq in allSeqs:
        print " ", seq
    Param = Parameters()
    Stats = Statistics()
    computeSumPairsScore(allSeqs, Param, Stats)
    Stats.dumpResult()
예제 #4
0
def main():
    argn = len(sys.argv)
    if argn != 2:
        usage()
        sys.exit()
    seqfile = sys.argv[1]
    allSeqs = parseSequences(seqfile)
    print "Input Sequence File:", seqfile
    for seq in allSeqs:
        print " ", seq
    Param = Parameters()
    Stats = Statistics()
    computeSumPairsScore(allSeqs, Param, Stats)
    Stats.dumpResult()
예제 #5
0
def main():
    argn = len(sys.argv)
    if argn != 2:
        usage()
        sys.exit()
    seqfile = sys.argv[1]
    allSeqs = parseSequences(seqfile)
    numSeqs = len(allSeqs)
    allModelSeqs = []
    allDataSeqs = []
    # odd-indexed sequence is model sequence
    for i in range(0, numSeqs,2):
        allModelSeqs.append(allSeqs[i])
    # even-indexed sequence is data sequence
    for i in range(1, numSeqs,2):
        allDataSeqs.append(allSeqs[i])
    allCOSeqs = convertPair2CO (allModelSeqs, allDataSeqs)
    for seq in allCOSeqs:
        print seq