This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Split the sequence into each SNP entry

Hi,

I have a file from SNP data and would like to add

  1. from each sequence split the sequence into two with each SNP entry - this is bit tricky to do

  2. Add fasta symbol before each sequence start then. I know how to do this one using shell - sed 's/^\([^acgt]\)/>\1/'.

     AAGGGTTTAGAAAAAAACCAAACAAACAATCGAAA[C/T]GAAATAGAAAAAGAAAAAGGGAAGGGGTTAAGTTC
 TTCATATAAAAATTGATATAGAATCTTTGAAAAAG[A/C]CCTTTCTTCCTAAGAAAGAAAAGGCTTACTGTCTT
 CCCAAATAAACAGGTATGGAAGCTATAATTGGAAA[C/T]CACGATCGAATTTATGGAAGCATTGGTTTATACAT
 GGATCCAAAAGAAACTTGGGCATTTATTACTTGGA[C/T]GATATTCGGGATTTATTTACATACTCGAACAAATA
 TATCAGTTAGTCTACCATATTTTTTTCTTGACAGA[A/C]AACTAAGGAAATGGCTCCATGTGCTCTAATTCATT
 ACTAACTCTAAAGTAAAGGATCTTTCCACCTTTTC[G/T]GATCCCATACCAATAGCTTTTTTTGATTCGTCCAT
 AGTTTACACACTTTTGTATTACCTCTTCTTACTGC[C/T]GTATTTATGTTAATGCATTTCCTAATGATACGTAA
 AATAGATCTGACAAGTCGCACTATATGTCAACCCA[A/C]GATGGATGCTTGTCCCCGGGACTTCGATAAGGTAC

Thank you

sequence snp
• 1,137 views • link
updated by cpad0112 cpad0112 • written by K • K

1 answer

like this?

$ awk -F "[][/]" -v OFS="\n" '{print ">seq_"NR"_"$2,$1$2$4,">seq_"NR"_"$3,$1$3$4}' test.txt 

>seq_1_C
AAGGGTTTAGAAAAAAACCAAACAAACAATCGAAACGAAATAGAAAAAGAAAAAGGGAAGGGGTTAAGTTC
>seq_1_T
AAGGGTTTAGAAAAAAACCAAACAAACAATCGAAATGAAATAGAAAAAGAAAAAGGGAAGGGGTTAAGTTC
>seq_2_A
TTCATATAAAAATTGATATAGAATCTTTGAAAAAGACCTTTCTTCCTAAGAAAGAAAAGGCTTACTGTCTT
>seq_2_C
TTCATATAAAAATTGATATAGAATCTTTGAAAAAGCCCTTTCTTCCTAAGAAAGAAAAGGCTTACTGTCTT
>seq_3_C
CCCAAATAAACAGGTATGGAAGCTATAATTGGAAACCACGATCGAATTTATGGAAGCATTGGTTTATACAT
>seq_3_T
CCCAAATAAACAGGTATGGAAGCTATAATTGGAAATCACGATCGAATTTATGGAAGCATTGGTTTATACAT
>seq_4_C
GGATCCAAAAGAAACTTGGGCATTTATTACTTGGACGATATTCGGGATTTATTTACATACTCGAACAAATA
>seq_4_T
GGATCCAAAAGAAACTTGGGCATTTATTACTTGGATGATATTCGGGATTTATTTACATACTCGAACAAATA
>seq_5_A
TATCAGTTAGTCTACCATATTTTTTTCTTGACAGAAAACTAAGGAAATGGCTCCATGTGCTCTAATTCATT
>seq_5_C
TATCAGTTAGTCTACCATATTTTTTTCTTGACAGACAACTAAGGAAATGGCTCCATGTGCTCTAATTCATT
>seq_6_G
ACTAACTCTAAAGTAAAGGATCTTTCCACCTTTTCGGATCCCATACCAATAGCTTTTTTTGATTCGTCCAT
>seq_6_T
ACTAACTCTAAAGTAAAGGATCTTTCCACCTTTTCTGATCCCATACCAATAGCTTTTTTTGATTCGTCCAT
>seq_7_C
AGTTTACACACTTTTGTATTACCTCTTCTTACTGCCGTATTTATGTTAATGCATTTCCTAATGATACGTAA
>seq_7_T
AGTTTACACACTTTTGTATTACCTCTTCTTACTGCTGTATTTATGTTAATGCATTTCCTAATGATACGTAA
>seq_8_A
AATAGATCTGACAAGTCGCACTATATGTCAACCCAAGATGGATGCTTGTCCCCGGGACTTCGATAAGGTAC
>seq_8_C
AATAGATCTGACAAGTCGCACTATATGTCAACCCACGATGGATGCTTGTCCCCGGGACTTCGATAAGGTAC
• 0 views • link
written by cpad0112 cpad0112

thank you so much. I really appreciate your response. It is working great!!!

• 0 views • link
written by K • K

Log in to answer this question.

More posts like this