def run(self, args): quality = gem._prepare_quality_parameter(args.quality) raw = False if args.input is not None: map_file = gem.files.open(args.input, quality=quality) else: map_file = gem.files.open(sys.stdin, quality=quality) raw = True cons = gem.extended_splice_consensus if args.no_xs: cons = None sam = gem.gem2sam(map_file, index=args.index, threads=args.threads, quality=args.quality, consensus=cons, raw=raw ) gem.sam2bam(sam, output=args.output, sorted=not args.no_sort, threads=args.threads, sort_memory=str(args.sort_memory)) if not args.no_index: gem.bamIndex(args.output)
def test_quality_parameter(): assert gem._prepare_quality_parameter(None) == 'ignore' assert gem._prepare_quality_parameter(33) == 'offset-33' assert gem._prepare_quality_parameter(64) == 'offset-64'