-
Notifications
You must be signed in to change notification settings - Fork 16
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Can not produce spliced alignments of whole transcript #47
Comments
Hi @LalitNM, thank you for trying Magic-BLAST and I am sorry you ran into problems. Do you use the same BLAST database in both runs? It is unexpected that you get a continuous alignment in your second run. Are you able to share the cDNA and genomic sequences? Is there a way that I can have a look at them? It would make it easier for me to narrow down the problem. |
Thank you for your response. Yes, I used same BLAST database in both runs. cDNA and genomic sequences are here - cDNA:
Genomic Sequence:
|
Thank you for your sequences. Magic-BLAST produces spliced alignments for all exons with your cDNA sequences. In SAM and BAM files all axons are shown in a single line. You can verify that these are spliced alignments by looking at the CIGAR string (6th column in the SAM file), for example: I am not too familiar with
|
Thank you for helping, though I could not produce desired |
Hi, I tried magic-blast to align a toy cDNA to a toy genome, I could either produce alignment of whole sequence (in blue) or the first exon only (in red).
Commands used to produce results
I want to know is there a way I can produce the expected alignment in a .bed file using magic-blast?
PS: I investigated respective .sam files as well, but could not find the information of all exons and introns there.
The text was updated successfully, but these errors were encountered: