caseywdunn / sharkmer

3 stars 0 forks source link

Issue with finding sPCR products #13

Open Deco313 opened 1 month ago

Deco313 commented 1 month ago

HI, I am testing Sharkmer's sPCR function on a library that I have previously retrieved a full mitogenome from using the standard primers for this group

sequence I am trying to recover GCTGCTAGTTGGTATTGGCATTTTGTAGATGTGGTTTGGTTGTTTTTATATGTATGTATTTACTGGTGGGGTTCTTAGCCCGGGCCATGGTTACATAATGCTTAGCTAAGCTTGTAGAGTCAACTAACGGTAAGTTTTCAGACTCATAATCTGATTATGCAGGTTCAACTCCTGCCTCTACCTATAGAATCGCCATGGGGAACTTATTCGCTCTGTATGCCCTTAGCAAGGGCCAAATACCTGTACAGTATTTTAACTTAATGGAGGAGAATTATTCTAAGTATGGGCTATCAGTAATCCAGCTTATCCAGATTGGTAAGTTCTATGAACTTTGGCAGAAGCCTGATGTTCCTAGTATCCAACAAGCATACTCTCAAGCTGAGTTATTAGTTGAGTCGTCCATACGAAGTAGGTCTTTAGGGGTAACACCTCCCATTGAACAAGTTGCTTCGTTACTTGATATGAGAATAACATCACCCGGTAAAAGATCCTTGCTTCAAATGGGGTTTCCTATTTATTCCCTAACTACTTATCTAAGCACCTTGTTGGATAAAGGTTGGACTATTATAGTTATTGATGAATTAGCCACTGATAAATCGGGGCCGAAACAACGCGCAGTATCTCAGGTTTATTCTCCTAGTTGTAATTTAGAAGATTGTCCGGAATTATCTTATGTGATATCAATCTATTTTAAAGATGACTTACTAGGTATTACTTTATTTTCTGCCATGAATGGGCACAGTATAATGTTTCCTGTCTATTGGGCCGACAGAGATAAAGTGGCTCGACTACTAATCAGTTATCGTATTAAAGAGGTAGTAATTAGGGCGGACTCGGAGGCCAATTATAGGGCACTAAATAAGATATATGATTTATTAATTGGTTGGAATTTATTCCCCTCTGAACCTAATGCCATAATTGAAGTTATGGAGGACCCGTGTTATTTATCTTATAGGTACGAAAATGATAATAAGGAGTGGCTTTTGCTTCA

Forward GCTGCTAGTTGGTATTGGCAT Reverse TSGAGCAAAAGCCACTCC

code:

zcat BAM_R*_P.fastq.gz | sharkmer -s BAM --max-reads 1000000 --pcr "GCTGCTAGTTGGTATTGGCAT,TSGAGCAAAAGCCACTCC,1200,mtMutS"

For gene mtMutS, no path was found from a forward primer binding site to a reverse binding site. Abandoning PCR. Suggested actions:

Is there any reason why this isn't working?