for key in genes: thisGeneName = key if len(genes[thisGeneName]) > totalReadLength: complete = makeRandomTestData.outputFastaSeq(filename + 'E1R100G100.fa', thisGeneName, genes[thisGeneName]) levelKey = random.randint(0, 2) if levelKey == 0: level = 10 elif levelKey == 1: level = 100 elif levelKey == 2: level = 1000 else: level = 0 complete = makeRandomTestData.makeRandomlyPlacedPairedEndReads(genes[thisGeneName], thisGeneName, readLength, insertSize, level, filename + 'E1R100G100RandomReadsPair1.fa', filename + 'E1R100G100RandomReadsPair2.fa', filename + 'E1R100G100RandomReads.txt', 1, 0) count += 1 genes2 = makeRandomTestData.inputFastaSeq(infileName) for i in range(0, r90): thisGeneName = random.choice(genes2.keys()) complete = makeRandomTestData.outputFastaSeq(filename + 'E1R90G100.fa', thisGeneName, genes2[thisGeneName]) genes2.pop(thisGeneName) genes2 = makeRandomTestData.inputFastaSeq(infileName) for i in range(0, r80): thisGeneName = random.choice(genes2.keys()) complete = makeRandomTestData.outputFastaSeq(filename + 'E1R80G100.fa', thisGeneName, genes2[thisGeneName]) genes2.pop(thisGeneName)