Skip to content

Instantly share code, notes, and snippets.

@mbk0asis
Last active August 7, 2017 05:55
Show Gist options
  • Save mbk0asis/9d8136ca6c18cbc0923a86f60002d6f3 to your computer and use it in GitHub Desktop.
Save mbk0asis/9d8136ca6c18cbc0923a86f60002d6f3 to your computer and use it in GitHub Desktop.
# replace 100th character (strings and/or numbers) to "M"
$ sed 's/./M/100' input
###############################
$ awk '{ if ($0 ~ /^>/) { print $0; } \ # if the line starts with ">", print all
else { printf("%s%c%s\n", substr($0, 1, 2), "X", substr($0, 4, length($0))); } }' \ # else
additional.primers.cosmic.ext200.gDNA.fa | less
# INPUT - sequence must be a one line fasta
>APC_c.637C>T_77
aacagacacaatgcagtctCTTAAAAACAAAATCCCTTAATGCTCACAAACCAGTAACCCAGATACAGCCCAGACACTTCAATTCTCTCCATCTTGGAGTATGGTCATGGTCTCATTTTCCATCCTCAGTTTCCAAGAAGGCTGAAGTACCCTTATTTATCTAAGTTCTGAATTGGGTTTAGGGCACCAGATCTAGGGATGTTTGGGAACATCATATACTTAAAATTTCAgaagtcctagggagagcaatcaggcaagagaaataagtaaaaggaatccagataggaaaagaagaagtcaaactatctctctttgctgacaatataattctatacataaagaatcctaaagactccaccagaaggctcccagaactgataaacaatttcagtaaagtttc
>APC_c.646C>T_52
aggttcacgccattctcctgcctcagcctctccaatagctgggactacagtcacccgccaccatgcccagctaatttttttttttttttttttttttttttttagtagagacgaggcttcaccatgttagccaggatggtctcgatctcctgacctcgtgatccacctgcctcggcctcccaaaatgctggtattacaggcgtgagccactgcgcccggccaatgagcactacttttttttttttaattatacttttaagttctagggtacatgtgcacaacgtgcactacttttgagtgtcatgtcagtgctcaaaaagttttagattttggaggatttcagattttcagatttggggtgcctaacctCTATCTCAAGGCATCTATGAAGCAGGGAGGG
# OUTPUT - every third string in each sequence is replced with "X"
>APC_c.637C>T_77
aaXagacacaatgcagtctCTTAAAAACAAAATCCCTTAATGCTCACAAACCAGTAACCCAGATACAGCCCAGACACTTCAATTCTCTCCATCTTGGAGTATGGTCATGGTCTCATTTTCCATCCTCAGTTTCCAAGAAGGCTGAAGTACCCTTATTTATCTAAGTTCTGAATTGGGTTTAGGGCACCAGATCTAGGGATGTTTGGGAACATCATATACTTAAAATTTCAgaagtcctagggagagcaatcaggcaagagaaataagtaaaaggaatccagataggaaaagaagaagtcaaactatctctctttgctgacaatataattctatacataaagaatcctaaagactccaccagaaggctcccagaactgataaacaatttcagtaaagtttc
>APC_c.646C>T_52
agXttcacgccattctcctgcctcagcctctccaatagctgggactacagtcacccgccaccatgcccagctaatttttttttttttttttttttttttttttagtagagacgaggcttcaccatgttagccaggatggtctcgatctcctgacctcgtgatccacctgcctcggcctcccaaaatgctggtattacaggcgtgagccactgcgcccggccaatgagcactacttttttttttttaattatacttttaagttctagggtacatgtgcacaacgtgcactacttttgagtgtcatgtcagtgctcaaaaagttttagattttggaggatttcagattttcagatttggggtgcctaacctCTATCTCAAGGCATCTATGAAGCAGGGAGGG
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment