Commit c01ccbfa authored by Mathieu Giraud's avatar Mathieu Giraud

core/kmerstore.h, vidjil.cpp: stores also in labels the affect for the negative strand

parent e598df7c
...@@ -213,6 +213,10 @@ void IKmerStore<T>::insert(Fasta& input, ...@@ -213,6 +213,10 @@ void IKmerStore<T>::insert(Fasta& input,
} }
labels.push_back(make_pair(T(label, 1), label)) ; labels.push_back(make_pair(T(label, 1), label)) ;
if (revcomp_indexed && ! T::hasRevcompSymetry()) {
labels.push_back(make_pair(T(label, -1), label)) ;
}
} }
template<class T> template<class T>
......
...@@ -910,6 +910,9 @@ int main (int argc, char **argv) ...@@ -910,6 +910,9 @@ int main (int argc, char **argv)
for (list< pair <KmerAffect, string> >::const_iterator it = index->labels.begin(); it != index->labels.end(); ++it) for (list< pair <KmerAffect, string> >::const_iterator it = index->labels.begin(); it != index->labels.end(); ++it)
{ {
if (it->first.getStrand() == -1)
continue ;
char key = affect_char(it->first.affect) ; char key = affect_char(it->first.affect) ;
cout << setw(12) << stats_max[key] << " " ; cout << setw(12) << stats_max[key] << " " ;
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment