Entering edit mode
6.5 years ago
Pranavathiyani G
▴
330
I have a list of genes (sequences) from bacteria and also have the whole genome sequence, how can i know the genes are encoded in positive strand or negative strand? Is there any possible way to find out?
If it's not available, is there any way to find it?
Perform a re-alignment to the genome using your sequences. Just use something simple like HISAT2 followed by StringTie, and then you can infer the strand to which the reads have aligned. I have done this before with bacteria and it works quite well.
Align (blast) gene sequences to reference sequence. If gene sequences align to reference sequence, then it is the on the same direction of your reference. If it aligns to complementary strand of reference sequence, then it is on the other strand.
Yeah a simple blast can do it.