strandsalgo.cc

Go to the documentation of this file.

#include "strandsalgo.h"
#include <qinputdialog.h>
#include <string>

using namespace std;

void StrandsAlgo::start()
{
  QStringList lst;
  lst << "Normal" << "Reverse";
  bool ok;
  
  QString res = QInputDialog::getItem ( "Select reads on strand", 
                                        "Select strand, normal or reverse", 
                                        lst, 
                                        0, FALSE, 
                                        &ok, 0, 0 );
  
  if ( !ok ) return;
  
  string strand;
  
  if ( res == "Normal" ) {
    strand = "U";
  }
  else {
    strand = "C";
  }
  
  
  for( size_t i = 0; i < getMAl()->get_num_seq(); i++ ) {
    
    if ( getMAl()->get_strand(i) != strand ) {
      getMAl()->select_read(i, false);
    }
  }
}

---