Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Feature/allele calling #17

Merged
merged 45 commits into from
Apr 11, 2024
Merged

Conversation

luissian
Copy link
Member

@luissian luissian commented Mar 17, 2024

Implementation of allele calling
In this PR, when it is searching the specif allele that match the blast sequence using grep, if grep match then this specific allele name will be set as the field of coding. This is to find the EXC allele, but could be set in conditions of ASM or ALM or PLOT, where part of the sequence match in the locus file.

  • save all coordinates in case NIPHEM/NIPH
  • if no blast result are with 80% query. it is considered from 75% and assigned to PLOT
  • check if programs like mash, prokka, etc exist before start execution
    • mash
    • prokka
    • blast
  • create SNP file all alleles compared with their Reference allele. including 3 DNA bases for ref and alt, Amino acid for ref and alt and change type of amino acid
  • Create multi alignment for each allele with the reference allele
  • create summary graphic for all samples and all allele classification
  • update code to process all samples
  • ensure that the get_inferred_allele method is accessed safely by multiple instances
  • Extend blast sequence for searching start and stop codon due to the variability in reference alleles
  • Parallel execution to create the multi alignment files

Important Note .-
This PR includes changes in analyze schema feature. This means it is mandatory to re-run the analyze schema to get the right annotation file

@luissian luissian marked this pull request as ready for review March 22, 2024 07:50
taranis/__main__.py Outdated Show resolved Hide resolved
taranis/__main__.py Outdated Show resolved Hide resolved
taranis/__main__.py Outdated Show resolved Hide resolved
taranis/allele_calling.py Show resolved Hide resolved
taranis/allele_calling.py Outdated Show resolved Hide resolved
taranis/allele_calling.py Outdated Show resolved Hide resolved
taranis/allele_calling.py Outdated Show resolved Hide resolved
taranis/allele_calling.py Outdated Show resolved Hide resolved
taranis/allele_calling.py Outdated Show resolved Hide resolved
taranis/allele_calling.py Outdated Show resolved Hide resolved
@luissian luissian force-pushed the feature/allele_calling branch from fe10c93 to 2821301 Compare April 11, 2024 11:28
@saramonzon saramonzon merged commit a2ffba4 into BU-ISCIII:develop Apr 11, 2024
3 checks passed
@luissian luissian deleted the feature/allele_calling branch April 11, 2024 21:04
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants