I am using maq's map program to map reads to a reference genome. I would like to use only the first 37bp of my reads (long story). I have used the -1 option, but I'm not sure I am using it correctly. I am using the following shell script:
maq map -1 37 . . .
Anyone have experience with the -1 option?
maq map -1 37 . . .
Anyone have experience with the -1 option?