Skip to content

Instantly share code, notes, and snippets.

@USMortality
Last active April 7, 2023 15:59
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save USMortality/6b22ef1bd147068687207e1a3f9c77a0 to your computer and use it in GitHub Desktop.
Save USMortality/6b22ef1bd147068687207e1a3f9c77a0 to your computer and use it in GitHub Desktop.
Drosten PCR Essay
#!/bin/bash
function check() {
echo "------------------------------------------------------------"
echo "$1"
seqkit amplicon -F "$1" "$2"
}
echo "------------------------------------------------------------"
echo "MN908947.1"
check "GTGARATGGTCATGTGTGGCGG" "out/MN908947.1.fasta"
check "CAGGTGGAACCTCATCAGGAGATGC" "out/MN908947.1.fasta"
check "CCAGGTGGWACRTCATCMGGTGATGC" "out/MN908947.1.fasta"
check "CARATGTTAAASACACTATTAGCATA" "out/MN908947.1.fasta"
check "ACAGGTACGTTAATAGTTAATAGCGT" "out/MN908947.1.fasta"
check "ACACTAGCCATCCTTACTGCGCTTCG" "out/MN908947.1.fasta"
check "ATATTGCAGCAGTACGCACACA" "out/MN908947.1.fasta"
check "CACATTGGCACCCGCAATC" "out/MN908947.1.fasta"
check "ACTTCCTCAAGGAACAACATTGCCA" "out/MN908947.1.fasta"
check "GAGGAACGAGAAGAGGCTTG" "out/MN908947.1.fasta"
echo "------------------------------------------------------------"
echo "MN908947.2"
check "GTGARATGGTCATGTGTGGCGG" "out/MN908947.2.fasta"
check "CAGGTGGAACCTCATCAGGAGATGC" "out/MN908947.2.fasta"
check "CCAGGTGGWACRTCATCMGGTGATGC" "out/MN908947.2.fasta"
check "CARATGTTAAASACACTATTAGCATA" "out/MN908947.2.fasta"
check "ACAGGTACGTTAATAGTTAATAGCGT" "out/MN908947.2.fasta"
check "ACACTAGCCATCCTTACTGCGCTTCG" "out/MN908947.2.fasta"
check "ATATTGCAGCAGTACGCACACA" "out/MN908947.2.fasta"
check "CACATTGGCACCCGCAATC" "out/MN908947.2.fasta"
check "ACTTCCTCAAGGAACAACATTGCCA" "out/MN908947.2.fasta"
check "GAGGAACGAGAAGAGGCTTG" "out/MN908947.2.fasta"
echo "------------------------------------------------------------"
echo "MN908947.3"
check "GTGARATGGTCATGTGTGGCGG" "out/MN908947.3.fasta"
check "CAGGTGGAACCTCATCAGGAGATGC" "out/MN908947.3.fasta"
check "CCAGGTGGWACRTCATCMGGTGATGC" "out/MN908947.3.fasta"
check "CARATGTTAAASACACTATTAGCATA" "out/MN908947.3.fasta"
check "ACAGGTACGTTAATAGTTAATAGCGT" "out/MN908947.3.fasta"
check "ACACTAGCCATCCTTACTGCGCTTCG" "out/MN908947.3.fasta"
check "ATATTGCAGCAGTACGCACACA" "out/MN908947.3.fasta"
check "CACATTGGCACCCGCAATC" "out/MN908947.3.fasta"
check "ACTTCCTCAAGGAACAACATTGCCA" "out/MN908947.3.fasta"
check "GAGGAACGAGAAGAGGCTTG" "out/MN908947.3.fasta"
echo "------------------------------------------------------------"
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment