This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to separate plasmid proteins from main chromosome proteins in a GenBank assembly record?

I'm working on a mega-project involving all proteobacterial proteomes.Recently NCBI relocated all assembly record here. Previously, records for plasmids were kept separately from main chromosomes, but now they are placed in one file. For example, GCA_000010825.1_ASM1082v1_protein.faa.gz.

Question: how in that record separate which proteins came from plasmids and which are from chromosomes?

If I worked on 1 genome, I could have traced each protein individually, but with more than 2,000 complete genomes it's not going to be feasible. Also, I cannot rely on sequence annotations such as,for example "plasmid backbone"since not all plasmid proteins are necessary "plasmid backbone" proteins.

Any ideas?

assembly genome genbank
• 2,795 views • link
updated by Ram Ram • written by svetlana.lockwood • svetlana.lockwood

1 answer

The feature table file..

cut -f5 -d $'\t' GCA_000010825.1_ASM1082v1_feature_table.txt | sort | uniq -c
 5400 chromosome
 844 plasmid
 1 seq_type

Parse the IDs from there and then extract from the faa file..

awk -F '\t' '{if($5=="chromosome")print $11}' GCA_000010825.1_ASM1082v1_feature_table.txt |\
    grep . > GCA_000010825.1_ASM1082v1_chromosome_protein_coding.acc
• 0 views • link
updated by Ram Ram • written by 5heikki 5heikki

Log in to answer this question.

More posts like this