Beefy Boxes and Bandwidth Generously Provided by pair Networks
Perl: the Markov chain saw

Re^3: Deconvolutinng FastQ files

by frozenwithjoy (Priest)
on Aug 07, 2012 at 15:35 UTC ( #986019=note: print w/replies, xml ) Need Help??

in reply to Re^2: Deconvolutinng FastQ files
in thread Deconvolutinng FastQ files

I took a look at and I think I've figured out a solution. I haven't tested it, but if you change line 161 from:
unless $barcode =~ m/^[AGCT]+$/;
unless $barcode =~ m/^[AGCTN]+$/;
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.

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 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
    Thank you frozenwithjoy. I should give this a go too. We are thinking about running our own Galaxy server in the EC2, so the revised fastx_barcode_splitter might come in handy. Browseruk's script below works super fast, I am not sure how might compare.

Log In?

What's my password?
Create A New User
Node Status?
node history
Node Type: note [id://986019]
and all is quiet...

How do I use this? | Other CB clients
Other Users?
Others chanting in the Monastery: (11)
As of 2018-01-23 14:19 GMT
Find Nodes?
    Voting Booth?
    How did you see in the new year?

    Results (247 votes). Check out past polls.