/*
* 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 |
Ensembl로 Gene symbol 알아내기 (0) | 2015.01.03 |
SAM/BAM 포맷, 구조 분석 (0) | 2014.12.12 |
NGS vs Sanger sequencing (10) | 2014.10.21 |