Is there a tool to make counting of base in each position of reference genome ?
Example:
1 2 3 4 5 6 7
A C T T C G A :reference
A C G T A C A :read1
A A T T G C A :read2
A C G T G G A :read3
A C T T G A A :read4
A C A A G T A :read5
the output of the tool that I'm searching :
5A 0C 0T 0G: position1
1A 4C 0T 0G: position2
1A 0C 2T 2G: position3
1A 0C 4T 0G: position4
1A 0C 0T 4G: position5
1A 2C 1T 1G: position6
5A 0C 0T 0G: position7
Example:
1 2 3 4 5 6 7
A C T T C G A :reference
A C G T A C A :read1
A A T T G C A :read2
A C G T G G A :read3
A C T T G A A :read4
A C A A G T A :read5
the output of the tool that I'm searching :
5A 0C 0T 0G: position1
1A 4C 0T 0G: position2
1A 0C 2T 2G: position3
1A 0C 4T 0G: position4
1A 0C 0T 4G: position5
1A 2C 1T 1G: position6
5A 0C 0T 0G: position7
Comment