SnapTools
SnapTools copied to clipboard
Bug in isBedQuerynameSorted()
Hi,
I tried to use snap_pre
to create snap
object from BED files.
However, it always reports as error: ./test.bed is not a sorted by read name!
, even I have sorted them.
I see that, in the code of the function isBedQuerynameSorted(), at line 166
cur_barcode = str(item).split()[3];
, perhaps this line is wrong, and it causes the error.
Should the string be split the same way as in the function getBarcodesFromBed(input_bed), line 256?
barcode = _read.split()[3].split(":")[0].upper();
Please check. Thank you very much.