aastat for trset2, trset3 (inprog) /bio/bio-grid/aabugs4/bugs/honbee/trset3,trset2,trset # ** hbee1 ab? or ads? has human/mouse contam, may be longest prots. #... ref,evg ............. #amel_OGSv32 nt=15314; average=1850; median=1498; min,max=1137,23421; nfull=; sum=1850164; gaps=1415,1.4 #evg1hbee.okay nt=68139; average=1953; median=1673; min,max=1329,16948; nfull=957; sum=1953068; gaps=7443,7.4 #evg1hbee.okalt nt=71984; average=2543; median=2225; min,max=1751,16940; nfull=916; sum=2543155; gaps=20541,20.5 #evg2hbee.okay nt=90131; average=2018; median=1716; min,max=1353,16948; nfull=949; sum=2018728; gaps=9302,9.3 #evg2hbee.okalt nt=119046; average=4003; median=3245; min,max=2665,16940; nfull=963; sum=4003910; gaps=10408,10.4 #... trset1 ............. #soaphbee1ab2n nt=727163; average=2371; median=2213; min,max=1863,5649; nfull=865; sum=2371718; gaps=12018,12 #velvhbee1ab2n nt=626033; average=4150; median=3623; min,max=2876,16948; nfull=894; sum=4150416; gaps=3510,3.5 #soaphbee1n2ads1 nt=547362; average=3009; median=2768; min,max=2210,5619; nfull=978; sum=3009598; gaps=30557,30.5 #velvhbee1f2ads1 nt=296142; average=1781; median=1615; min,max=1281,4836; nfull=668; sum=1781437; gaps=32,0 ### soaphbee1nb1 is set1, below #... trset2 ............. #... fatbody fatc has shortish genes #soaphbee1fatc143 nt=147937; average=1134; median=1053; min,max=900,2875; nfull=870; sum=1134936; gaps=7391,7.3 #velvhbee1fatc141 nt=98368; average=1035; median=936; min,max=771,2569; nfull=667; sum=1035840; gaps=3662,3.6 #soaphbee1fatnorm4n nt=215761; average=1335; median=1205; min,max=1049,4305; nfull=882; sum=1335027; gaps=6319,6.3 #velvhbee1fatnorm5n nt=274943; average=1658; median=1436; min,max=1270,4122; nfull=881; sum=1658645; gaps=2504,2.5 #.. #soaphbee1nb1 nt=546976; average=2871; median=2502; min,max=2105,9443; nfull=960; sum=2871260; gaps=21226,21.2 #velvhbee1nb11 nt=71519; average=2402; median=2034; min,max=1618,15633; nfull=906; sum=2402115; gaps=3175,3.1 #velvhbee1brainorm14 nt=786151; average=2742; median=2459; min,max=2057,6683; nfull=813; sum=2742801; gaps=3764,3.7 #soaphbee1brainorm : FAIL, retry #... UofI data set... #velvhbee1ui455a nt=388428; average=5427; median=4615; min,max=3785,16918; nfull=918; sum=5427725; gaps=2852,2.8 #velvhbee1uinorm5n nt=423918; average=4561; median=3926; min,max=3315,16920; nfull=940; sum=4561389; gaps=1978,1.9 #soaphbee1uinorm4n nt=393000; average=3269; median=2928; min,max=2388,8941; nfull=917; sum=3269047; gaps=13900,13.9 #trinhbee1uinorm3 nt=253206; average=2256; median=2061; min,max=1740,5080; nfull=761; sum=2256621; gaps=0,0 #trinhbee1uinorm4 nt=664488; average=2835; median=2602; min,max=2301,5107; nfull=873; sum=2835458; gaps=0,0 #... trset3 ongoing ... #soaphbee1bf3p3 nt=578829; average=3006; median=2646; min,max=2189,11159; nfull=964; sum=3006669; gaps=11913,11.9 #trinhbee1bf3p4 nt=462257; average=2580; median=2407; min,max=1953,13304; nfull=905; sum=2580350; gaps=0,0 #velvhbee1bf3p6 nt=225526; average=3057; median=2588; min,max=2113,10168; nfull=927; sum=3057842; gaps=3126,3.1 #soaphbee1bn5p4 nt=500702; average=2966; median=2618; min,max=2132,11209; nfull=945; sum=2966154; gaps=9702,9.7 #velvhbee1bn5p6 nt=118235; average=2531; median=2200; min,max=1761,5709; nfull=890; sum=2531292; gaps=2902,2.9 #velvhbee1bf6p5a nt=207196; average=3160; median=2775; min,max=2276,9269; nfull=918; sum=3160493; gaps=3827,3.8 #... trset4 .. env main=0 $evigene/scripts/prot/aastat.sh ../trset4/aaset/*.qual #velvhbee1brainorm8n nt=146598; average=2138; median=1816; min,max=1469,12408; nfull=913; sum=2138370; gaps=45,0 #velvhbee1f3brainm8t nt=429606; average=3376; median=2935; min,max=2400,16513; nfull=927; sum=3376774; gaps=0,0 #velvhbee1f3brainm9t nt=529161; average=3538; median=3049; min,max=2497,16513; nfull=936; sum=3538712; gaps=0,0 #velvhbee1n3brainm8n nt=402675; average=3450; median=2937; min,max=2357,16513; nfull=951; sum=3450608; gaps=10,0 #velvhbee1nf1brnorm28n nt=203834; average=2077; median=1810; min,max=1480,9903; nfull=911; sum=2077611; gaps=0,0 #velvhbee1nf1brnorm28t nt=243554; average=2568; median=2224; min,max=1806,9904; nfull=893; sum=2568556; gaps=0,0 #velvhbee1nf1brnorm28u nt=300494; average=2669; median=2350; min,max=1913,9905; nfull=936; sum=2669011; gaps=2,0 #soaphbee1ui454n nt=192639; average=3202; median=2887; min,max=2237,12149; nfull=971; sum=3202302; gaps=15837,15.8 #velvhbee1uinorm8n nt=407666; average=4326; median=3728; min,max=3056,16923; nfull=956; sum=4326620; gaps=252,0.2 #... f3brainorm cmp env main=t $evigene/scripts/prot/aastat.sh aaset/*bee1f3brai*.qual *.aa.qual #trinhbee1f3brainm4 nt=343598; average=1543; median=1395; min,max=1180,5025; nfull=815; sum=1543438; gaps=0,0 #velvhbee1f3brainm8t nt=267449; average=2387; median=2097; min,max=1678,16513; nfull=925; sum=2387992; gaps=19,0 #velvhbee1f3brainm9t nt=325388; average=2427; median=2110; min,max=1702,16513; nfull=937; sum=2427604; gaps=26,0 #soaphbee1f3brainm4n nt=246437; average=1967; median=1764; min,max=1366,5645; nfull=906; sum=1967213; gaps=4391,4.3 #.. #.... several velv sets are hi-kmer only, puke out for k<63 or higher.. not mem problem? killed in oases, velvg w/o err msg #.... need to find sol to velv lo-kmer fails. 1 bug: bad inslen=300, vel-ins-len measured=150 #.... note abysmal trset3/trinity3 stat vs velv,soap same data BUT hi score for trinhbee1bf3p4 #.. redo trinhbee1uinorm4, same data, newer trinvers, using stranded SS flag, about same aa result but 2x trasm #.. trin4 is modest improvement on 1000longaa stats, but at cost of >2x trasm, not sure it will be better after nr filter