Dear colleagues,
I have the fastq files of paired end sequences of 42 animals but when I run bwa mem to generate the .sam files i did not get the correct output, it is empty.
bwa mem -M CHR6.fa file_1.fastq file_2.fastq > out_2R.sam.
Can anyone help me?
Thanks in advance
2 answers
Both your previous thread (BBmerge output question) and this one suffer from the same problem: the question leaves out important information. In both questions, the error probably arose from a previous step in your workflow (processing of the raw fastq files), but as you didn't provide details, the errors were not easily found - and still haven't been. You will have to describe in detail - providing all commands - how the fastq files have been processed so far. Also, stating your final goals may help - for example, in the other thread, it was pointed out merging reads is not necessary for calling SNPs.
You said you are a beginner, so, for the time being, the best course of action would be to follow one tutorial, from start to end. Try to find a tutorial which clearly states the versions of the programs being used, then you can use the same versions. Of course, it would be better if the tutorial is recent. You can certainly read several tutorials, but don't mix commands from them at first. If you do so, there will be plenty of additional hurdles along the way, such as different tutorials using different versions of the same programs, with incompatible command-line formats between these versions.
sorry, the link to the image:
Log in to answer this question.
-
bwa mem hangs after a few thousand reads
written by cee28 •
-
paired reads have different names
written by daffodil •
-
problem to aligning FASTQ through BWA aligner
written by modarzi
-
How to loop over multiple files for BWA MEM in bash for loop?
written by DNAngel
-
Skip orientation as there are not enough pairs by bwa mem
written by seta
-
Aligning Multiple paired end files together
written by David_emir
-
why HISAT2 and BWA BAM files differ in number of entries?
written by marongiu.luigi
-
BWA MEM vs BOWTIE2 , which is best?
written by John
-
Help in writing NGS pipeline (or: the road from fasta to vcf)
written by Hervé •
-
bwa mem mismatch option
written by Varun Gupta
what are the error/messages displayed ?
Why don't you also put in the command you used to index CHR6.fa, just to cover all the bases.
the fastq files has only sequences of chromosome 6
thanks¡
By the way, it's always a bad idea to align the read on only one chromosome rather than on the whole reference.
The commands run, this is the message that I can see in the screen, but the output file only has three lines
The image doesn't work. Can't you just cut and paste the output?
sorry, the link to the image:
That's an error. Did these fastq come from the sequencer like that or did you do so something to it?
This is a problem with the way your reads are ordered in the two fastq. Both files should have the read same names and in the very same order.