Re: [Seqan-dev] Trouble Hashing Kmers into Qgram Index

[Rahn, René <rene.maerker@fu-berlin.de>]

Hey,

the problem lies in the hash(indexShape(index), kmer) function.

The hash function only takes an iterator or pointer as second argument. 

So in your case the simple solution would be to call the function hash(indexShape(index), it), while using an iterator to iterate over the text instead of using the unsigned integer. This little behaviour does not seem obvious, as in the tutorial you pass a string literal. But in fact a string literal is just a const char *, so it works.

Here you see an example code to work with a) using the finder interface and b) directly operating on the index.

IHTH.

cheers,

René 

    typedef Index<DnaString, IndexQGram<UngappedShape<12> > > TQGramIndex;

    typedef Finder<TQGramIndex> TFinder;

    typedef Iterator<DnaString, Standard>::Type TIterator;

    typedef Fibre<TQGramIndex, QGramShape>::Type TShape;

    TQGramIndex index(ref);

    indexRequire(index, QGramSADir());  // On-demand index creation.

    // a) Using the finder interface.

    TFinder finder(index);

    for (TIterator it = begin(query, Standard()); it != end(query, Standard()) - 12; ++it)

    {

        std::cout << "Occ at: ";

        while(find(finder, infix(query, it, it+12)))

            std::cout << position(finder) << " ";

        std::cout << std::endl;

        clear(finder);  // Clear finder for next search.

    }

    // b) Using the index interface.

    TShape & shape = indexShape(index);

    hashInit(shape, begin(query, Standard()));

    for (TIterator it = begin(query, Standard()); it != end(query, Standard()) - 12; ++it)

    {

        std::cout << "Occ at: ";

        hashNext(shape, it);

        for (unsigned i = 0; i < length(getOccurrences(index, shape)); ++i)

            std::cout << getOccurrences(index, shape)[i] << " ";

        std::cout << std::endl;

    }

On 19 Nov 2014, at 07:34, Brett Bowman <bnbowman@gmail.com> wrote:

Sorry to be a bother, but I decided to go back to the docs and start fresh and I'm still having issues.  

We're using Qgrams, so I trying to write a simple loop to search for all 12mers from a query sequence "seq" in the index formed from a StringSet of references "refSeq", building on the example code in the Q-Gram Index tutorial:

"""

  TFinder qgramFinder(index);

 for (size_t i = 0; i < length(seq)-12; i++)

 {

     TInfix::Type kmer = infix(seq, i, i+12);

     std::cout << kmer << std::endl;

     seqan::hash(indexShape(index), kmer);

     for (unsigned i = 0; i < length(getOccurrences(index, indexShape(index))); ++i)

     {

         std::cout << getOccurrences(index, indexShape(index))[i] << std::endl;

     }

     std::cout << std::endl;

 }

"""

but I'm getting a fairly obtuse error message I don't quite understand, and doesn't point me back to any specific lines in my code.

"""

/usr/include/seqan/index/shape_base.h:535:55: error: indirection requires pointer operand ('seqan::Segment<const seqan::String<seqan::SimpleType<unsigned

      char, seqan::Dna5_>, seqan::Alloc<void> >, seqan::InfixSegment>' invalid)

                        hash * ValueSize<TValue>::VALUE + ordValue((TValue)*it),

                                                                           ^~~

/usr/include/seqan/index/shape_base.h:548:22: note: in instantiation of function template specialization 'seqan::_hashFixedShape<unsigned long,

      seqan::SimpleType<unsigned char, seqan::Dna5_>, seqan::Segment<const seqan::String<seqan::SimpleType<unsigned char, seqan::Dna5_>, seqan::Alloc<void> >,

      seqan::InfixSegment>, 12>' requested here

                return me.hValue = _hashFixedShape(me.hValue, it, TValue(), UngappedShape<q>());

                                   ^

/home/bbowman/git/SRSLI/src/C++/SparseAlignment.hpp:102:16: note: in instantiation of function template specialization 'seqan::hash<seqan::SimpleType<unsigned

      char, seqan::Dna5_>, 12, seqan::Segment<const seqan::String<seqan::SimpleType<unsigned char, seqan::Dna5_>, seqan::Alloc<void> >, seqan::InfixSegment>

      >' requested here

        seqan::hash(indexShape(index), kmer);

               ^

/home/bbowman/git/SRSLI/src/C++/main.cpp:86:9: note: in instantiation of function template specialization 'FindSeeds2<FindSeedsConfig<12,

      seqan::UngappedShape<12>, seqan::IndexQGram<seqan::UngappedShape<12>, seqan::Tag<seqan::Default_> > > >' requested here

        FindSeeds2(querySeedHits, refSetIndex, idxAndRecord.second.Seq);

"""

Suggestions?

Sincerely,

Brett

_______________________________________________
seqan-dev mailing list
seqan-dev@lists.fu-berlin.de
https://lists.fu-berlin.de/listinfo/seqan-dev

---

René Rahn

Ph.D. Student

--------------------------------

Tel:  (+49) 30 838 75277

Mail: rene.rahn@fu-berlin.de

--------------------------------

Institute of Computer Science

Algorithmic Bioinformatics (ABI)

--------------------------------

Freie Universität Berlin

Takustraße 9

14195 Berlin

--------------------------------