Re: [Seqan-dev] SAM issue?


It should: An insertion in the read regarding the reference inserts a gap which shifts all positions in the alignment right of it by one.

Am 21.03.2011 um 12:37 schrieb Mat:

> Ahh okay - well that might be a problem. But still i think the output shouldn't change with the number of lines in the sam file?!
> 
> thanks!
> 
> Am 3/21/11 12:14 PM, schrieb Holtgrewe, Manuel:
>> ontigStore                        TContigStore;
>> typedef typename Value<TContigStore>::Type                           TContig;
>> typedef typename TFragmentStore::TContigSeq                          TContigSeq;
>> typedef Gaps<TContigSeq, AnchorGaps<typename TContig::TGapAnchors> > TContigGaps;
>> 
>> typedef typename TFragmentStore::TAlignedReadStore                   TAlignedReadStore;
>> typedef typename Value<TAlignedReadStore>::Type                      TAlignedRead;
>> typedef typename TAlignedRead::TPos                                  TAlignedReadPos;
>> 
>> TContigGaps contigGaps(contigStore[contigId].seq, contigStore[contigId].gaps);
>> TAlignedRead const & alignedRead = alignedReadStore[idx];
>> // Translate end position from aligne
>> 
> <ATT00001..txt>