script to analyse QUASAR sequence