[maker-devel] Odd missing label bug
Fields, Christopher J
cjfields at illinois.edu
Wed Nov 13 20:51:36 MST 2013
I’m seeing an odd bug in GFF output where the label is being sporadically left off some child ‘match_part' features. This is when using maker-2.30p. I found this while splitting out the data by source as individual tracks for IGV viewing; I have several sets of ESTs from individuals of the same species, so I have them labeled (from maker_opts):
est=./ESTs/Trinity_1128_v2_combined_pr2.clean.fasta:strain_1128
protein=./ESTs/xeno/Taeniopygia_guttata.taeGut3.2.4.73.pep.all.fa:TaeGut,./ESTs/xeno/Ficedula_albicollis.FicAlb_1.4.73.pep.all.fa:FicAlb
I planned on separating these into separate tracks, one for each BLASTX, BLASTN, etc.
After a test run on a few scaffolds, a small number of the BLASTX and BLASTN look like this:
KB913263.1 blastx:TaeGut protein_match 374736 401929 330 - . ID=KB913263.1:hit:9:3.10.0.3;Name=ENSTGUP00000000074
KB913263.1 blastx match_part 401756 401929 312 - . ID=KB913263.1:hsp:55:3.10.0.3;Parent=KB913263.1:hit:9:3.10.0.3;Gap=M58;Target=ENSTGUP00000000074 1 58
KB913263.1 blastx match_part 394486 394665 186 - . ID=KB913263.1:hsp:56:3.10.0.3;Parent=KB913263.1:hit:9:3.10.0.3;Gap=M42 D9 M9;Target=ENSTGUP00000000074 56 106
KB913263.1 blastx match_part 392059 392178 204 - . ID=KB913263.1:hsp:57:3.10.0.3;Parent=KB913263.1:hit:9:3.10.0.3;Gap=M40;Target=ENSTGUP00000000074 97 136
KB913263.1 blastx match_part 387479 387547 116 - . ID=KB913263.1:hsp:58:3.10.0.3;Parent=KB913263.1:hit:9:3.10.0.3;Gap=M23;Target=ENSTGUP00000000074 151 173
KB913263.1 blastx match_part 383644 383742 156 - . ID=KB913263.1:hsp:59:3.10.0.3;Parent=KB913263.1:hit:9:3.10.0.3;Gap=M33;Target=ENSTGUP00000000074 172 204
KB913263.1 blastx match_part 382860 383063 242 - . ID=KB913263.1:hsp:60:3.10.0.3;Parent=KB913263.1:hit:9:3.10.0.3;Gap=M68;Target=ENSTGUP00000000074 193 260
KB913263.1 blastx match_part 382538 382648 186 - . ID=KB913263.1:hsp:61:3.10.0.3;Parent=KB913263.1:hit:9:3.10.0.3;Gap=M37;Target=ENSTGUP00000000074 247 283
KB913263.1 blastx match_part 382072 382236 268 - . ID=KB913263.1:hsp:62:3.10.0.3;Parent=KB913263.1:hit:9:3.10.0.3;Gap=M55;Target=ENSTGUP00000000074 282 336
KB913263.1 blastx match_part 379270 379458 330 - . ID=KB913263.1:hsp:63:3.10.0.3;Parent=KB913263.1:hit:9:3.10.0.3;Gap=M63;Target=ENSTGUP00000000074 336 398
KB913263.1 blastx match_part 374736 374900 147 - . ID=KB913263.1:hsp:64:3.10.0.3;Parent=KB913263.1:hit:9:3.10.0.3;Gap=M55;Target=ENSTGUP00000000074 4 58
Note the missing label in the source column for the ‘match_part’ types. Strangely, only a small # have this problem, but the exact same ones appear to be consistently popping up on repeated runs (this is from a test prior to a full launch using openmpi). It does seem like only child ‘match_part’ appears affected.
Any ideas?
chris
More information about the maker-devel
mailing list