Hi everyone!
I'm trying to efficiently align a fasta file with a human sequence to panTro6 sequence but I'm not sure on how to avoid the gaps between both species. Are there tools that take these gaps already into account? Thanks in advance for the help!
1 answer
And another question @JC : I have done the lift over and I'm having a look at the sequences with UGene (similar to Bioedit). I see now that a big part of panTro6 sequence is displaced just by 1 or 2 nucleotides in some regions. My next step will be to calculate the divergence between both species, so I'm wondering whether this will affect a lot my results... should I eliminate all those "unmatching" regions from both species?? Or should I add/remove this nucleotides to make the sequences match?
Thanks a lot again!
Log in to answer this question.
-
Multi Alignment - How to deal with uncertainty?
written by andrebolerbarros •
-
Sorting and writing multifasta entries to new fasta files
written by lachiemck •
-
Issues with BedToIntervalList Picard tool
written by Elisa •
-
mitochondrial DNA in vcf for haplogrep2
written by biosol
-
mean Tajima's D value
written by biosol
-
Picard tools duplicate removal
written by blur
-
crossover operator for evolutionary algorithm
written by Moses
-
Best way to compare rate of amino acid change between proteins?
written by Rubal
-
What tools to perform alignment without adding gaps in?
written by biolab
-
Parsing Protein Trees to determine orthologs and paralogs
written by lchau91 •
what does it mean ? if the sequences are large, there will be gaps. you cannot avoid it.
I would like to eliminate from my human fasta those positions that don't have a match with the panTro fasta.