Question: BLAST geneate Output SAM file , Missing aligned sequences.
1
gravatar for clear.choi
3.3 years ago by
clear.choi30
United States
clear.choi30 wrote:

I have tried to use ncbi blast 2.3.0 version which support output sam file.

I did it like below,

/program/ncbi-blast-2.3.0+/bin/blastn -db /db/class2.db -max_target_seqs 1 -outfmt 17 -query ./S4--S4.fasta -out ./S4--S4.sam

Then I got strange results that first all reference Name and Query Name is wrongly generated.

For example.expected name

queryname = m160213_235908_42273_c100966372550000001823205304301641_s1_p0/32769/0_2213 0.89 28
Reference name = HLA_HLA00688_DRB1_04_03_01_15246_bp

real name,

Query name = BL_ORD_ID:16
Reference Name : Query_1

And also there is no sequences are generated only Cigar , position, Map Quality, I want to put sequences also,

Do you have any idea how to handle blast with output as SAM file?

rna-seq output sam bam blast • 1.9k views
ADD COMMENTlink modified 8 months ago by RamRS21k • written 3.3 years ago by clear.choi30

There is the blast2sam.pl script from samtools. IIRC, this requires the standard pairwise format, which is the least stable, so it may or may not work with the latest blast version.

ADD REPLYlink modified 8 months ago by RamRS21k • written 3.3 years ago by SES8.2k
1
gravatar for Pierre Lindenbaum
3.3 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum120k wrote:

SAM in blast 2.3.0 is still beta

FYI: my student wrote a blast2sam tool

ADD COMMENTlink modified 8 months ago by RamRS21k • written 3.3 years ago by Pierre Lindenbaum120k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1034 users visited in the last hour