I am trying to use fastX barcode splitter perl script on my fastq file. However, every time I run the program I get this message Error: barcode file not specified (use '--bcfile [FILENAME]'). All I am doing is dragging the perl script into terminal, and dragging the fastq file into the terminal and hit enter and that what i get. I also read I had to make mybarcode text file which is BC1 and then tab space then the barcode where does that go. Any help would be appreciated
You need to use the actual UNIX command structure instead of dragging and dropping, and you need to provide the barcode text file.
Make a barcode text file just as you specified, either in textedit (sounds like you are on a Mac) or within the terminal, using something like
nano. If you want to make a text file like this using
vi, follow these direction exactly:
## This opens the built-in text editor
## To type, you need to switch to "insert mode" by pressing the
## Type the following, as you specified:
BC1 ATCGTT ## BC1 is the name you are giving the barcode followed by the barcode ID, separated by a tab
## To save and close the document, press
esc to exit "insert mode", and type
:wq to write the file and quit
Now you will have the barcode file you need! So then run your command like this:
cat my_raw_read_file.fastq | fastx_barcode_splitter.pl --bcfile my_barcode.txt
Please keep in mind there are many parameters that you may want to use, so please thoroughly check out the explanation and usage examples found in the documentation: http://hannonlab.cshl.edu/fastx_toolkit/commandline.html#fastx_barcode_splitter_usage
I also have similar question. I have 56 sample and barcodes. How should i make a barcode file, should i make a file for individual barcode or i can make a single file for 56 sample using tab delimited text formate. I made a single barcode file with 56 barcode and try to split the original FASTQ file into 56 small FASTQ file based on barcode but everytime i am getting only one file with first barcode and all the reads match with first barcode. Could you please suggest me is there any away to split it into 56 file.