Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

About RepeatMasker's search engine #494

Closed
ken0-1n opened this issue May 14, 2021 · 1 comment
Closed

About RepeatMasker's search engine #494

ken0-1n opened this issue May 14, 2021 · 1 comment

Comments

@ken0-1n
Copy link

ken0-1n commented May 14, 2021

Hi, I'm trying to use GRIDSS_annotate_vcf_repeatmasker.sh. I'm thinking of using HMMER as the search engine for RepeatMasker.
When I run RepeatMasker with HMMER as the search engine, I get bit score instead of SW in rm.fa.out.

Does GRIDSS support the output with HMMER as the search engine?

@ken0-1n
Copy link
Author

ken0-1n commented May 14, 2021

I ran GRIDSS_annotate_vcf_repeatmasker.sh (2.11.1) and got an error.

Exception in thread "main" htsjdk.tribble.TribbleException$MalformedFeatureFile: Error parsing LineIteratorImpl(SynchronousLineReader) at the first record, for input source: ./gridss.vcf.fa.out
        at htsjdk.tribble.TribbleIndexedFeatureReader$WFIterator.readNextRecord(TribbleIndexedFeatureReader.java:396)
        at htsjdk.tribble.TribbleIndexedFeatureReader$WFIterator.<init>(TribbleIndexedFeatureReader.java:342)
        at htsjdk.tribble.TribbleIndexedFeatureReader.iterator(TribbleIndexedFeatureReader.java:309)
        at gridss.repeatmasker.AnnotateVariantsRepeatMasker.doWork(AnnotateVariantsRepeatMasker.java:64)
        at picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:305)
        at gridss.repeatmasker.AnnotateVariantsRepeatMasker.main(AnnotateVariantsRepeatMasker.java:100)
Caused by: java.lang.NumberFormatException: For input string: "bit"
        at java.base/java.lang.NumberFormatException.forInputString(NumberFormatException.java:65)
        at java.base/java.lang.Integer.parseInt(Integer.java:652)
        at java.base/java.lang.Integer.parseInt(Integer.java:770)
        at au.edu.wehi.idsv.repeatmasker.RepeatMaskerCodec.decodeSummary(RepeatMaskerCodec.java:94)
        at au.edu.wehi.idsv.repeatmasker.RepeatMaskerCodec.decode(RepeatMaskerCodec.java:52)
        at au.edu.wehi.idsv.repeatmasker.RepeatMaskerCodec.decode(RepeatMaskerCodec.java:17)
        at htsjdk.tribble.TribbleIndexedFeatureReader$WFIterator.readNextRecord(TribbleIndexedFeatureReader.java:375)
        ... 5 more
"java -Xmx2g $jvm_args -cp $GRIDSS_JAR gridss.repeatmasker.AnnotateVariantsRepeatMasker --INPUT "$input_vcf" --OUTPUT "$output" --REPEAT_MASKER "$rmfile"" command completed with exit code 1.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants