moritzbuck / mOTUlizer

Utility to analyse a group of closely related MAgs/Genomes/bins/SUBs of more or less dubious origin
GNU General Public License v3.0
12 stars 4 forks source link

Your analysis on Prochlorococcus clade A #2

Open ggautreau opened 3 years ago

ggautreau commented 3 years ago

Hi @moritzbuck,

Thank you for your soft mOTUpan, looks interesting.

I got this issue while I tried to reproduce your analysis on Prochlorococcus clade A.

Do you have any idea about the reason for this error?

mOTUpan.py --faas annot/*/*.faa -o output.tsv
387  genomes in your clade
Creating mOTU for mOTUpan
concatenating all faas
running mmseqs easy-cluster with params --min-seq-id 0.8 --cov-mode 0 -c 0.8 in /tmp/CJGIVWYA_mmseqsCluster_3twi_jpa
For 425577 CDSes we got  25264  gene-clusters
iteration 1 :  0 sum_abs_LLHR: 11857620.542042103
iteration 2 :  0 sum_abs_LLHR: 1378284.0516655238

Your CJGIVWYA_-run for 387 genomes (with mean initial completeness 95.00) resulted
in a core of 0 traits with a total sum of loglikelihood-ratios 1378284.05 and a corrected 
mean completness of 0.01, resulting to a estimated mean traits per genome count of 10968087.86

Traceback (most recent call last):
  File "[HIDDEN_PATH]/mOTUpan.py", line 163, in <module>
    motupan(args)
  File "mOTUpan.py", line 123, in motupan
    print(motu.pretty_pan_table(), file = out_handle)
  File "[HIDDEN_PATH]//lib/python3.7/site-packages/mOTUlizer/classes/mOTU.py", line 313, in pretty_pan_table
    out_dict[k]['type'] = 'core' if k in stats['core'] else 'accessory'
TypeError: argument of type 'NoneType' is not iterable

My .faa files are provided by prokka.

$ls annot/*/*.faa
annot/GCA_000634135/GCA_000634135.faa  annot/GCA_003278925/GCA_003278925.faa  annot/GCA_902571575/GCA_902571575.faa
annot/GCA_000634155/GCA_000634155.faa  annot/GCA_003278965/GCA_003278965.faa  annot/GCA_902572075/GCA_902572075.faa
annot/GCA_000634175/GCA_000634175.faa  annot/GCA_003279085/GCA_003279085.faa  annot/GCA_902572235/GCA_902572235.faa
annot/GCA_000634215/GCA_000634215.faa  annot/GCA_003280565/GCA_003280565.faa  annot/GCA_902572285/GCA_902572285.faa
annot/GCA_000634235/GCA_000634235.faa  annot/GCA_003280615/GCA_003280615.faa  annot/GCA_902572355/GCA_902572355.faa
annot/GCA_000634255/GCA_000634255.faa  annot/GCA_003280645/GCA_003280645.faa  annot/GCA_902573135/GCA_902573135.faa
annot/GCA_000634275/GCA_000634275.faa  annot/GCA_003280655/GCA_003280655.faa  annot/GCA_902573685/GCA_902573685.faa
annot/GCA_000634295/GCA_000634295.faa  annot/GCA_003280685/GCA_003280685.faa  annot/GCA_902573925/GCA_902573925.faa
annot/GCA_000634315/GCA_000634315.faa  annot/GCA_003280735/GCA_003280735.faa  annot/GCA_902574315/GCA_902574315.faa
annot/GCA_000634335/GCA_000634335.faa  annot/GCA_003280775/GCA_003280775.faa  annot/GCA_902574545/GCA_902574545.faa
annot/GCA_000634355/GCA_000634355.faa  annot/GCA_003280805/GCA_003280805.faa  annot/GCA_902574585/GCA_902574585.faa
annot/GCA_000634375/GCA_000634375.faa  annot/GCA_003280815/GCA_003280815.faa  annot/GCA_902576595/GCA_902576595.faa
annot/GCA_000634415/GCA_000634415.faa  annot/GCA_003280835/GCA_003280835.faa  annot/GCA_902577315/GCA_902577315.faa
annot/GCA_000634435/GCA_000634435.faa  annot/GCA_003280895/GCA_003280895.faa  annot/GCA_902577325/GCA_902577325.faa
annot/GCA_000634455/GCA_000634455.faa  annot/GCA_003280935/GCA_003280935.faa  annot/GCA_902577725/GCA_902577725.faa
annot/GCA_000634475/GCA_000634475.faa  annot/GCA_003280995/GCA_003280995.faa  annot/GCA_902577735/GCA_902577735.faa
annot/GCA_000634495/GCA_000634495.faa  annot/GCA_003281145/GCA_003281145.faa  annot/GCA_902578065/GCA_902578065.faa
annot/GCA_000634535/GCA_000634535.faa  annot/GCA_003281205/GCA_003281205.faa  annot/GCA_902578075/GCA_902578075.faa
annot/GCA_000634555/GCA_000634555.faa  annot/GCA_003281225/GCA_003281225.faa  annot/GCA_902578835/GCA_902578835.faa
annot/GCA_000634575/GCA_000634575.faa  annot/GCA_003281235/GCA_003281235.faa  annot/GCA_902578875/GCA_902578875.faa
annot/GCA_000634595/GCA_000634595.faa  annot/GCA_003281565/GCA_003281565.faa  annot/GCA_902578895/GCA_902578895.faa
annot/GCA_000634635/GCA_000634635.faa  annot/GCA_902518225/GCA_902518225.faa  annot/GCA_902579805/GCA_902579805.faa
annot/GCA_000634655/GCA_000634655.faa  annot/GCA_902519155/GCA_902519155.faa  annot/GCA_902580075/GCA_902580075.faa
annot/GCA_000634675/GCA_000634675.faa  annot/GCA_902519425/GCA_902519425.faa  annot/GCA_902580965/GCA_902580965.faa
annot/GCA_000634695/GCA_000634695.faa  annot/GCA_902519535/GCA_902519535.faa  annot/GCA_902581105/GCA_902581105.faa
annot/GCA_000634715/GCA_000634715.faa  annot/GCA_902519565/GCA_902519565.faa  annot/GCA_902581135/GCA_902581135.faa
annot/GCA_000634735/GCA_000634735.faa  annot/GCA_902519705/GCA_902519705.faa  annot/GCA_902581165/GCA_902581165.faa
annot/GCA_000634755/GCA_000634755.faa  annot/GCA_902519915/GCA_902519915.faa  annot/GCA_902581705/GCA_902581705.faa
annot/GCA_000634775/GCA_000634775.faa  annot/GCA_902520005/GCA_902520005.faa  annot/GCA_902581725/GCA_902581725.faa
annot/GCA_000634795/GCA_000634795.faa  annot/GCA_902520305/GCA_902520305.faa  annot/GCA_902582185/GCA_902582185.faa
annot/GCA_000634815/GCA_000634815.faa  annot/GCA_902520615/GCA_902520615.faa  annot/GCA_902582195/GCA_902582195.faa
annot/GCA_000634835/GCA_000634835.faa  annot/GCA_902521105/GCA_902521105.faa  annot/GCA_902582565/GCA_902582565.faa
annot/GCA_000634855/GCA_000634855.faa  annot/GCA_902521185/GCA_902521185.faa  annot/GCA_902582995/GCA_902582995.faa
annot/GCA_000634875/GCA_000634875.faa  annot/GCA_902521235/GCA_902521235.faa  annot/GCA_902583165/GCA_902583165.faa
annot/GCA_000634895/GCA_000634895.faa  annot/GCA_902521405/GCA_902521405.faa  annot/GCA_902583795/GCA_902583795.faa
annot/GCA_000634915/GCA_000634915.faa  annot/GCA_902522105/GCA_902522105.faa  annot/GCA_902583915/GCA_902583915.faa
annot/GCA_000634935/GCA_000634935.faa  annot/GCA_902522905/GCA_902522905.faa  annot/GCA_902583935/GCA_902583935.faa
annot/GCA_000634955/GCA_000634955.faa  annot/GCA_902523415/GCA_902523415.faa  annot/GCA_902584605/GCA_902584605.faa
annot/GCA_000634975/GCA_000634975.faa  annot/GCA_902524975/GCA_902524975.faa  annot/GCA_902584875/GCA_902584875.faa
annot/GCA_000634995/GCA_000634995.faa  annot/GCA_902525715/GCA_902525715.faa  annot/GCA_902585675/GCA_902585675.faa
annot/GCA_000635015/GCA_000635015.faa  annot/GCA_902525905/GCA_902525905.faa  annot/GCA_902585975/GCA_902585975.faa
annot/GCA_000635035/GCA_000635035.faa  annot/GCA_902526035/GCA_902526035.faa  annot/GCA_902586365/GCA_902586365.faa
annot/GCA_000635055/GCA_000635055.faa  annot/GCA_902526575/GCA_902526575.faa  annot/GCA_902586725/GCA_902586725.faa
annot/GCA_000635075/GCA_000635075.faa  annot/GCA_902526685/GCA_902526685.faa  annot/GCA_902588025/GCA_902588025.faa
annot/GCA_000635095/GCA_000635095.faa  annot/GCA_902527465/GCA_902527465.faa  annot/GCA_902588035/GCA_902588035.faa
annot/GCA_000635115/GCA_000635115.faa  annot/GCA_902528015/GCA_902528015.faa  annot/GCA_902588405/GCA_902588405.faa
annot/GCA_000635135/GCA_000635135.faa  annot/GCA_902528325/GCA_902528325.faa  annot/GCA_902588875/GCA_902588875.faa
annot/GCA_000635155/GCA_000635155.faa  annot/GCA_902528615/GCA_902528615.faa  annot/GCA_902588905/GCA_902588905.faa
annot/GCA_000635175/GCA_000635175.faa  annot/GCA_902528695/GCA_902528695.faa  annot/GCA_902589645/GCA_902589645.faa
annot/GCA_000635195/GCA_000635195.faa  annot/GCA_902529015/GCA_902529015.faa  annot/GCA_902589905/GCA_902589905.faa
annot/GCA_000635215/GCA_000635215.faa  annot/GCA_902529355/GCA_902529355.faa  annot/GCA_902591985/GCA_902591985.faa
annot/GCA_000635235/GCA_000635235.faa  annot/GCA_902530275/GCA_902530275.faa  annot/GCA_902592175/GCA_902592175.faa
annot/GCA_000635255/GCA_000635255.faa  annot/GCA_902530375/GCA_902530375.faa  annot/GCA_902592295/GCA_902592295.faa
annot/GCA_000635275/GCA_000635275.faa  annot/GCA_902530845/GCA_902530845.faa  annot/GCA_902592885/GCA_902592885.faa
annot/GCA_000635295/GCA_000635295.faa  annot/GCA_902531115/GCA_902531115.faa  annot/GCA_902593615/GCA_902593615.faa
annot/GCA_000635315/GCA_000635315.faa  annot/GCA_902531355/GCA_902531355.faa  annot/GCA_902593725/GCA_902593725.faa
annot/GCA_000635335/GCA_000635335.faa  annot/GCA_902531745/GCA_902531745.faa  annot/GCA_902593785/GCA_902593785.faa
annot/GCA_000635355/GCA_000635355.faa  annot/GCA_902532075/GCA_902532075.faa  annot/GCA_902594065/GCA_902594065.faa
annot/GCA_000635375/GCA_000635375.faa  annot/GCA_902532605/GCA_902532605.faa  annot/GCA_902595095/GCA_902595095.faa
annot/GCA_000635395/GCA_000635395.faa  annot/GCA_902533015/GCA_902533015.faa  annot/GCA_902595615/GCA_902595615.faa
annot/GCA_000635415/GCA_000635415.faa  annot/GCA_902534475/GCA_902534475.faa  annot/GCA_902595995/GCA_902595995.faa
annot/GCA_000635435/GCA_000635435.faa  annot/GCA_902534515/GCA_902534515.faa  annot/GCA_902596865/GCA_902596865.faa
annot/GCA_000635455/GCA_000635455.faa  annot/GCA_902534615/GCA_902534615.faa  annot/GCA_902596885/GCA_902596885.faa
annot/GCA_000635475/GCA_000635475.faa  annot/GCA_902534835/GCA_902534835.faa  annot/GCA_902596945/GCA_902596945.faa
annot/GCA_000635495/GCA_000635495.faa  annot/GCA_902534925/GCA_902534925.faa  annot/GCA_902597005/GCA_902597005.faa
annot/GCA_000635515/GCA_000635515.faa  annot/GCA_902535525/GCA_902535525.faa  annot/GCA_902597615/GCA_902597615.faa
annot/GCA_000635555/GCA_000635555.faa  annot/GCA_902535845/GCA_902535845.faa  annot/GCA_902597665/GCA_902597665.faa
annot/GCA_000635595/GCA_000635595.faa  annot/GCA_902536165/GCA_902536165.faa  annot/GCA_902597705/GCA_902597705.faa
annot/GCA_000635615/GCA_000635615.faa  annot/GCA_902536235/GCA_902536235.faa  annot/GCA_902598035/GCA_902598035.faa
annot/GCA_000635635/GCA_000635635.faa  annot/GCA_902536265/GCA_902536265.faa  annot/GCA_902598635/GCA_902598635.faa
annot/GCA_000635655/GCA_000635655.faa  annot/GCA_902536415/GCA_902536415.faa  annot/GCA_902599035/GCA_902599035.faa
annot/GCA_000635675/GCA_000635675.faa  annot/GCA_902537015/GCA_902537015.faa  annot/GCA_902599095/GCA_902599095.faa
annot/GCA_000635695/GCA_000635695.faa  annot/GCA_902537125/GCA_902537125.faa  annot/GCA_902599895/GCA_902599895.faa
annot/GCA_000635715/GCA_000635715.faa  annot/GCA_902537375/GCA_902537375.faa  annot/GCA_902600275/GCA_902600275.faa
annot/GCA_000635735/GCA_000635735.faa  annot/GCA_902537495/GCA_902537495.faa  annot/GCA_902600335/GCA_902600335.faa
annot/GCA_000635755/GCA_000635755.faa  annot/GCA_902537585/GCA_902537585.faa  annot/GCA_902600805/GCA_902600805.faa
annot/GCA_000635775/GCA_000635775.faa  annot/GCA_902537825/GCA_902537825.faa  annot/GCA_902603085/GCA_902603085.faa
annot/GCA_000635795/GCA_000635795.faa  annot/GCA_902537995/GCA_902537995.faa  annot/GCA_902606595/GCA_902606595.faa
annot/GCA_000635815/GCA_000635815.faa  annot/GCA_902538325/GCA_902538325.faa  annot/GCA_902610025/GCA_902610025.faa
annot/GCA_000635835/GCA_000635835.faa  annot/GCA_902538675/GCA_902538675.faa  annot/GCA_902610505/GCA_902610505.faa
annot/GCA_000635855/GCA_000635855.faa  annot/GCA_902539255/GCA_902539255.faa  annot/GCA_902611725/GCA_902611725.faa
annot/GCA_000635875/GCA_000635875.faa  annot/GCA_902539785/GCA_902539785.faa  annot/GCA_902612025/GCA_902612025.faa
annot/GCA_000760035/GCA_000760035.faa  annot/GCA_902540005/GCA_902540005.faa  annot/GCA_902612425/GCA_902612425.faa
annot/GCA_000760115/GCA_000760115.faa  annot/GCA_902540205/GCA_902540205.faa  annot/GCA_902613285/GCA_902613285.faa
annot/GCA_001180245/GCA_001180245.faa  annot/GCA_902540935/GCA_902540935.faa  annot/GCA_902613485/GCA_902613485.faa
annot/GCA_001180285/GCA_001180285.faa  annot/GCA_902541385/GCA_902541385.faa  annot/GCA_902614045/GCA_902614045.faa
annot/GCA_001180325/GCA_001180325.faa  annot/GCA_902541565/GCA_902541565.faa  annot/GCA_902614095/GCA_902614095.faa
annot/GCA_001627435/GCA_001627435.faa  annot/GCA_902541645/GCA_902541645.faa  annot/GCA_902614255/GCA_902614255.faa
annot/GCA_001627445/GCA_001627445.faa  annot/GCA_902541735/GCA_902541735.faa  annot/GCA_902614505/GCA_902614505.faa
annot/GCA_001627475/GCA_001627475.faa  annot/GCA_902542205/GCA_902542205.faa  annot/GCA_902614535/GCA_902614535.faa
annot/GCA_001627525/GCA_001627525.faa  annot/GCA_902542885/GCA_902542885.faa  annot/GCA_902614565/GCA_902614565.faa
annot/GCA_001989415/GCA_001989415.faa  annot/GCA_902544545/GCA_902544545.faa  annot/GCA_902615195/GCA_902615195.faa
annot/GCA_001989435/GCA_001989435.faa  annot/GCA_902544775/GCA_902544775.faa  annot/GCA_902615475/GCA_902615475.faa
annot/GCA_001989455/GCA_001989455.faa  annot/GCA_902545685/GCA_902545685.faa  annot/GCA_902615955/GCA_902615955.faa
annot/GCA_002025865/GCA_002025865.faa  annot/GCA_902545745/GCA_902545745.faa  annot/GCA_902616045/GCA_902616045.faa
annot/GCA_002025905/GCA_002025905.faa  annot/GCA_902546505/GCA_902546505.faa  annot/GCA_902616285/GCA_902616285.faa
annot/GCA_002026005/GCA_002026005.faa  annot/GCA_902546555/GCA_902546555.faa  annot/GCA_902616505/GCA_902616505.faa
annot/GCA_002026085/GCA_002026085.faa  annot/GCA_902546595/GCA_902546595.faa  annot/GCA_902616515/GCA_902616515.faa
annot/GCA_002026095/GCA_002026095.faa  annot/GCA_902547805/GCA_902547805.faa  annot/GCA_902616815/GCA_902616815.faa
annot/GCA_002026175/GCA_002026175.faa  annot/GCA_902548105/GCA_902548105.faa  annot/GCA_902617165/GCA_902617165.faa
annot/GCA_002026185/GCA_002026185.faa  annot/GCA_902550355/GCA_902550355.faa  annot/GCA_902617455/GCA_902617455.faa
annot/GCA_002026205/GCA_002026205.faa  annot/GCA_902551545/GCA_902551545.faa  annot/GCA_902617855/GCA_902617855.faa
annot/GCA_002026245/GCA_002026245.faa  annot/GCA_902551645/GCA_902551645.faa  annot/GCA_902618265/GCA_902618265.faa
annot/GCA_002812945/GCA_002812945.faa  annot/GCA_902552495/GCA_902552495.faa  annot/GCA_902618405/GCA_902618405.faa
annot/GCA_003209615/GCA_003209615.faa  annot/GCA_902552795/GCA_902552795.faa  annot/GCA_902618755/GCA_902618755.faa
annot/GCA_003209675/GCA_003209675.faa  annot/GCA_902552955/GCA_902552955.faa  annot/GCA_902618795/GCA_902618795.faa
annot/GCA_003210155/GCA_003210155.faa  annot/GCA_902555385/GCA_902555385.faa  annot/GCA_902619035/GCA_902619035.faa
annot/GCA_003212125/GCA_003212125.faa  annot/GCA_902555735/GCA_902555735.faa  annot/GCA_902619715/GCA_902619715.faa
annot/GCA_003212215/GCA_003212215.faa  annot/GCA_902557525/GCA_902557525.faa  annot/GCA_902619985/GCA_902619985.faa
annot/GCA_003212245/GCA_003212245.faa  annot/GCA_902557665/GCA_902557665.faa  annot/GCA_902620395/GCA_902620395.faa
annot/GCA_003213045/GCA_003213045.faa  annot/GCA_902558515/GCA_902558515.faa  annot/GCA_902620435/GCA_902620435.faa
annot/GCA_003213335/GCA_003213335.faa  annot/GCA_902559205/GCA_902559205.faa  annot/GCA_902620535/GCA_902620535.faa
annot/GCA_003213795/GCA_003213795.faa  annot/GCA_902559455/GCA_902559455.faa  annot/GCA_902620595/GCA_902620595.faa
annot/GCA_003213815/GCA_003213815.faa  annot/GCA_902561415/GCA_902561415.faa  annot/GCA_902620875/GCA_902620875.faa
annot/GCA_003213875/GCA_003213875.faa  annot/GCA_902563005/GCA_902563005.faa  annot/GCA_902621315/GCA_902621315.faa
annot/GCA_003213995/GCA_003213995.faa  annot/GCA_902565285/GCA_902565285.faa  annot/GCA_902622105/GCA_902622105.faa
annot/GCA_003214415/GCA_003214415.faa  annot/GCA_902565665/GCA_902565665.faa  annot/GCA_902622295/GCA_902622295.faa
annot/GCA_003217005/GCA_003217005.faa  annot/GCA_902565705/GCA_902565705.faa  annot/GCA_902622345/GCA_902622345.faa
annot/GCA_003278505/GCA_003278505.faa  annot/GCA_902567815/GCA_902567815.faa  annot/GCA_902622415/GCA_902622415.faa
annot/GCA_003278545/GCA_003278545.faa  annot/GCA_902567915/GCA_902567915.faa  annot/GCA_902622435/GCA_902622435.faa
annot/GCA_003278585/GCA_003278585.faa  annot/GCA_902568655/GCA_902568655.faa  annot/GCA_902622895/GCA_902622895.faa
annot/GCA_003278605/GCA_003278605.faa  annot/GCA_902568755/GCA_902568755.faa  annot/GCA_902623015/GCA_902623015.faa
annot/GCA_003278625/GCA_003278625.faa  annot/GCA_902569085/GCA_902569085.faa  annot/GCA_902623115/GCA_902623115.faa
annot/GCA_003278635/GCA_003278635.faa  annot/GCA_902570065/GCA_902570065.faa  annot/GCA_902623355/GCA_902623355.faa
annot/GCA_003278665/GCA_003278665.faa  annot/GCA_902570325/GCA_902570325.faa  annot/GCA_902623705/GCA_902623705.faa
annot/GCA_003278685/GCA_003278685.faa  annot/GCA_902570975/GCA_902570975.faa  annot/GCA_902623805/GCA_902623805.faa
annot/GCA_003278705/GCA_003278705.faa  annot/GCA_902571145/GCA_902571145.faa  annot/GCA_902633145/GCA_902633145.faa
annot/GCA_003278845/GCA_003278845.faa  annot/GCA_902571165/GCA_902571165.faa  annot/GCA_902634715/GCA_902634715.faa
annot/GCA_003278855/GCA_003278855.faa  annot/GCA_902571355/GCA_902571355.faa  annot/GCA_902635915/GCA_902635915.faa

Thx. Guillaume.