Algorithm
Here is the algorithm by which Metagenomic Thermometer calculates MPT
Quality filtering
Adapter trimming and removal of low quality reads were performed using Fastp v0.20.0. The lower limit of the read length was set to 50 bp, and the allowable limit of the number of ambiguous nucleotides (N) per read was set to 1
Converting fastq to fna
Fastq file is converted to DNA fasta file using SeqKit v0.16.1.
Gene prediction
ORF regions were identified from DNA fasta files using Prodigal v2.6.3, and amino acid fasta files of ORF regions were obtained. Prodigal is run with "-p meta" options, otherwise with default parameters.
MPT calculation
First, calculate the proportion of IVYWREL in the total amino acids contained in the amino acid sequence.Then, Metagenomic Predicted Temperature is calculated by following equation.
