The description of RNA-seq analysis spans nearly the entire paragraph, by my interpretation and limited understanding of the methods. > RNA-Seq analysis was pe

The description of RNA-seq analysis spans nearly the entire paragraph, by my interpretation and limited understanding of the methods. > RNA-Seq analysis was performed using the latest version of the human transcriptome (GRCh38_latest_rna.fna, 160,062 transcripts to which we appended the SARS-CoV-2 reference genome, MN908947). Mapping parameters were set with a mismatch cost of two, insertion and deletion cost of three, and both length and similarity fraction were set to 0.985. TPMs were generated for all 160,063 transcripts for the nine COVID-19 samples and the 40 controls (Supplementary file 2). The resulting transcript mappings for genes of interest were manually inspected to account for any expression artifacts, such as reads mapping solely to repetitive elements such as the Alu transposable element or all reads mapping to a UTR or pseudogene therein. Transcripts whose counts came solely from (or were dominated by) reads at repetitive elements were removed from the analysis. For the controls cases we ran an outlier analysis using the prcomp function in the R package factoextra. Input data were TPM for transcripts that averaged greater than one across all samples (30,102, Supplementary file 2).