Demultiplex on i7 outer barcodes

Outer barcodes/indicies (e.g., i7 barcodes) can be used to label samples for demultiplexing, or, often to label different libraries. An example could be you and a friend both prepared separate libraries for different organisms but used the same set of 48 internal barcodes. You then both add a unique outer barcode to your samples so that you can pool them together for sequencing. The resulting data would look like below, where the i7 index for the library is in red and the internal barcode is in blue.

@NB551405:60:H7T2GAFXY:1:11101:6611:1038 1:N:0: CGAACTGT+NNNNNCNA CCGAATTCCTATCGGAAACGGATCAATAGCCCAATTGAAAATCAAGATAAGTTGAGGAAGACCAAGTCTGAAGAATTATCAAAT \+ AAAAA#EEEAEEEEE/EEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEE6EEEEA/EEEAAEEEEEEEAEEEEAEEEEAEE

i7 indices

Even if there are few samples to be demultiplexed (e.g., here we have just two, one for your library and one for your friend's) you still generally need to use multiple barcodes to create sufficient base diversity for Illumina sequencing. In this case we have 12 i7 barcodes to indicate two different libraries, so our barcodes file for i7 demultiplexing looks like the following:


In [2]:
pwd


Out[2]:
'/home/deren/Documents/toytree/manuscript/code'

demultiplex on i7 indices


In [ ]:
# demux outer i7s to plate1 and plate2
outer1 = ip.Assembly("demux_run2_i7s")
outer1.params.project_dir = "/moto/eaton/projects/PAG3"
outer1.params.raw_fastq_path = "/moto/eaton/projects/RAW_DATA/Amaranthus/Amaranthus_COL_3RAD_run2_R*.gz"
outer1.params.barcodes_path = "/moto/eaton/projects/RAW_DATA/Amaranthus/Amaranthus_COL_3RAD_run2_i7_barcodes.txt"
outer1.params.datatype = 'pairddrad'

# important: set hackers params to demux on i7 
outer1.hackersonly.demultiplex_on_i7_tags = True
outer1.hackersonly.merge_technical_replicates = True

# run step 1 to demux
outer1.run("1", ipyclient=ipyclient, force=True)

In [ ]: