mojaveazure / snp_utils Goto Github PK
View Code? Open in Web Editor NEWExtract SNP positions from BLAST XML and SAM formats
Extract SNP positions from BLAST XML and SAM formats
There is a bug, which can not take care of the gap in the reference. So some of the physical positions are wrong.
The examples I found are here: /panfs/roc/scratch/llei/wrong_position_SNPuti
I'm trying to run snp_utils on Ubuntu 22.04 and python 3.10.6.
It generates the xml output file but later it stop with errors.
Below are the last lines on the shell. I cannot figure out what is happening.
No query found in GMI_ES02_c16875_166:8.60053e-64
hmm GMI_ES03_c17575_451:0.0
No reference allele for GMI_ES03_c17575_451:0.0, have you run Hsp.get_snp_position() yet?
hmm GMI_ES03_c17575_451:1.12148e-21
No query found in GMI_ES03_c17575_451:1.12148e-21
hmm GMI_ES03_c17575_451:3.55129e-111
No query found in GMI_ES03_c17575_451:3.55129e-111
hmm GMI_ES03_c17575_451:5.21775e-20
No query found in GMI_ES03_c17575_451:5.21775e-20
hmm GMI_ES03_c17575_451:5.21775e-20
No query found in GMI_ES03_c17575_451:5.21775e-20
hmm GMI_ES03_c17370_63:1.92383e-137
Traceback (most recent call last):
File "/home/davide/Programs/SNP_Utils-master/snp_utils.py", line 285, in
main()
File "/home/davide/Programs/SNP_Utils-master/snp_utils.py", line 262, in main
snp_list, no_snps, ref_gen, bconf = blast_based(args, lookup_dict)
File "/home/davide/Programs/SNP_Utils-master/snp_utils.py", line 103, in blast_based
hsp.add_snp(lookup=lookup)
File "/home/davide/.local/lib/python3.10/site-packages/overload.py", line 181, in f
return callable(*usable_args, **_kw)
File "/home/davide/Programs/SNP_Utils-master/Objects/blast.py", line 267, in add_snp
self.add_snp(this_snp=snp.SNP(lookup=lookup, hsp=self))
File "/home/davide/.local/lib/python3.10/site-packages/overload.py", line 181, in f
return callable(*usable_args, **_kw)
File "/home/davide/Programs/SNP_Utils-master/Objects/snp.py", line 103, in init
self._position = hsp.get_snp_position(lookup=lookup)
File "/home/davide/.local/lib/python3.10/site-packages/overload.py", line 181, in f
return callable(*usable_args, **_kw)
File "/home/davide/Programs/SNP_Utils-master/Objects/blast.py", line 233, in get_snp_position
return self.get_snp_position(query_snp=lookup.get_code(), expected=adj)
File "/home/davide/.local/lib/python3.10/site-packages/overload.py", line 181, in f
return callable(*usable_args, **_kw)
File "/home/davide/Programs/SNP_Utils-master/Objects/blast.py", line 198, in get_snp_position
raise NoSNPError("Cannot find the SNP in " + self) # Error out to avoid infinite loops
TypeError: can only concatenate str (not "Hsp") to str
A declarative, efficient, and flexible JavaScript library for building user interfaces.
๐ Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
An Open Source Machine Learning Framework for Everyone
The Web framework for perfectionists with deadlines.
A PHP framework for web artisans
Bring data to life with SVG, Canvas and HTML. ๐๐๐
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
Some thing interesting about web. New door for the world.
A server is a program made to process requests and deliver data to clients.
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
Some thing interesting about visualization, use data art
Some thing interesting about game, make everyone happy.
We are working to build community through open source technology. NB: members must have two-factor auth.
Open source projects and samples from Microsoft.
Google โค๏ธ Open Source for everyone.
Alibaba Open Source for everyone
Data-Driven Documents codes.
China tencent open source team.