i have a transcriptome dataset of an insect, and it's must longest transcript is 24081b, the shortest one is 101b, but its not a model species. so i want to align the dataset with the genome of Tribolium castaneum, because they belong to the same family.i am hoping by aligning my specimen' transcriptome with the genome of Tribolium castaneum to find out how many transcripts could be aligned and what r those genes function in Tribolium castaneum ?
so, by using what kind of tool could i get what i want?
looking forward tips
sincerely.
3 answers
I would use Tophat for starters. The manual doesn't specify a top end for read length, which makes sense because Tophat chops up each read into segments and tries to map those segments early in the process.
You might consider gmap, similar to blat, but much faster.
Log in to answer this question.
-
bacterial genome assembly with nano pore sequence data
written by rthapa •
-
Difference between NCBI mRNA sequences and RNA-seq assembled sequences
written by Kurban
-
MitoSeek to extract animal mtDNA
written by mpjuers •
-
how to do GOSlim-Generic
written by Kurban
-
how to identify same protien sequences in one fasta file
written by Kurban
-
can I get a heat map from two groups of transcriptome data?
written by Kurban
-
how do I check the sequencing reads
written by Kurban
-
transcription factor database
written by Kurban
-
how to annotate transcriptome
written by Kurban
-
align transcriptome to genome
written by Kurban
By "transcriptome", do you mean NGS data (short reads) or do you mean FASTA sequences of your transcriptome?
Maybe I was not clear before, sorry. I meant FASTA sequences of transcriptome, not NGS data.
Thanks for the clarification. In that case, I would probably not recommend an aligner that was written for short reads (such as STAR or TopHat). If someone has evidence that these aligners do as well as more general-purpose aligners like blat and gmap, I'd love to hear it (because I don't have any data either way).
i got what u meant. recently I have read a paper"STAR: ultrafast universal RNA-seq aligner", and it said STAR can align a transcript to a genome when it's length reached to around 6k, but after u recommended BLAT, I found a paper "Using BLAT to Find Sequence Similarity in Closely Related Genomes", by reading it I guess that might be the tool. thank you.