I am trying to use fread to read the genome alignment in data.tablein R. Here is a snapshot of the alignment file:
USI-EAS28:1:100:1786:674
USI-EAS28:1:100:1786:940
USI-EAS28:1:100:1786:705
USI-EAS28:1:100:1786:1168
USI-EAS28:1:100:1787:582
USI-EAS28:1:100:1787:62
USI-EAS28:1:100:1788:573
USI-EAS28:1:100:1788:854
USI-EAS28:1:100:1788:185
USI-EAS28:1:100:1788:1339
The second column reports that the line that it is reading is aligned ( +- forward, -- reverse). Unfortunately, fread tries to read this column as an integer, always assigning a value of 0. This column should be read as a character or even logical, for that matter. Trying to play with arguments sepand sep2does not help.
source
share