Extracting Split read alignments
0
0
Entering edit mode
14 months ago

I am trying to extract the split read alignments with bwa-mem for pair end sequencing. This is my code.

Extract the split-read alignments

samtools view -h rscm2.bam \
    | scripts/extractSplitReads_BwaMem -i stdin \
    | samtools view -Sb - \
    > rscm2.splitters.unsorted.bam

But I am getting this error message. /cm/local/apps/slurm/var/spool/job16557014/slurm_script: line 18: scripts/extractSplitReads_BwaMem: No such file or directory Can anyone please help to fix this error? Thank you in advance.

alignment • 486 views
ADD COMMENT
2
Entering edit mode

does the directory 'scripts' exists in the current directory ? try to use a full path for scripts/extractSplitReads_BwaMem instead of a relative path ?

ADD REPLY

Login before adding your answer.

Traffic: 1025 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6