I am converting the fq.gz. files (which are the results of the mgi study) to bam files to view on igv.
Hey everyone, before i start apologies for the inconvenience cause of my wrong or inappropriate use of terms. I take some fails of bwa mem lately. As i mentioned in the post title i'm working with fastq files. And i get:
[E::main_mem] fail to open file '/mnt/MGI-R4/L01/(sample_name)_(barcode_number)_1.fq.gz'.
and
[main_samview] fail to read the header from "-".
And here is my command:
bwa mem -M -t 8 ~/genome/hg19.fa /mnt/MGI-R4/L01/(samplename)_(lane)_(barcodenumber)_1.fq.gz /mnt/MGI-R4/L01/(samplename)_(barcodenumber)_2.fq.gz | samtools view -o (samplename)_(lane)_(barcodenumber).bam -
i know that my command works because i used to convert. Can you please help me?
• 4,156 views
•
link
updated
by
jkbonfield
•
written
by
canavar •
•
written
by
canavar •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
slurm batch
written by bestone •
-
For loop to go over all files in a folder to execute bwa mem
written by ravigenomeinfo •
-
BWA alignment/Samtools; Fail to read the header
written by Qboy •
-
[E::bwa_idx_load_from_disk] fail to locate the index files
written by melissachua90 •
-
transfer parameter to reader group in bwa mem command
written by Chenfei Zheng
-
How to pass arguments in a bash script to perform read alignment to the reference?
written by Karma
-
samtools sort error
written by evelyn
-
Sambamba-sort not enough data in stream
written by fiona.newberry
-
Running BWA parallel on 1000 samples
written by ebrown1955
-
BWA aligning help
written by lcc1844 •
is it the real command ? I mean, are the files literally named "(samplename)_(lane)_(barcodenumber)_1.fq.gz" (?!!)
Even if they are not canavar seems to be using them as such in command line (unless they have gone to the trouble of obfuscating the file name even in the error message).
canavar :
( )are place holders. You need to use real file names for the command to work. For example with a sample namedSample1your command will look like this:I am sorry i didn't really know that file names matters. There's nothing important to hide. I thought this makes it easier. Sorry again. Let me post it all again:
There is no need to
view. Following should work to create a sorted BAM directly.Thank you so much, i'll try this command.
As i mentioned before i am really a beginner. I had no such purpose. And thank you all for your answers already.
by the way i am using this command to sort:
and it goes like this:
to index the file:
and after that there comes
rmandcpphasesPlease use
ADD REPLY/ADD COMMENTwhen replying to existing posts. Do not useSUBMIT ANSWERunless you are posting a new answer for original question.You do not need to sort the BAM file if you use the command above since it will already be sorted.
samtoolsis smart about that. You will just need to Index the resulting sorted BAM file.Thank you and apologies.. Yes, you are right there's no need to sort for bam files. I am sorting files for bai files. And i have fix the error btw i guess that was related to my 1tb storage. Thanks alot again. Have a healty day kind sir.
You can also now add
--write-indexto modern Samtools versions so the index is written at the same time as the BAM (etc) file. It needs to be modern though! We've fixed a couple bugs in that code when combined with multi-threading. (It can be quite challenging to get indices to work there due to migrating reads to the next block if they don't fit.)