WHERE <filter_expression> ''')
parser.add_argument('-s', dest='startidx',
default=0,
help='Starting base index of DNA sequences that are written to file, used to miss out cutsite if desired.')
parser.add_argument('-f', dest='format',
default='fasta',
help='Format of file written to output.')
print sys.argv
args = parser.parse_args()
# Write records to output
db = Reads_db(args.input, recbyname=True)
clusters_list = db.get_cluster_by_size()
fastafile_handle = db.write_reads(args.output, format=args.format,
filter_expression=args.filter_expression,
startidx=args.startidx)
if args.query_expression:
fastafile_handle = db.write_reads(args.pattern, args.output,
use_type_column=args.typeflag, format='fasta')
