1

Modifying a sequence using positions!

view story
linux-howto

http://www.unix.com – Hi.. i have two files one with positions information and another is sequence information. Now i need to read the positions and take the snps at the positions and replace that position base with the snp information in the sequence and write it in the snp information file.. for example Snp file contains 10 A C A/C Sequence file contains ATCGAACTCTACATTAC Here 10th element is T so i will replace T with [A/C] so the final output should be 10 A C A/C ATCGAACTC[A/C]ACATTAC Example files are Snp file SNPRefAlt10AC19GC30CT42AG While replacing the snps here from Ref and Alt column, we need to (HowTos)