/*

 http://sosal.kr/
 * made by so_Sal
 */



서점에서 지나가는듯 Bioinformatics 책을 봤는데

역상보서열 만들기 주제로 단원이 하나 있길래 심심해서 저도 만들어봅니다 ㅎㅎ.


커맨드 인수 argv를 통해서 input을 받고, 해당 서열의 역상보 서열을 출력합니다.




#include<iostream>

#include<string.h>

#include<stdlib.h>

using namespace std;

 

string convert(string str);

 

int main(int argc, char* argv[]){

 

    string str(argv[1]);

    string revcomp = convert(str);

 

    cout<<"Input  sequence: "<<str<<endl;

    cout<<"Output sequence: "<<revcomp<<endl;

}

 

string convert(string str){

    for(int i=0;i<str.length();i++){

        switch( str[i] ){

            case 'A':

                str[i] = 'T';

                break;

            case 'C':

                str[i] = 'G';

                break;

            case 'G':

                str[i] = 'C';

                break;

            case 'T':

                str[i] = 'A';

                break;

            default:

                cout<<"Argument error."<<endl;

                exit(-1);

        }

    }

 

    return str;

}






뭐 너무 쉬워서 부가설명 할게 없네요..


'Major Study. > Bioinformatics' 카테고리의 다른 글

Transition vs Transversion  (0) 2015.02.12
DNAI1 - ENSG, Variant Interpretation  (0) 2015.01.28
주어진 DNA의 역상보서열 만들기  (0) 2015.01.09
Ensembl로 Gene symbol 알아내기  (0) 2015.01.03
SAM/BAM 포맷, 구조 분석  (0) 2014.12.12
NGS vs Sanger sequencing  (10) 2014.10.21
Posted by sosal sosal

댓글을 달아 주세요