Extracting Split read alignments
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.
• 1,803 views
•
link
updated
by
Pierre Lindenbaum
•
written
by
prameelaawale •
•
written
by
prameelaawale •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
samtools view: failed to open for reading , samtools sort :can't open file
written by maria •
-
Samtools indexing problem during bwa mem alignment
written by GlouGlou •
-
Getting errors trying to run rmats
written by Omnia •
-
Samtools SORT or bedtools bamToFastq giving error when converting reads to fastq.
written by DNAngel
-
Error while running Jupyter script on HPC
written by Bioinfonext
-
Samtools segmentation fault when extracting exons from bam file.
written by jstevenson2256 •
-
Error in extracting Counts from the bam file using htseq-count
written by Biologist
-
Issues running RNAseqc on my bam file
written by vineelagangalapudi •
-
Extract reads where both pairs maps properly and none of the pairs are split from a bam file
written by rimjhim.roy.ch
-
Can't convert pair-end bam to bedpe using bedtools bamtobed
written by izzy.yichao.cai
does the directory 'scripts' exists in the current directory ? try to use a full path for scripts/extractSplitReads_BwaMem instead of a relative path ?