Bioinformatics Asked on April 5, 2021
I’m trying to use the code specified in this link to create a .bam file without the need for a .sam file.
Here is the code I’m using:
bwa mem -R "@RGtID:$SAMPLEtPL:illuminatSM:$SAMPLE"
$REF_PATH/hg38.fasta
mapping/R1.fastq.gz
mapping/R2.fastq.gz |
samtools sort - |
tee mapping/$SAMPLE.bam |
samtools index - mapping/$SAMPLE.bam.bai
2>> log
and I get the following error:
[M::bwa_idx_load_from_disk] read 0 ALT contigs
`[E::main_mem] fail to open file ` '.
The ref path and file locations work when creating and converting a .sam to a .bam file.
Any ideas?
This is an error from bwa
saying that it can't find your input files:
Here's a checklist:
newline breaks. I.e.
is the very last character$REF_PATH
-> "$REF_PATH"
set -x
at the top of your bash script to debug.Correct answer by Bastian Schiffthaler on April 5, 2021
Why make this so complicated with the tee command?
I don't think you even can index an unsorted bam, so you have to wait until sorting is done before you index. Indexing is really short compared to mapping and sorting, you aren't saving much by trying to do both at once.
Answered by swbarnes2 on April 5, 2021
1 Asked on March 24, 2021 by timd1
2 Asked on March 23, 2021 by whateversclever
1 Asked on March 22, 2021 by swa_mi
1 Asked on March 22, 2021 by nitha
1 Asked on March 20, 2021
2 Asked on March 19, 2021 by lazer-guided-lazerbeam
2 Asked on March 19, 2021 by celinedion
1 Asked on March 19, 2021 by user3390486
1 Asked on March 16, 2021 by maxno3
0 Asked on March 13, 2021 by mendel
1 Asked on March 13, 2021 by ryan-ward
0 Asked on March 12, 2021 by user257566
1 Asked on March 11, 2021
Get help from others!
Recent Questions
Recent Answers
© 2023 AnswerBun.com. All rights reserved. Sites we Love: PCI Database, MenuIva, UKBizDB, Menu Kuliner, Sharing RPP, SolveDir