in reply to Parsing multispecies database BLAST output

A oneliner is enough
#cat subspecies.txt qspecies1 subjectspecies1 subjsp1.seqA 30 qspecies1 subjectspecies1 subjsp1.seqB 90 qspecies1 subjectspecies1 subjsp1.seqC 100 qspecies1 subjectspecies1 subjsp1.seqD 40 qspecies1 subjectspecies2 subjsp2.seqA 100 qspecies1 subjectspecies2 subjsp2.seqB 12 qspecies2 subjectspecies1 subjsp1.seqD 30 qspecies2 subjectspecies1 subjsp1.seqE 90 qspecies2 subjectspecies2 subjsp2.seqC 10 qspecies2 subjectspecies2 subjsp2.seqD 70 qspecies2 subjectspecies3 subjsp3.seqA 80 # # #perl -MData::Dump -lanE "$hash{$F[0]}{$F[1]}=[$F[2],$F[3]] unless $ha +sh{$F[0]}{$F[1]}[1]>$F[3];END{dd%hash}" subspecies.txt ( "qspecies2", { subjectspecies1 => ["subjsp1.seqE", 90], subjectspecies2 => ["subjsp2.seqD", 70], subjectspecies3 => ["subjsp3.seqA", 80], }, "qspecies1", { subjectspecies1 => ["subjsp1.seqC", 100], subjectspecies2 => ["subjsp2.seqA", 100], }, )

L*

There are no rules, there are no thumbs..
Reinvent the wheel, then learn The Wheel; may be one day you reinvent one of THE WHEELS.