4
$\begingroup$

We have a problem trying to demultiplex MinION sequences with custom barcodes. Do you have any software recommendations we can try for demultiplexing or how to demultiplex these custom barcodes with Albacore? We have tried using albacore but it only recognizes barcodes provided by Nanopore Technologies.

Improve this question
$\endgroup$
0
3
$\begingroup$

My semi-automatic workflow for demultiplexing nanopore reads using LAST can be found here - pay attention to the comments, where I've identified bits in the steps which might trip up people trying to follow them:

https://dx.doi.org/10.17504/protocols.io.7vmhn46

While this protocol is designed around using the ONT barcodes, the barcode fasta file can be replaced by another file [using barcode ID in the header] without changing the protocol. Make sure that this file only includes the unique parts of the barcodes (with a few bases of context if necessary to exclude random matches), so that reads won't be matched to the non-unique part of the wrong barcode sequence.

Improve this answer
$\endgroup$
3
  • $\begingroup$ gringer, thanks for the workflow. It would tremendously help if you included a link to where to find LAST. Thanks! $\endgroup$
    – Lev Yampolsky
    Dec 2 '19 at 0:00
  • $\begingroup$ Thanks for the comment; I've added a link in the comments of the protocol, and in the answer for this question. I'll put a LAST link in the method once I do my next update of the protocol. $\endgroup$
    – gringer ♦
    Dec 2 '19 at 20:10
  • $\begingroup$ Thank you! I am using minibar, works great and very easy to add custom barcodes. L.Y. $\endgroup$
    – Lev Yampolsky
    Dec 4 '19 at 2:28
3
$\begingroup$

There are two other ways in which you could demultiplex reads if guppy doesn't have the config file for the kit you are using.

  1. Minibar: This works for dual index barcodes specifically
  2. Porechop: This is intended for the ONT sets, but you can also replace them with your own. However, it is designed to work with single index barcodes. That means, you will have to run it once with your forward indices, then take that output and run again with the reverse, and then do some gymnastics to combine it all together
Improve this answer
$\endgroup$
1
  • $\begingroup$ What do you mean with run it twice, with forward and reverse indices? $\endgroup$
    – F.Lira
    Apr 20 at 7:35

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.