BUStools / bustools

Tools for working with BUS files
https://bustools.github.io/
BSD 2-Clause "Simplified" License
92 stars 23 forks source link

Issue with bus capture RNA velocity #20

Closed mcmccon closed 5 years ago

mcmccon commented 5 years ago

Sorry just saw very similar closed issue, installed previous version and it worked fine. Please ignore/delete!

##################################################################### I am following the RNA velocity tutorial and got stumped on the bustools capture part.

First of all, I think it is missing the -s flag and I am getting this error:

Error: capture list type must be specified (one of -s, -u, or -b)
Usage: bustools text [options] bus-files

Then I added the -s flag (hopefully that is correct).

bustools capture -o introns_capture/ -s -c introns_tx_to_capture.txt -e matrix.ec -t transcripts.txt output.correct.sort.bus
Parsing transcripts .. done
Parsing ECs .. done
Parsing capture list .. done
Read in 259019682 BUS records, wrote 95361447 BUS records

But then introns_capture is empty and does not contain the split.ec needed for bustools count.

Please help me figure this out.