oushujun / LTR_retriever

LTR_retriever is a highly accurate and sensitive program for identification of LTR retrotransposons; The LTR Assembly Index (LAI) is also included in this package.
https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5813529/
GNU General Public License v3.0
177 stars 40 forks source link

Some errors about :LAI is not applicable on the current genome assembly #57

Closed 123guoxiaoyu closed 4 years ago

123guoxiaoyu commented 4 years ago

Hello, I am running LTR retriever, and I am coming across these errors: 1.when i run /LTR_retriever/LTR_retriever,It displays the following error, but it can continue to run. Fri Oct 25 06:19:16 Local time zone must be set--see zic manual page 2019 1 truncated LTR-RTs found Warning: LOC list /dellfsqd2/ST_OCEAN/USER/zhangrui/test/guo/Flow_test/rawdata/smdlt_test.part-1.fa.ltrTE.veryfalse is empty. 2.finally,it displays the following error,I want to know whether the reason for this error is that there are too few LTRs in my input file, or there is a problem with my operation. Fri Oct 25 06:31:10 Local time zone must be set--see zic manual page 2019 The identity of LTR sequences: 88.7317261904762%

                            【Warning】 The identity drops below the safe limit. Instead, identity of 92% will be used for LAI adjustment.

Fri Oct 25 06:31:10 Local time zone must be set--see zic manual page 2019 Calculate LAI:

                            【Error】Intact LTR-RT content (0.02%) is too low for accurate LAI calculation (min 0.1% required)
                            【Error】 Total LTR sequence content (0.3%) is too low for accurate LAI calculation (min 5% required)

                            Sorry, LAI is not applicable on the current genome assembly.

Fri Oct 25 06:31:13 Local time zone must be set--see zic manual page 2019 All analyses were finished! 3.(1)Here are example commands that lead to the error: ./ltr_finder -D 15000 -d 1000 -L 7000 -l 100 -p 20 -C -M 0.8 ./rawdata/smdlt_test.fa > ./rawdata/smdlt.finder.scn (2)./LTR_retriever -threads 4 -genome ./rawdata/1.fa -infinder ./rawdata/1.finder.scn my input file is a whole genome.fa Any ideas as to what the issue is? Thank you!

oushujun commented 4 years ago

Hello,

If you have not, please update LTR_retriever to v2.7 and try again.

Please use the Arabidopsis genome to test on the pipeline and make sure it's correctly installed.

If you are testing the pipeline with a small subset of your data, try to run it on the full set. It's pretty quick. The wheat genome took 8 days.

If it works properly, and still says limited intact LTR and total LTR, that's what it literally means. Not all species has abundant and avtive LTR retrotransposons.

Hope these helps!

Best, Shujun

On Fri, Oct 25, 2019, 2:38 AM 123guoxiaoyu notifications@github.com wrote:

Hello, I am running LTR retriever, and I am coming across these errors: 1.when i run /LTR_retriever/LTR_retriever,It displays the following error, but it can continue to run. Fri Oct 25 06:19:16 Local time zone must be set--see zic manual page 2019 1 truncated LTR-RTs found Warning: LOC list /dellfsqd2/ST_OCEAN/USER/zhangrui/test/guo/Flow_test/rawdata/smdlt_test.part-1.fa.ltrTE.veryfalse is empty. 2.finally,it displays the following error,I want to know whether the reason for this error is that there are too few LTRs in my input file, or there is a problem with my operation. Fri Oct 25 06:31:10 Local time zone must be set--see zic manual page 2019 The identity of LTR sequences: 88.7317261904762%

                        【Warning】 The identity drops below the safe limit. Instead, identity of 92% will be used for LAI adjustment.

Fri Oct 25 06:31:10 Local time zone must be set--see zic manual page 2019 Calculate LAI:

                        【Error】Intact LTR-RT content (0.02%) is too low for accurate LAI calculation (min 0.1% required)

                        【Error】 Total LTR sequence content (0.3%) is too low for accurate LAI calculation (min 5% required)

                        Sorry, LAI is not applicable on the current genome assembly.

Fri Oct 25 06:31:13 Local time zone must be set--see zic manual page 2019 All analyses were finished! 3.(1)Here are example commands that lead to the error: ./ltr_finder -D 15000 -d 1000 -L 7000 -l 100 -p 20 -C -M 0.8 ./rawdata/smdlt_test.fa > ./rawdata/smdlt.finder.scn (2)./LTR_retriever -threads 4 -genome ./rawdata/1.fa -infinder ./rawdata/1.finder.scn my input file is a whole genome.fa Any ideas as to what the issue is? Thank you!

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/oushujun/LTR_retriever/issues/57?email_source=notifications&email_token=ABNX4NGQKWVSEYR62DZX3XTQQKO5PA5CNFSM4JE7W22KYY3PNVWWK3TUL52HS4DFUVEXG43VMWVGG33NNVSW45C7NFSM4HUJ5JKQ, or unsubscribe https://github.com/notifications/unsubscribe-auth/ABNX4NGZWGR6MGI7IZLNHP3QQKO5PANCNFSM4JE7W22A .

123guoxiaoyu commented 4 years ago

Hi Shujun, Thank you for your reply! After I tried the way you said, I was very happy and solved the problem. Thanks again! Regreet xiaoyu