This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Filter assembly contigs bei depth/coverage

Hi all,

is there an easy way to filter a contig fasta file which looks like this by the depth stated in the identifier e.g. make a new file with contigs with depth <2.0? Thanks in advance :)

>1 length=282707 depth=1.00x
TACGAAAAAAAAGCAAAAAAAAGAGGAGGAGGCTCACGCCTCCCCCTCTCCCTTTCAGCTTACGCCACCA
GCACGTCAGTCCGCGTCAACAGCATGCCCGTCAATCGGCGCCCTTATGTTGGGACATGATCAATGACAGG
ACACCCGGAACTTGTCTTTGTGCAGGCCTGTCCTCATACGTAAGGGGCTTTGCATCCCGCATCGCAATGA

>2 length=26815 depth=3.45x
CCAGGGAAGCAGCAGTCAACACGCTAAAAACCTGAAATCAGGCACAAGAGAAGTCAATCAAAGGTTTTTC
GAACCCTCCAGGTTTGACCCCGCAAAAGACCTCGCCAACCAAAAAAAGCACAAACTGTCTCTGGCCTTTG
GCGACCGTATTTTCGAGGACGATAATCACCTGATCTTGCCCACCATCCGCATTGAGGATGAGGAGGATCG

>3 length=24243 depth=1.02x
CACTTTTTTCGAAAAAACAATCTCCGCCTGAATTCGGATCGCGGTTAACTGCCGATTTCGACGGGCATCT
TTTTATCGTCCCCTTACCAGCAGCACCAGACGAAGTGGTATAATATAACAATAATCCGTTGGGCAGGATA
ACTGCTATTATCAAAATCAGCGCAGACCCTGCCAAAGACGGCCCGTAATATGCGGCCATCGGGGAGCGGC
CGCCGTAGCTGTTACGCCGTACTGCCCCAATCCATGTTCCATTCAAGGCAAGAGACGCACAAACCGATGC
AACCTGCTCCCAAGGCAACACTTCTCACGCGGCTCAATGAACAGATCGCCGTCAAGATGACCTTCATTTT
TGGAAGCATCTGGTGCGTTTACGCCTTCCTTGTCTTCTCCCTGATCCCGCTTGCCGCACCGGGATGGCAG
coverage fasta filter assembly contigs
• 1,073 views • link
written by thorbensieksmeyer • thorbensieksmeyer

1 answer 

awk -F '[=x ]' '/^>/ {P=($5 < 2.0)}  {if(P) print}' in.fasta
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this

No similar posts found.