[Seqan-dev] problem to use BamRecordToAlignment
- From: marion dubarry <dubarry.ma@gmail.com>
- To: seqan-dev@lists.fu-berlin.de
- Date: Tue, 21 Apr 2015 12:46:16 +0200
- Reply-to: SeqAn Development <seqan-dev@lists.fu-berlin.de>
- Subject: [Seqan-dev] problem to use BamRecordToAlignment
Hello,
I think I have a basic question, sorry in advance. I don't understand what is suppose to contain "references" of a record from the methods BamRecordToAlignment().Reference is suppose to contain : "
I start with the example and I modify it :
The string with the reference that record lies on."BamAlignmentRecord record;
StringSet<Dna5String> references;
Align<Dna5String> align;
while (!atEnd(bamFileIn))
{
readRecord(record, bamFileIn);
appendValue(references,record.seq);
if (record.rID != BamAlignmentRecord::INVALID_REFID)
bamRecordToAlignment(align, references[record.rID], record);
}
seqan-src/include/seqan/sequence/string_base.h:429 Assertion failed : static_cast<TStringPos>(pos) < static_cast<TStringPos>(length(me)) was: 1 >= 1 (Trying to access an element behind the last one!)
Aborted (core dumped)
I think the problem come from my variable references, could you explain me what is it suppose to contain ?
Thanks,
Sorry for my (stupid) question,- Follow-Ups:
- Re: [Seqan-dev] problem to use BamRecordToAlignment
- From: Jordan Willis <jwillis0720@gmail.com>
- Re: [Seqan-dev] problem to use BamRecordToAlignment
-
seqan-dev - April 2015 - Archives indexes sorted by:
[ thread ] [ subject ] [ author ] [ date ] - Complete archive of the seqan-dev mailing list
- More info on this list...