I converted my BlasyN to SAM using:
blast2sam.pl <blastN file> > <output file>
I then tried to generate a "bam" file using:
samtools view -bT <reference.fasta> <sam file> -o <bam output>
I got the following output:
[sam_header_read2] 0 sequences loaded.
[samopen] no @SQ lines in the header.
[sam_read1] missing header? Abort!
Any help?
Thanks in advance
blast2sam.pl <blastN file> > <output file>
I then tried to generate a "bam" file using:
samtools view -bT <reference.fasta> <sam file> -o <bam output>
I got the following output:
[sam_header_read2] 0 sequences loaded.
[samopen] no @SQ lines in the header.
[sam_read1] missing header? Abort!
Any help?
Thanks in advance
Comment