This is a test version of Biostars. For the public version, visit https://www.biostars.org.
how to remove specific sequences from multi-fasta file which contains N ?

Hi EveryOne,

I have a multifasta file which contains 11000 (30kb each) genomes. Now i want to remove all the reads(whole genome) which contains N (minimum atleast one N ). How can I do this with sed or awk? Thanks in advance.

I have input like this :

Genome1 ATCGTCGTACAGATACAGATACANNNcGATAGACATAGACA

Genome2 AGTCGATCAGTACAGATACAGATACAGATACAGATAC

I want output like this

Genome2 AGTCGATCAGTACAGATACAGATACAGATACAGATAC

genome sequencing sequence alignment

Thanks @Pierre Lindenbaum, I have gone through that thread you mentioned, but it was not working fine with my large data, coz after executed that code, still the genome have Ns. Since that thread was 4 yrs old, i created my own thread. Can you help with this? Thanks

have you found a solution?

0 answers

No answers yet.

Log in to answer this question.