We present a statistical method to estimate reliability thresholds for conserved functional regions. Single stranded positions and base pair positions are collected that have a high evolutionary reliability. We write down only the base pair part. Single-stranded positions are treated analogously. For this purpose, do the following 1. Generate shuffled alignment A shuffle by shuffling the alignment columns. Then, we generate again the most likely structure under the shuffled alignment, i.e., we generate