ana scripts

Carn_test/run_findDups.sh

@@ -0,0 +1,11 @@
+#!/bin/bash
+#SBATCH -n 2
+#SBATCH -p fast
+#SBATCH -t 02:00:00
+#SBATCH --mem=60GB
+#SBATCH -e Carn_acyltransf.err
+#SBATCH -o Carn_acyltransf.out
+
+
+
+python3.6 ../findDups_countLoss.py ../Carn_acyltransf.faa.nw ../Carn_acyltransf.faa ./23072021_1/ midpoint

README.md

@@ -5,6 +5,8 @@ A simple program to find main OG groups out of a tree
 Example usage: 
 
 ```
-python find_tree_groups.py example_trees/Laminin_G_1.faa.nw.rooted Opisthokonta
+ python3.6 find_groups_get_seqs.py tree_file raw_fasta outdir  midpoint
+
 ```
 
+

example_trees/Antimicrobial12.faa.nw.rooted

@@ -0,0 +1 @@
+(((((144197.ENSSPAP00000002626:1.2182,((32473.ENSXCOP00000025453:1.0829,(8175.ENSSAUP00010037166:1.3838,56716.ENSCGOP00000011067:0.85710):0.48703):0.42382,(8078.ENSFHEP00000012525:0.67069,8128.ENSONIP00000041407:1.1173):0.22263):0.14095):0.21278,((((((8090.ENSORLP00000000172:0.34762,(30732.ENSOMEP00000019115:0.26751,8090.ENSORLP00000038614:0.11898):0.15552):0.10667,(8175.ENSSAUP00010028332:0.089990,(8175.ENSSAUP00010028337:0.0094700,8175.ENSSAUP00010028338:0.13847):0.21781):0.38352):0.28196,(8175.ENSSAUP00010028347:0.40293,64144.ENSATEP00000029251:0.61045):0.071080):0.096180,((8078.ENSFHEP00000012618:0.50727,(28743.ENSCVAP00000018691:0.094110,(8078.ENSFHEP00000029344:0.23069,(48701.ENSPMEP00000007586:0.015770,8083.ENSXMAP00000029961:0.16238):0.053580):0.033650):0.10917):0.25641,(8078.ENSFHEP00000029384:1.2185,((30732.ENSOMEP00000019102:0.51344,(28743.ENSCVAP00000018671:0.15228,(48699.ENSPLAP00000019403:0.0000,48701.ENSPMEP00000025955:0.0000):0.22026):0.069750):0.11290,(181472.ENSSFAP00005001971:0.87261,28743.ENSCVAP00000018657:1.0541):0.53126):0.19300):0.012750):0.12759):0.028910,((173247.ENSENLP00000022339:0.51515,((37003.ENSKMAP00000027613:0.47749,(8078.ENSFHEP00000012531:0.52817,(((48701.ENSPMEP00000005183:0.073440,8083.ENSXMAP00000024724:0.37225):0.32292,((37003.ENSKMAP00000027633:0.13668,((37003.ENSKMAP00000010307:0.0000,37003.ENSKMAP00000027642:0.0000):0.20323,37003.ENSKMAP00000027649:0.057440):0.24675):0.11344,37003.ENSKMAP00000010293:0.51055):0.081440):0.16297,(((((37003.ENSKMAP00000027621:0.10090,52670.A0A2I4CVB8:0.31083):0.12999,8090.ENSORLP00000037391:0.83093):0.17060,(((48699.ENSPLAP00000000696:0.0000,48701.ENSPMEP00000007676:0.0000):0.027710,8081.ENSPREP00000025181:0.056080):0.12614,(32473.ENSXCOP00000025455:0.38396,8083.ENSXMAP00000028942:0.021840):0.19769):0.28245):0.00054000,(28743.ENSCVAP00000018642:0.27423,(8083.ENSXMAP00000026750:0.32297,((48701.ENSPMEP00000005159:0.0075400,(48701.ENSPMEP00000001636:0.00055000,48699.ENSPLAP00000009672:0.023120):0.00054000):0.090410,(8081.ENSPREP00000025164:0.066250,(33528.ENSGAFP00000007805:0.062850,32473.ENSXCOP00000025451:0.093920):0.042720):0.078450):0.055680):0.050130):0.081190):0.040750,(28743.ENSCVAP00000018667:0.41473,(32473.ENSXCOP00000025449:0.81478,8078.ENSFHEP00000029333:0.62234):0.29235):0.43047):0.024700):0.097530):0.021400):0.00055000,((((61819.ENSACIP00000015856:0.00055000,61819.ENSACIP00000015865:0.16567):0.21536,(((144197.ENSSPAP00000002629:0.23396,(((80972.ENSAOCP00000014038:0.037980,80966.ENSAPOP00000032781:0.10515):0.0092100,80972.ENSAOCP00000028780:0.0098600):0.0098100,161767.ENSAPEP00000014562:0.00054000):0.13091):0.18163,(((8187.ENSLCAP00010023299:0.0000,8187.ENSLCAP00010023302:0.0000):0.040910,8187.ENSLCAP00010023306:0.087680):0.27877,(1841481.ENSSLDP00000006899:0.054700,41447.ENSSDUP00000007516:0.080470):0.38233):0.058420):0.13809,(32507.ENSNBRP00000013545:0.26514,(47969.ENSOABP00000047297:0.00055000,8128.ENSONIP00000023045:0.00055000):0.0018200):0.10725):0.034710):0.092670,(((8154.ENSACLP00000040821:0.090880,(32507.ENSNBRP00000013568:0.023900,(47969.ENSOABP00000047295:0.061690,8128.ENSONIP00000023044:0.0077100):0.057320):0.0068100):0.037440,(106582.ENSMZEP00005007278:0.015530,(303518.ENSPNYP00000029227:0.028430,(8153.ENSHBUP00000013818:0.017440,8154.ENSACLP00000040825:0.018210):0.0087600):0.018430):0.014960):0.030100,106582.ENSMZEP00005007282:0.017640):0.082190):0.31688,205130.ENSMAMP00000031607:0.28124):0.044870):0.18442):0.052310,((215358.ENSLCRP00005030485:0.66217,28743.ENSCVAP00000032147:1.2011):0.33269,(303518.ENSPNYP00000029232:0.95356,(((8175.ENSSAUP00010028344:0.54001,8175.ENSSAUP00010028340:0.43207):0.40072,(((409849.ENSPMGP00000022033:0.0013200,(409849.ENSPMGP00000004459:0.36557,409849.ENSPMGP00000004457:0.69364):0.14633):1.1326,33528.ENSGAFP00000007813:0.80514):0.099970,(158456.ENSBSLP00000036699:1.0461,181472.ENSSFAP00005001973:0.95967):0.16568):0.10015):0.099940,(441366.ENSGWIP00000009416:1.4620,(161767.ENSAPEP00000014572:0.094570,80972.ENSAOCP00000028782:0.38227):0.62089):0.49203):0.00054000):0.12862):0.16980):0.064380):0.17481,(((((210632.ENSPRNP00000031955:0.51360,(181472.ENSSFAP00005001968:0.29849,181472.ENSSFAP00005001970:0.21734):0.15841):0.13168,94237.ENSMMOP00000020091:0.43185):0.15107,((375764.ENSSORP00005035214:0.18807,375764.ENSSORP00005035216:0.079580):0.32570,((((80966.ENSAPOP00000012388:0.055440,80966.ENSAPOP00000012399:0.20758):0.30771,(80972.ENSAOCP00000028762:0.30922,(61819.ENSACIP00000010149:0.35357,((47969.ENSOABP00000047299:0.0063500,8128.ENSONIP00000068072:0.0050500):0.17048,((((106582.ENSMZEP00005007314:0.0000,303518.ENSPNYP00000029219:0.0000):0.0000,8153.ENSHBUP00000030535:0.0000):0.0000,8154.ENSACLP00000040828:0.0000):0.055710,32507.ENSNBRP00000013503:0.00054000):0.042660):0.19972):0.22333):0.14522):0.054620,64144.ENSATEP00000029245:0.41053):0.045060,((109280.ENSHCOP00000010297:0.45424,215358.ENSLCRP00005030487:0.31081):0.28349,(41447.ENSSDUP00000007509:0.38788,(205130.ENSMAMP00000031601:0.10716,8187.ENSLCAP00010023292:0.26473):0.13087):0.041590):0.17244):0.034420):0.085410):0.013500,61819.ENSACIP00000010307:0.45649):0.14552,(375764.ENSSORP00005035217:0.30794,28743.ENSCVAP00000018573:1.3718):0.30778):0.37902):0.073600):0.13583,(56716.ENSCGOP00000011071:0.67789,(43700.ENSMALP00000009271:0.59611,43700.ENSMALP00000023037:0.093990):1.0815):0.48601):0.057230,56723.ENSLBEP00000033569:0.64881):0.078604,((((((8175.ENSSAUP00010028361:0.14875,8175.ENSSAUP00010028368:0.18926):0.34711,(205130.ENSMAMP00000000449:0.10514,(205130.ENSMAMP00000031608:0.52148,205130.ENSMAMP00000000456:0.28177):0.29643):0.40352):0.23110,((173247.ENSENLP00000022325:1.0262,((144197.ENSSPAP00000002638:0.94882,(80966.ENSAPOP00000032259:0.55723,(161767.ENSAPEP00000014469:0.36204,80972.ENSAOCP00000014050:0.61180):0.19300):0.24754):0.10588,8187.ENSLCAP00010023327:1.2864):0.12792):0.00055000,(43700.ENSMALP00000023023:0.63116,(43700.ENSMALP00000005643:0.094030,((((8187.ENSLCAP00010023318:0.0000,8187.ENSLCAP00010027336:0.0000):0.53210,(8187.ENSLCAP00010023322:0.0000,8187.ENSLCAP00010027332:0.0000):0.10051):0.11970,(8187.ENSLCAP00010023330:0.20342,(8187.ENSLCAP00010023333:0.18390,(8187.ENSLCAP00010023337:0.0000,8187.ENSLCAP00010027330:0.0000):0.14618):0.10872):0.19056):0.32101,(8187.ENSLCAP00010023342:0.24286,8187.ENSLCAP00010027327:0.12823):0.23705):0.21861):0.22011):0.24937):0.17329):0.062640,(56716.ENSCGOP00000011076:0.96227,(48701.ENSPMEP00000021983:0.65639,(161767.ENSAPEP00000014584:1.1354,210632.ENSPRNP00000031964:0.86650):0.087390):0.18688):0.28383):0.072350,(((375764.ENSSORP00005035212:0.53022,(375764.ENSSORP00005035197:1.0568,(8175.ENSSAUP00010037175:2.0125,215358.ENSLCRP00005030483:0.0013300):1.0558):0.26127):0.21317,(8153.ENSHBUP00000030543:0.013760,(((8128.ENSONIP00000051109:0.032560,(303518.ENSPNYP00000029210:1.1170,47969.ENSOABP00000047300:0.00055000):0.083140):0.10806,106582.ENSMZEP00005007543:0.15834):0.036790,8154.ENSACLP00000040835:0.00055000):0.013290):0.90643):0.20678,(64144.ENSATEP00000029296:0.71574,((181472.ENSSFAP00005001966:0.83726,(375764.ENSSORP00005035195:0.68965,375764.ENSSORP00005035211:0.89480):0.35184):0.13848,((8175.ENSSAUP00010028349:0.40693,8175.ENSSAUP00010028358:0.16510):0.17468,(61819.ENSACIP00000010280:0.74080,(32507.ENSNBRP00000013596:0.68644,((303518.ENSPNYP00000029245:0.020330,(106582.ENSMZEP00005007249:0.0097500,8154.ENSACLP00000040813:0.0096100):0.00055000):0.39055,(47969.ENSOABP00000047291:0.00055000,8128.ENSONIP00000031398:0.00053000):0.57885):0.14709):0.18838):0.72231):0.40613):0.14438):0.29337):0.023160):0.022680,((1841481.ENSSLDP00000006886:0.58658,41447.ENSSDUP00000007537:0.34327):0.33658,((((((((32473.ENSXCOP00000025446:0.17745,64144.ENSATEP00000029274:1.3666):0.052010,(28743.ENSCVAP00000018596:0.22708,8081.ENSPREP00000025125:0.74242):0.16163):0.047090,(((47308.ENSNMLP00000038437:1.7761,8081.ENSPREP00000025091:0.11053):0.34862,(48699.ENSPLAP00000025427:0.016400,48701.ENSPMEP00000005141:0.043890):0.10132):0.12295,8078.ENSFHEP00000029298:0.51595):0.11590):0.21114,((61819.ENSACIP00000010290:0.41069,(47969.ENSOABP00000047293:0.71401,(32507.ENSNBRP00000013579:0.058550,((106582.ENSMZEP00005007252:0.0000,303518.ENSPNYP00000029241:0.0000):0.00054000,(8154.ENSACLP00000040816:0.27175,8153.ENSHBUP00000011585:0.38389):0.17612):0.10669):0.093960):0.026520):0.11302,8128.ENSONIP00000037231:0.66205):0.22613):0.11383,((161767.ENSAPEP00000014529:0.19889,80966.ENSAPOP00000032272:0.16601):0.22272,109280.ENSHCOP00000010301:0.38675):0.00054000):0.17733,43700.ENSMALP00000023052:0.49577):0.14472,(8167.A0A484C9P7:0.15869,8175.ENSSAUP00010028352:0.75690):0.29103):0.12353,(((94237.ENSMMOP00000020090:1.2438,(((48699.ENSPLAP00000025441:0.20096,48701.ENSPMEP00000021986:0.37396):0.15396,(48701.ENSPMEP00000005174:0.00055000,(48699.ENSPLAP00000025431:0.045260,(32473.ENSXCOP00000025448:0.00055000,8083.ENSXMAP00000032811:0.10090):0.16030):0.29311):0.25730):0.37549,(28743.ENSCVAP00000018623:0.72141,8078.ENSFHEP00000029323:0.37728):0.13487):0.14817):0.19018,(((48699.ENSPLAP00000009639:0.33804,48701.ENSPMEP00000009757:0.74885):0.25422,(32473.ENSXCOP00000025445:0.11838,(33528.ENSGAFP00000007780:0.77497,8083.ENSXMAP00000024221:0.22702):0.26547):0.31828):0.14589,8078.ENSFHEP00000029311:0.67403):0.26022):0.076870,(((8081.ENSPREP00000025102:0.59819,((210632.ENSPRNP00000021818:0.13306,210632.ENSPRNP00000031957:0.50195):0.36793,(8187.ENSLCAP00010023310:0.23967,8187.ENSLCAP00010023315:0.066060):0.45908):0.087810):0.12828,(28743.ENSCVAP00000032101:0.64622,43700.ENSMALP00000007881:1.2242):0.16513):0.075220,8078.ENSFHEP00000012484:0.53958):0.093300):0.080440):0.14919):0.12337):0.11629);

