Statistics The Number Of Identical Sequence Between Different Samples
I have three samples,they are small RNAs and these three or more samples were sequenced by Solexa high-thoughput sequencing.They are fasta format,the length of reads are 18-28nt.I want to statistics the number of identical sequence between different samples . their names are a.fa,b.fa,c.fa, For example.the a.fa,b.fa,c.fa there is a similar quention,but I don't know how modify it, I hope you can do it .Statistics The Number Of Identical Sequence Between Different Samples
a.fa
>1_x2
ATCG
>2_x3
ACTG
>3_x1
GAAG
b.fa
>a_x5
GAAG
>b_x3
ATCG
c.fa
>c_x1
ACTG
>2_x2
AAG
The result is like below.
sequence a b c
ATCG 2 3 0
ACTG 3 0 1
GAAG 1 5 0
AAG 0 0 2
• 2,324 views
•
link
0 answers
No answers yet.
Log in to answer this question.
Can you explain how you derive the matrix values in your result, given the sample input? It's unclear how you are parsing input to get to that matrix.
So far as I can tell, the "similar question" is in fact identical to this one. So this one should be closed.
I agree--same question as before.
I am sorry, but it is not our problem that you cannot figure out how to edit your question. Also, I have asked the user before to avoid putting stuff on this good<dot>gd and use a different service, (one that doesn't include flashy porn ads to host a zip file, and that most people will able to read) for example pastebin. I might want to put a ban on the user for reapeatedly ignoring this.
So the slight modification for this new question seems to be the "x2" or "x5" multiplier in the read names that should increase the read count by e.g. 2 or 5 respectively.
It's a terribly basic programming question and from the last question it seems that the poster just wants someone to do their coding for them since they cannot...
Just for the fun (thanks to matted for pointing the modification):