in reply to Re^2: Deconvolutinng FastQ files
in thread Deconvolutinng FastQ files
to:unless $barcode =~ m/^[AGCT]+$/;
then you should be able to prefix your barcodes w/ 3 N's as long as you set --mismatches to at least 3 on the command line when running the script.unless $barcode =~ m/^[AGCTN]+$/;
One caveat is that you will want to toss out any reads that have any Ns in the first X bases (where X = 3+ barcode length). Have you run FastQC? If so, this will tell you the per base N content. It probably won't be an issue if you've already done preliminary filtering based on Illumina's Y/N flags (assuming Illumina sequencing, of course).
Also, (depending on your computer, of course) I suspect fastx_barcode_splitter.pl will run a lot faster at the CLI than on Galaxy (at least if you are using the public galaxy server).
Edit: to avoid the potential problem w/ Ns, just use some other non-nucleotide character!
|
|---|
| Replies are listed 'Best First'. | |
|---|---|
|
Re^4: Deconvolutinng FastQ files
by snakebites (Initiate) on Aug 08, 2012 at 02:53 UTC |