This is a test version of Biostars. For the public version, visit https://www.biostars.org.
RNA Seq

I have a fasta file and I need to count the number of exact matched sequences in it. I tried

cat my.fasta | tr -d "\n" | grep "MySequence"

but it's not working for me.

rna fasta
• 1,812 views • link
updated by seidel seidel • written by far.zi • far.zi

How is this question related to RNAseq?

• 0 views • link
written by Ram Ram

How is it not working for you? Does grep not find the pattern that you can see is there? Do you need to create 1 line per fasta entry? (Multiline Fasta To Single Line Fasta)? If you have bowtie available, you could create an alignment index with your fasta, then supply "MySequence" on the command line, and count up the results.

• 0 views • link
written by seidel seidel

1 answer

I'd do

grep ">" myfile.fa | wc -l

That should just count the number of lines with ">", so line breaks won't matter.

• 0 views • link
written by swbarnes2 swbarnes2

Going off the grep in their command and the phrase "exact sequence" in their post, I think OP is trying to count number of exact matches to their query sequence.

• 0 views • link
written by Ram Ram

Yes, you're right. I need to count how many of an exact short sequence are there in my data. I could open it in TextWrangler and find but the file is huge and my mac is not helping.

• 0 views • link
written by far.zi • far.zi

Please edit your question and rephrase it so this is clear to everyone reading the question. This question has also been addressed multiple times on the forum, please use the search bar. If you want a quick lead, try exploring seqkit grep.

• 0 views • link
written by Ram Ram

Log in to answer this question.

More posts like this