Re: [Seqan-dev] SAM issue?
- From: "Holtgrewe, Manuel" <manuel.holtgrewe@fu-berlin.de>
- To: SeqAn Development <seqan-dev@lists.fu-berlin.de>
- Date: Mon, 21 Mar 2011 13:12:45 +0100
- Acceptlanguage: en-US, de-DE
- Reply-to: SeqAn Development <seqan-dev@lists.fu-berlin.de>
- Subject: Re: [Seqan-dev] SAM issue?
Your problem is that you use an alignment id as the read id. See my proposal for a fixed version attached that yields consistent results for me (regarding the positions, there still is the bug using alignment ids for read ids so the printed read names are wrong). $ ./demos/SAMtest -s bug2.sam Loading SAM file bug2.sam Reads cached... n:2825 id: 1022 illumina_80bp_3kb.000000199 19625,19705 $ ./demos/SAMtest -s bug3.sam Loading SAM file bug3.sam Reads cached... n:3253 id: 1022 illumina_80bp_3kb.000000170 19625,19705 $ ./demos/SAMtest -s bug4.sam Loading SAM file bug4.sam Reads cached... n:3648 id: 1022 illumina_80bp_3kb.000000148 19625,19705 Am 21.03.2011 um 12:52 schrieb Mat: > Hi! > > So if i adapt the code i still get the weired results: > > test.sam: > > illumina_80bp_3kb.000000571 > 2128,2048 > > > bug2.sam (70000 lines of test.sam): > > illumina_80bp_3kb.000000571 > 237,157 > > (code attached...) > > cheers > <SAMtest.cpp><ATT00001..txt>
Attachment:
SAMtest.cpp
Description: SAMtest.cpp
- Follow-Ups:
- Re: [Seqan-dev] SAM issue?
- From: Mat <matthias.dodt@mdc-berlin.de>
- Re: [Seqan-dev] SAM issue?
- References:
- [Seqan-dev] SAM issue?
- From: Mat <matthias.dodt@mdc-berlin.de>
- Re: [Seqan-dev] SAM issue?
- From: Manuel Holtgrewe <manuel.holtgrewe@fu-berlin.de>
- Re: [Seqan-dev] SAM issue?
- From: Mat <matthias.dodt@mdc-berlin.de>
- Re: [Seqan-dev] SAM issue?
- From: Manuel Holtgrewe <manuel.holtgrewe@fu-berlin.de>
- Re: [Seqan-dev] SAM issue?
- From: Mat <matthias.dodt@mdc-berlin.de>
- Re: [Seqan-dev] SAM issue?
- From: "Holtgrewe, Manuel" <manuel.holtgrewe@fu-berlin.de>
- Re: [Seqan-dev] SAM issue?
- From: Mat <matthias.dodt@mdc-berlin.de>
- [Seqan-dev] SAM issue?
-
seqan-dev - March 2011 - Archives indexes sorted by:
[ thread ] [ subject ] [ author ] [ date ] - Complete archive of the seqan-dev mailing list
- More info on this list...