Re: [Seqan-dev] Dna5Q - how do you access qual or seq ?

[Leon Kuchenbecker <lkuchenb@inf.fu-berlin.de>]

Hi Stephen!

> //where seqs is the StringSet of Dna5Q read in above
> seqan::Dna5Q seq= seqs[0];
> 
> Rcpp::Rcout << getQualityValue(seq[0]) << std::endl;    // error
> [...]

seqan::Dna5Q is a type for storing a single character from the ACGTN
alphabet including the associated quality. If you want to store a
sequence of such characters, you need String<Dna5Q>.

Thus, the correct types for your variables would be

> seqan::StringSet<seqan::CharString> ids;
> seqan::StringSet<seqan::String<seqan::Dna5Q> > seqsQ;

Your usage of getQualityValue() was in principle correct, just that you
were trying to subscript a single character.

Cheers
  Leon

-- 
You can use my PGP pubkey (0x010E29F5) to encrypt mail addressed to me
()  ascii ribbon campaign - against html e-mail
/\  www.asciiribbon.org   - against proprietary attachments

Attachment: smime.p7s
Description: S/MIME Cryptographic Signature