The Frontline Support team will be offline February 18 for President's Day but will be back February 19th. Thank you for your patience as we get to all of your questions!
picard MergeBamAlignment function require paired reads
Does picard MergeBamAlignment function require paired reads as ALIGNED_BAM? This is not mentioned in the picard documents
However, we use trimmomatic to do reads QC, which generates three outputs: paired reads, unique reads after QC for 1st reads and 2nd reads. Then we use bwa to map three files individually to generate three sam files. Next, we merged bwa results into one bam file using picard MergeSamFiles. At last, we try to create a clean-up mapping file using picard MergeBamAlignment with a unmapped_bam file created from input fastq reads files using picard FastqToSam. The sort order of files are correct.
Does this mean picard and/or GATK does not work with a mixture of paired-reads and single-reads maps?