combine_kreports.py: Latest version missing unclassified reads for last sample

[sfredr]

Hi,
It looks like the unclassified read counts are always lost for the last sample in the most recent version of combine_kreports.py?

python combine_kreports.py --report-file $WD/kraken_report/* --output $WD/kraken_table.txt

Previous version (50b081e), everything looks good:

#perc	tot_all	tot_lvl	1_all	1_lvl	2_all	2_lvl	3_all	3_lvl	lvl_type	taxid	name
3.3133	6507	6507	713	713	190	190	5604	5604	U	0	unclassified
96.6867	189885	17	41811	1	75544	0	72530	16	R	1	root
96.676	189864	803	41810	74	75544	109	72510	620	D	3	  Bacteria

Current version (8e009c2), read counts missing and replaced by values belonging in taxid and name columns:

#perc	tot_all	tot_lvl	S1_all	S1_lvl	S2_all	S2_lvl	S3_all	S3_lvl	lvl_type	taxid	name
3.3133	6507	6507	713	713	190	190	U	0	unclassified		
96.6867	189885	17	41811	1	75544	0	72530	16	R	1	root
96.676	189864	803	41810	74	75544	109	72510	620	D	3	  Bacteria

[jenniferlu717]

Thanks for letting me know. I'll fix it now

[sfredr]

Thank you for the quick response! It works as expected now.