Subsample.txt: Difference between revisions
Jump to navigation
Jump to search
>Noi mNo edit summary |
>Noi mNo edit summary |
||
Line 1: | Line 1: | ||
'''[[Noi/NOTES/2014-7-20| Back]]''' | |||
import sys, random | import sys, random | ||
import HTSeq | import HTSeq |
Latest revision as of 22:41, 20 July 2014
import sys, random import HTSeq fraction = float( sys.argv[1] ) in1 = iter( HTSeq.FastqReader( sys.argv[2] ) ) in2 = iter( HTSeq.FastqReader( sys.argv[3] ) ) out1 = open( sys.argv[4], "w" ) out2 = open( sys.argv[5], "w" ) while True: read1 = next( in1 ) read2 = next( in2 ) if random.random() < fraction: read1.write_to_fastq_file( out1 ) read2.write_to_fastq_file( out2 ) out1.close() out2.close()