aligmnet of mutli-fasta fie to reference genome
3
0
Entering edit mode
6.1 years ago
anna • 0

Hi, I want to check similarity of the draft genome (250 contigs in total) to the reference genome using BRIG. I made a run but it did not work. On terminal I had some complains in regard to java such as:

java.lang.IllegalArgumentException: Comparison method violates its general contract!
    at java.util.TimSort.mergeHi(TimSort.java:868)
    at java.util.TimSort.mergeAt(TimSort.java:485)
    at java.util.TimSort.mergeCollapse(TimSort.java:408)
    at java.util.TimSort.sort(TimSort.java:214)
    at java.util.TimSort.sort(TimSort.java:173)
    at java.util.Arrays.sort(Arrays.java:659)
    at brig.Main.reverseBlast(Main.java:206)
    at brig.Main.BlastMaster(Main.java:1494)
    at brig.Main.ParseBlast(Main.java:1352)
    at brig.Three2.compute(Three2.java:337)
    at brig.Three2.access$300(Three2.java:33)
    at brig.Three2$4.run(Three2.java:234)

I made a test where I selected one and two contigs and ran that against reference. That worked quite ok. However, when I selected more than 10 contigs I see the same message as in case of 250 contigs. I will appreciate any suggestions.

software error • 1.4k views
ADD COMMENT
0
Entering edit mode

Run it on batches of 10 contigs... 1-10, 11-20, 21-30... see where the error is. Clearly it's one of the contigs, or maybe amount of contigs.

ADD REPLY
1
Entering edit mode
6.1 years ago

This is a programming error. http://stackoverflow.com/questions/11441666

There is nothing you can do unless modifying the java code.

ADD COMMENT
0
Entering edit mode
6.1 years ago

You can use MUMmer tool (http://mummer.sourceforge.net/) for Genome Alignment.

ADD COMMENT
0
Entering edit mode
6.1 years ago

if your data is rna the use bowtie2 or even try shrimp to do it

ADD COMMENT

Login before adding your answer.

Traffic: 1970 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6