receptor-utils 0.0.1

Some tools I find useful for working with Ig receptor sequences

receptor_utils

Some tools I find useful for working with Ig receptor sequences.

Installation

git clone https://github.com/williamdlees/receptor_utils
pip install receptor_utils

The module requires Biopython.

(will be on PyPi soon)

Overview

Please refer to the files themselves for slightly more detailed documentation.

simple_bio_seq

Contains some convenience functions that are backed by BioPython but simplified for my use case. It uses the following approach to keep things simple (at the expense of some flexibility/scalability):

• store sequences as strings, use dicts for collections
• convert sequences to upper case on input
• coerce iterators into lists for ease of debugging

from receptor_utils import simple_bio_seq as simple
seqs = simple.read_fasta('seqfile.fasta')  # read sequences into a dict with names as keys
seq = simple.read_single_fasta('seqfile.fasta')  # reads the first or only sequence into a string
seq = simple.reverse_complement(seq)

See the file for other functions.

novel_allele_name

Contains the function name_novel(), which will generate a name for a 'previously undocumented' allele, given its sequence. The name will consist of the name of the nearest allele in a reference set provided to the function, suffixed by the SNPs that differentiate it, for example:

IGHV1-69*01_a29g_c113t

Numbering of V-sequences uses the IMGT alignment. The naming convention follows that used by Tigger and VDJbase.

number_ighv

Contains various functions for working with V-sequences according to the IMGT numbering scheme. The most useful is gap_sequence() which will gap the provided V-sequence by using the closest sequence in a reference set as a template.

Example scripts

These may be useful in their own right, but also show how to use some of the functions mentioned above.

extract_refs.py

A script which uses simple_bio_seq to extract files for particular loci and species from an IMGT reference file.

gap_inferred.py

A script which will gap a set of sequences listed in a FASTA file, using the closest sequences discovered from a reference set.

identical_seqs.py

A script which uses simple_bio_seq to list identical sequences and sub-sequences in a fasta file.