Enter a genome
Function: band pass(sqrt(10 *(numbers('yxr') if 100 *(numbers('yxr')) > numbers("xxy") else numbers('rep'))))
Enter another genome