Skip to content

Commit 454446b

Browse files
committed
BUG: properly slices qual scores
1 parent 099278b commit 454446b

File tree

1 file changed

+4
-4
lines changed

1 file changed

+4
-4
lines changed

qiita_ware/demux.py

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -438,7 +438,7 @@ def to_ascii(demux, samples=None):
438438
for samp, idx, seq, qual, bc_ori, bc_cor, bc_err in fetch(demux, samples):
439439
seq_id = id_fmt % {'sample': samp, 'idx': idx, 'bc_ori': bc_ori,
440440
'bc_cor': bc_cor, 'bc_diff': bc_err}
441-
if qual is not None:
441+
if qual != []:
442442
qual = qual.astype(np.uint8)
443443

444444
yield formatter(seq_id, seq, qual)
@@ -516,7 +516,7 @@ def fetch(demux, samples=None, k=None):
516516
seqs = demux[pjoin(dset_paths['sequence'])][indices]
517517

518518
# only yield qual if we have it
519-
quals = repeat(None)
519+
quals = repeat([])
520520
if demux.attrs['has-qual']:
521521
if len(indices) == 1:
522522
if indices[0]:
@@ -531,8 +531,8 @@ def fetch(demux, samples=None, k=None):
531531
iter_ = zip(repeat(sample), np.arange(indices.size)[indices], seqs,
532532
quals, bc_original, bc_corrected, bc_error)
533533

534-
for item in iter_:
535-
yield item
534+
for samp, idx, seq, qual, bc_ori, bc_cor, bc_err in iter_:
535+
yield (samp, idx, seq, qual[:len(seq)], bc_ori, bc_cor, bc_err)
536536

537537

538538
def stats(demux):

0 commit comments

Comments
 (0)