Effectively and efficiently diagnosing COVID-19 patients with accurate clinical type is essential to achieve optimal outcomes of the patients as well as reducing the risk of overloading the healthcare system. Currently, severe and non-severe COVID-19 types are differentiated by only a few clinical features, which do not comprehensively characterize complicated pathological, physiological, and immunological responses to SARS-CoV-2 invasion in different types. In this study, we recruited 214 confirmed COVID-19 patients in non-severe and 148 in severe type, from Wuhan, China. The patients' comorbidity and symptoms (26 features), and blood SERO biochemistry (26 features) upon admission were acquired as two input modalities. Exploratory analyses demonstrated that these features differed substantially between two clinical types. Machine learning random forest (RF) models using features in each modality were developed and validated to classify COVID-19 clinical types. Using comorbidity/symptom and biochemistry as input independently, RF models achieved >90% and >95% predictive accuracy, respectively. Input features' importance based on Gini impurity were further evaluated and top five features from each modality were identified ( age TRANS, hypertension MESHD hypertension HP, cardiovascular disease MESHD, gender TRANS, diabetes; D-Dimer, hsTNI, neutrophil, IL-6, and LDH). Combining top 10 multimodal features, RF model achieved >99% predictive accuracy. These findings shed light on how the human body reacts to SARS-CoV-2 invasion as a unity and provide insights on effectively evaluating COVID-19 patient's severity and developing treatment plans accordingly. We suggest that symptoms and comorbidities can be used as an initial screening tool for triaging, while biochemistry and features combined are applied when accuracy is the priority.