astrorenales / blasttosam Goto Github PK
View Code? Open in Web Editor NEWTool for converting BLASTN query results to the SAM format
Home Page: https://astror.pavo.uberspace.de
License: MIT License
Tool for converting BLASTN query results to the SAM format
Home Page: https://astror.pavo.uberspace.de
License: MIT License
Hi,
thanks for this nice tool!
I am looking forward to using it.
Could you please specify (maybe also in the README) what BLAST output format is supposed to be the input (ASN.1, XML, etc.)?
Thank you very much.
Best,
Cedric
Hi,
I get the following error message from SAMtools-0.1.19 with the most recent version of BlastToSam when applying it to an NCBI-BLAST+-2.2.30-based alignment of two reconstructed bacterial genomes, i.e., both represented by multiple contigs.
Unfortunately, I can not share the files but the command lines used:
blastn -num_threads 12 -query query.fa -db ./target -outfmt 11 -out ./query-blastn_against-target.asn
blast_formatter -archive ./query-blastn_against-target.asn -out ./query-blastn_against-target.out -outfmt 0
java -jar BlastToSam.jar -i query-blastn_against-target.out -o query-blastn_against-target.out.sam
samtools view -bS query-blastn_against-target.out.sam > query-blastn_against-target.out.bam
Do you have any idea what might be the problem?
[samopen] SAM header is present: 29 sequences.
Line 129, sequence length 129 vs 117 from CIGAR
Parse error at line 129: CIGAR and sequence length are inconsistent
Aborted
Thank you very much for your support.
Best,
Cedric
UPDATE:
I tried SAMtools-1.2 which returns the following
[E::sam_parse1] CIGAR and query sequence are of different length
[W::sam_read1] parse error at line 129
[main_samview] truncated file.
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.