Use bam seq len, for bounded sequence.

natir · natir · commit 26d7d3aed9de · 2017-04-23T16:37:24.000+02:00
diff --git a/sequence/bamreader.cpp b/sequence/bamreader.cpp
@@ -33,8 +33,8 @@ bool BamReader::next()
             qseq[i] = seq_nt16_str[bam_seqi(q,i)]; //gets nucleotide id and converts them into IUPAC id.
         }
 
-        seq.setSequence(QByteArray(qseq));
-        seq.setQuality(QByteArray(qseq));
+        seq.setSequence(QByteArray(qseq, len));
+        seq.setQuality(QByteArray(qseq, len));
 
 
         // update sequence

Original file line number	Diff line number	Diff line change
`@@ -33,8 +33,8 @@ bool BamReader::next()`
`33`	`33`	`qseq[i] = seq_nt16_str[bam_seqi(q,i)]; //gets nucleotide id and converts them into IUPAC id.`
`34`	`34`	`}`
`35`	`35`
`36`		`- seq.setSequence(QByteArray(qseq));`
`37`		`- seq.setQuality(QByteArray(qseq));`
	`36`	`+ seq.setSequence(QByteArray(qseq, len));`
	`37`	`+ seq.setQuality(QByteArray(qseq, len));`
`38`	`38`
`39`	`39`
`40`	`40`	`// update sequence`