X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=kmer.cpp;h=50574f4a5ee0709daf2ebb67b8f0e49971e20ef0;hp=6ee84d3da3e8917c8c02d9c2d37deb68e08397b5;hb=16abd6271c455bd01b34ff89a2e3641bef0fa128;hpb=896a4f281982a3c2889f6ce6d73be997072aceae diff --git a/kmer.cpp b/kmer.cpp index 6ee84d3..50574f4 100644 --- a/kmer.cpp +++ b/kmer.cpp @@ -131,6 +131,27 @@ string Kmer::getKmerBases(int kmerNumber){ } return kmer; } +/**************************************************************************************************/ + +int Kmer::getReverseKmerNumber(int kmerNumber){ + + string kmerString = getKmerBases(kmerNumber); + + //get Reverse + string reverse = ""; + for(int i=kmerString.length()-1;i>=0;i--){ + if(kmerString[i] == 'A') { reverse += 'T'; } + else if(kmerString[i] == 'T'){ reverse += 'A'; } + else if(kmerString[i] == 'G'){ reverse += 'C'; } + else if(kmerString[i] == 'C'){ reverse += 'G'; } + else { reverse += 'N'; } + } + + int reverseNumber = getKmerNumber(reverse, 0); + + return reverseNumber; + +} /**************************************************************************************************/