example_trees/DUF4776.faa.nw.rooted

@@ -0,0 +1 @@
+(((((66420.A0A194QHM8:0.015280,76193.A0A194RMM9:0.029550):0.25253,(168631.A0A437BLB0:0.20822,7091.H9JAU8:0.24832):0.071560):0.58562,(168631.A0A3S2M4J9:0.65309,(104452.A0A0L7LBR1:0.40950,(7102.A0A2A4J5U4:0.43229,(278856.A0A212FEV8:0.34087,(66420.A0A194QNE9:0.031490,76193.A0A194RNV0:0.072420):0.45585):0.060340):0.051200):0.11566):0.32390):0.90786,(104452.A0A0L7L0U5:1.3075,(((((151549.A0A4C1TVJ8:0.55505,(66420.A0A194PFM2:0.013670,76193.A0A0N1PGP9:0.019170):0.090450):0.058150,278856.A0A212ETH7:0.082880):0.034900,7091.H9J5S3:0.076780):0.045840,168631.A0A437BTV1:0.11879):0.13542,7102.A0A2A4JV95:0.20092):0.65441):0.45104):0.41900,(((((195883.A0A482XA68:0.29217,195883.A0A482XKD5:0.76554):0.46018,121845.A0A3Q0JH94:2.6890):0.24816,(((543379.A0A232F1F6:0.078470,7425.K7JP06:0.22953):0.66926,((2015173.A0A026W4L2:0.20960,((300112.A0A4S2KEP0:0.099930,(64791.A0A151XCR8:0.038060,(456900.A0A195CMP9:0.081760,((471704.A0A151JC81:0.020560,(520822.A0A151I2X0:0.0042400,12957.A0A158NCR6:0.10897):0.057150):0.0052700,34720.A0A195FAG3:0.077930):0.010020):0.013660):0.066180):0.081930,(104421.E2AG90:0.10598,67767.A0A0J7L4C1:0.078060):0.072970):0.044840):0.20300,(178035.A0A154P3K8:0.36141,(597456.A0A0L7QY65:0.21794,(166423.A0A0M9A4D3:0.16807,7460.A0A088AB63:0.15199):0.18830):0.14011):0.19165):0.10817):0.54766,13249.T1HSG9:2.6867):0.15668):0.26650,121224.E0VSI4:1.4812):0.17655,((6973.A0A2P8YAS2:1.4544,(77166.U4TVT2:0.62267,(1661398.A0A482VH85:0.49864,7070.D6WSZ6:0.24266):0.33427):0.95724):0.42657,(((7176.B0WWA5:0.79355,((((69004.A0A182QWD6:0.22871,7168.A0A182NA91:0.13582):0.16548,((199890.A0A182PA65:0.15941,(((30066.A0A182URZ5:0.019150,34690.A0A182TGX1:0.017590):0.0040900,34691.A0A182XI89:0.011000):0.0022200,(1518534.A0A1C7CZW0:0.0056000,7165.F5HK67:0.0079300):0.0014400):0.11148):0.12376,(30069.A0A182YKH7:0.12669,((112268.A0A182VX62:0.10725,139723.A0A182MRP0:0.095880):0.059430,62324.A0A1Y9HF40:0.11282):0.026430):0.11501):0.11391):0.21282,(41427.A0A182ISM6:0.43231,74873.A0A084W1P5:0.36156):0.11852):0.088190,(43151.W5J2R1:0.24973,7167.A0A182FL27:1.0930):0.38053):0.20710):0.83950,(568069.A0A1J1HLM5:0.88155,568069.A0A1J1HMY0:0.85999):0.69517):0.22458,(36166.T1GC78:0.99894,(201502.A0A1A9XPB8:1.0387,((((30019.A0A0M5J1Q1:0.33994,(7222.B4JES9:0.19579,(7232.A0A484BJN0:0.21488,7244.B4M6J2:0.11846):0.040710):0.18757):0.094480,(7260.B4NGS0:0.37044,((7217.B3LZC0:0.34605,(7227.FBpp0084953:0.018590,(7238.B4HZM2:0.022200,7240.B4R1I5:0.00055000):0.014390):0.21689):0.21562,(7266.A0A3B0KMT8:0.10941,(46245.Q29C15:0.0047500,7234.B4GNQ9:0.00055000):0.048020):0.19707):0.088250):0.085300):0.54311,(7260.B4NGS1:0.43782,((30019.A0A0M4EPX2:0.35273,(7232.A0A484BKY8:0.25747,(7222.B4JET0:0.26206,7244.B4M6J3:0.15334):0.049230):0.15111):0.25840,((46245.Q29C16:0.0077300,7234.B4GNQ8:0.0050800):0.32157,(7217.B3LZC1:0.20302,(7227.FBpp0084951:0.014950,(7238.B4HZM3:0.0083900,7240.B4R1I6:0.0054100):0.012630):0.19474):0.16961):0.098720):0.054390):0.56209):0.46398,(((37001.A0A1A9WN19:0.086100,((201502.A0A1A9XNU8:0.00079000,67801.A0A240SXF3:0.0014900):0.017130,(37546.A0A1B0GE25:0.0084400,(7395.A0A1A9VPW0:0.046390,7398.A0A1A9ZFA3:0.013340):0.0037400):0.016240):0.055190):0.24633,(7375.A0A0L0C7G4:0.55683,(35570.A0A1I8NZ08:0.22578,7370.A0A1I8MKX5:0.18374):0.19226):0.11241):0.32727,(((((7266.A0A3B0KS60:0.10574,((7234.B4G4U5:0.010560,(46245.B5DV94:0.0018900,46245.B5DYL9:0.0015500):0.0084400):0.099070,(46245.B5DYM0:0.015740,7234.B4G4U4:0.0085200):0.019260):0.054540):0.14994,(7217.B3LYR0:0.19573,(7227.FBpp0082368:0.018970,(7238.B4HF50:0.0041900,7240.B4R0H2:0.0024800):0.0048700):0.15033):0.21111):0.090760,(7260.B4N9I5:0.25342,7260.B4N9I6:0.48163):0.12846):0.060530,(7260.B4N9I7:0.48844,((((7227.FBpp0082366:0.027950,(7238.B4HF46:0.0081800,7240.B4R0G9:0.0048600):0.015510):0.26698,(7227.FBpp0082365:0.030730,(7238.B4HF45:0.014780,7240.B4R0G8:0.0081700):0.037920):0.33105):0.064400,(7217.B3LYQ6:0.36499,(7217.B3LYQ7:0.024130,7217.B3LYQ8:0.23314):0.26689):0.20041):0.31535,(7266.A0A3B0KHY2:0.16352,(46245.Q298X9:0.0071800,7234.B4G4U6:0.0064900):0.11419):0.29097):0.076340):0.25785):0.065700,(((7222.B4JTC4:0.0088000,7222.B4K022:0.0039500):0.29495,(7232.A0A484AXC7:0.76570,(7232.A0A484AW88:0.28929,7244.B4M046:0.28282):0.029270):0.051930):0.30906,((((7222.B4JTC7:0.011280,7222.B4K018:0.012860):0.39318,(7232.A0A484AUC4:0.59351,7244.A0A0Q9WGW7:0.37450):0.10489):0.079970,(((7222.B4JTC5:0.0012200,7222.B4K021:0.0023700):0.19467,(7244.B4M045:0.19983,(7232.A0A484AV33:0.17355,7232.A0A484AV69:0.36772):0.19398):0.061440):0.071310,(30019.A0A0M5IZF5:0.34146,30019.A0A0M4EH19:0.50642):0.071460):0.037680):0.034070,((7244.A0A0Q9WGU4:0.33004,(7222.B4K432:0.028650,(7222.B4JTC6:0.16389,7222.B4K020:0.0049200):0.050980):0.31896):0.16678,(7244.A0A0Q9WH61:0.47442,((7232.A0A484AUZ2:0.63407,7244.B4M043:0.26869):0.066440,(7232.A0A484AXJ4:0.62612,7244.B4M044:0.29673):0.17337):0.044060):0.045890):0.054920):0.038370):0.063630):0.45106):0.11224):0.12048):0.039080):0.17957):0.12516):0.16109):0.13383);