Skip to content

Influenza data pipeline to automate genotyping assignment

License

Notifications You must be signed in to change notification settings

USDA-VS/GenoFLU

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

17 Commits
 
 
 
 
 
 
 
 
 
 

Repository files navigation

GenoFLU

This tool uses BLAST to identify North American H5NX genomes in the 2.3.4.4b clade from a curated database. Pre-defined genotypes are cross-referenced with the top segment identifications, and a genotype is assigned. A cutoff of 2% difference from the closest curated sequence identifies new reassortment. New reassortment is reviewed using segment-based phylogenetic trees. If appropriate, new segment sequences will be added to the curated database and new genotype assignments updated.

Installation

conda create -c conda-forge -c bioconda -n genoflu genoflu

Usage

FASTA file containing a single segmented influenza genome, with each segment having its own individually named header.

genoflu.py -f <*.fasta>

Output

Genotype summary as Excel and tab delimited text file.

Test

Test genome available at test/test-genome-A1.fasta

genoflu.py -f test-genome-A1.fasta

test-genome-A1 Genotype --> A1: PB2:ea1, PB1:ea1, PA:ea1, HA:ea1, NP:ea1, NA:ea1, MP:ea1, NS:ea1

About

Influenza data pipeline to automate genotyping assignment

Resources

License

Stars

Watchers

Forks

Packages

No packages published

Languages