Comments (3)
Hello,
Yes, I realized that LTR elements like Dasheng may contain tandem repeats. They are currently filtered. We are still working on a new way to process the tandem repeat information.
LTR_retriever does not count the copy number of flanking sequences, but compares the flanking with internal sequences and discards those with flanking sequences similar to internal sequences. This usually works effectively but may be sabotaged by centromeric repeats where Dasheng is located. It will be helpful to tell more exactly how it is filtered from information in the .defalse file.
Thanks,
Shujun
from ltr_retriever.
The contents of the .defalse file are presented below:
Even when disabling trf to retain Dasheng LTR, it still gets filtered out by the Identifier function in LTR.identifier.pl:
During the self-alignment of Dasheng using blastn, the number of alignment records it generates far exceeds 8:
from ltr_retriever.
Hey sorry for the delay. From the defalse output, the program cannot identify the LTR region, which is because of multiple self-alignment entries reported by blastn as you shown. You can pick out the LTR region by eyes which is the line with 1...442 alignments, but I have not implement a programmatic scheme to find such regions out of multiple self alignments. Can you come up with a way?
Shujun
from ltr_retriever.
Related Issues (20)
- an error in Calculate LAI from EDTA HOT 2
- fa.mod.nmtf.pass.list is empty HOT 2
- LTR_retriever failed to generate a file - data passing in a multithreaded run HOT 2
- Are there any changes in Module 1 between versions 2.9.0 and 2.9.9? HOT 3
- cleanup.pl Bug? HOT 3
- same bug before about LTR.identifier.pl HOT 3
- A reviewer's comment related to the setting of mutaion rate HOT 2
- Dependency checker looking in wrong directory for makeblastdb HOT 1
- Assessing Chrysanthemum lavandulifolium genome assembly HOT 3
- Exploring the transposition profile of specific LTRs HOT 2
- LTR retriever is not compatible with RepeatModeler2 since v2.9.8? HOT 5
- Kimura2Paramter to ks value/Time HOT 4
- Calculated LAI is too large HOT 6
- LTR_retriever produce different result in different data combine ways HOT 1
- low LAI when using ltrharvest result as input HOT 1
- TSDs are not identical HOT 1
- ID Conversion HOT 1
- Is a writable `$genome` necessary? HOT 3
- CP: cannot stat 'DM8.fa.retriever.scn.adj':No such file or directory HOT 2
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from ltr_retriever.