code: def stats(args):
      
      logger.info("Reading sequeces")
      data = parse_ma_file(args.ma)
      logger.info("Get sequences from sam")
      is_align = _read_sam(args.sam)
      is_json, is_db = _read_json(args.json)
      res = _summarise_sam(data, is_align, is_json, is_db)
      _write_suma(res, os.path.join(args.out, "stats_align.dat"))
      logger.info("Done")
docs: Create stats from the analysis