This is a test version of Biostars. For the public version, visit https://www.biostars.org.
overlap sequences in two assemblies

Dear All

I want to find overlapping sequences (overlapping genes ) in two assemblies. For instance;

transcriptome A:

gene1 AGTAGTACTG

transcriptome B:

gene1 AGTAGCTGAT

first I want to compare all genes in two assemblies and find overlapping genes. Later, I want to do a venn diagram to show overlapping genes (number of genes in transcriptome A only, number of genes in transcriptome B only, and shared genes between two assemblies). My point is to find out how similar two transcriptomes.

genome gene alignment next-gen sequencing
• 1,565 views • link
written by Mehmet Mehmet

A homebrew way to do that would be to align the multifasta containing the transcriptome A transcript sequences to the one of transcriptome B, perhaps using BLAT which outputs an easily processable PSL file (although old as format). You can set the -minIdentity of sequence in BLAT, which helps in having a clean output, and you can also set other interesting parameters. Beware, it can be very slow.

• 0 views • link
written by Matteo Schiavinato Matteo Schiavinato

0 answers

No answers yet.

Log in to answer this question.

More posts like this