grenaud / leeHom

Bayesian reconstruction of ancient DNA fragments
GNU General Public License v3.0
23 stars 14 forks source link

Hi #18

Open sultandromeda opened 1 month ago

sultandromeda commented 1 month ago

I apologize for such a question, but should we use all the generated FASTQ files (including failed ones to merge) for the alignment? Or do you suggest using only the successfully merged fastq files?

grenaud commented 1 month ago

We generally map the merged as single-end and the unmerged as paired-end, sort each and merge. You can do this in one line using file descriptors actually.