Need help with blast...

[Guigra]

Hi,
Does anyone know if I need to make a assembly of my alignment before attempting to make a blast?
I'm asking because after I do an alignment on BWA, I created a file in fasta format to make a blast. But the output file is similar to a sequence without anything aligned.

Someone can help me?

[Cofactor Genomics]

Hello Guigra,
Can you share your alignment parameters with the group here... as well as confirm you expect the two to align because you established this with the prior BWA alignments? There are a number of very experienced BLAST users on the board here (as well as on our team) that can help if you provide more information.

My best.

Jarret Glasscock
Cofactor Genomics

http://www.cofactorgenomics.com

[Guigra]

Hello,
After doing the alignment with BWA, I need to pass the generated file to fasta format, so I can make a blast.
But when I do this, apparently the fasta file generated seems only the reads. Lose alignment. You know what's happening?
I use te samtools to pass the file in sam format for bam format, and after I use the blast2fastx to convert the alignment in bam format to fasta format.
I don't know if have some influence, but I'm work with the paired-end.

[westerman]

FastA *is* only the sequence; aka, reads. The alignment and other information will be lost.