This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Counting Contigs in Fasta/GBK File

I wanted to check the no of contigs present in either a FASTA or GBK file, I am aware of algorithms such as CheckM that will allow for this process, however is there a direct code to check no of contigs in a sequence directly with python or biopython?

checkm contigs genome
• 6,399 views • link
updated by Joe Joe • written by biohacker_tobe • biohacker_tobe

you can try with basic utilities in *nix.

• 0 views • link
written by cpad0112 cpad0112

An easy grep solution to count entries in a genbank, is the number of LOCUS lines:

grep -c "LOCUS" multigenbank.gb

For a multifasta, you can use ^> instead of LOCUS as you have noted.

• 0 views • link
written by Joe Joe

2 answers

grep, sed, awk etc. Something like this:

$ cat test.fa 
>a
atgc
>b
atgc
>c
atgc

$ awk '/>/ {a++} END {print "number of sequences in this file: " a}' test.fa
number of sequences in this file: 3
• 0 views • link
written by cpad0112 cpad0112

yeah, I just tried this command this helped for determining the number of contigs per file (Just change the extension file for both cases):

Individual File:

$grep -c "^>" Streptomyces_sp_12.fna

Multiple Files :

$ grep -c "^>" *.fna

• 0 views • link
written by biohacker_tobe • biohacker_tobe

If an answer was helpful, you should upvote it; if the answer resolved your question, you should mark it as accepted. You can accept more than one if they work.
Upvote|Bookmark|Accept

• 0 views • link
written by GenoMax GenoMax

Easy in BioPython.

 from Bio import SeqIO
 recs = list(SeqIO.parse('genbank.gbk', 'genbank'))
 len(recs)

This could be more memory efficient with an iterator, but this is a quick and easy way.

This is likely a more robust solution too, since *nix solutions require that you know your files very well, such that they don't have any nasty surprises in them.

• 0 views • link
written by Joe Joe

Log in to answer this question.

More posts like this