I have a Multiple Sequence Alignment, where there are around 5000 sequences. There also exists many sequences where, there are so much of non-sequenced regions (for instance, AU----CGGGCA--NNNNNNNNNN).
So, by trial and error, I found a threshold value of 25%, so I removed the sequences which contained more than 25% of non-sequenced region (N's here), retaining around 400 sequences in the alignment, by writing a simple R script.
But, I would like to know, if there is any other standard method to do this, because I do not think, removing sequences based on finding a threshold for the percentage of non-sequenced regions, by trial and error is more efficient.If there is, can you please tell me, how to do?
Thanks in advance.