From b746564993c8e5635793641278c72a0e9ce314a1 Mon Sep 17 00:00:00 2001 From: Zengtudor Date: Wed, 18 Sep 2024 16:44:17 +0800 Subject: [PATCH] update --- src/main.cpp | 56 +++++++++++++++++++++++++++++++++++----------------- 1 file changed, 38 insertions(+), 18 deletions(-) diff --git a/src/main.cpp b/src/main.cpp index c757e76..4ced0dd 100644 --- a/src/main.cpp +++ b/src/main.cpp @@ -12,6 +12,7 @@ #include"tools.hpp" #include #include +#include // std::string reversedComplement(std::string DNAsequence); @@ -33,22 +34,38 @@ void reverseComplement(std::array &DNAsequence,const size_t buf_ for (i = 0; i < q; i++) { char &l = DNAsequence[i]; - if (DNAsequence[i] == 'A') - l = l + 'T'; - if (DNAsequence[i] == 'a') - l = l + 'T'; - if (DNAsequence[i] == 'c') - l = l + 'G'; - if (DNAsequence[i] == 'C') - l = l + 'G'; - if (DNAsequence[i] == 'g') - l = l + 'C'; - if (DNAsequence[i] == 'G') - l = l + 'C'; - if (DNAsequence[i] == 't') - l = l + 'A'; - if (DNAsequence[i] == 'T') - l = l + 'A'; + if (DNAsequence[i] == 'A'){ + // l = l + 'T'; + l = 'T'; + } + if (DNAsequence[i] == 'a'){ + // l = l + 'T'; + l = 'T'; + } + if (DNAsequence[i] == 'c'){ + // l = l + 'G'; + l = 'G'; + } + if (DNAsequence[i] == 'C'){ + // l = l + 'G'; + l = 'G'; + } + if (DNAsequence[i] == 'g'){ + // l = l + 'C'; + l = 'C'; + } + if (DNAsequence[i] == 'G'){ + // l = l + 'C'; + l = 'C'; + } + if (DNAsequence[i] == 't'){ + // l = l + 'A'; + l = 'A'; + } + if (DNAsequence[i] == 'T'){ + // l = l + 'A'; + l = 'A'; + } } } @@ -74,12 +91,15 @@ int main() { int m = lines%2; const auto buf_len = strlen(buf.data()); - + const std::string_view suffix("\n"); if (m == 1){ // output_file_stream << reverseComplement(buf) << endl; reverseComplement(buf,buf_len); } - output_file_stream.write(buf.data(), buf_len); + for(int i=0;i