Hello everyone,
I have recently started to use bwa. As I was using this tool I tried to change the orders of the fastq files. So instead of
R1.fastq R2.fastq > result1.sam
I used
R2.fastq R1.fastq > result2.sam
vcf files generated from those sam files were different than each other. They had a little bit of difference but still I am wondering why they had different result.
Thank You
I do not agree with this. Reads being properly-paired is not influenced by the order of the two mate files towards each other. See e.g. with a selection of 100k human reads:
There is no information about first/second in pair in fastq files, this is assigned after alignment. Switching order should only change the strand information. Rather than that OP should give information about the variant calling, please provide the full command line. Maybe the differences come from low-quality alignments and/or low-complexity regions.
Edit: I would still recommend that you keep the correct fastq order at all times to avoid any potential issues that could happen. There is no reason to mess around with this.
thank you so much. can I ask what is what is OP. What commandline do you want me to provide.
First thanks a lot for your quick response. I'll try to learn about proper pairs. I have checked a little bit but if you know how does it change the outcome?