Skip to content
Open
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 3 additions & 1 deletion docs-site/src/gen-reads/fastq-shuffling.md
Original file line number Diff line number Diff line change
Expand Up @@ -3,7 +3,9 @@
real sequencer's FASTQ is in effectively random order. Variant calling does not depend on this,
since the BAM is coordinate-sorted before calling. Tools that read the FASTQ in batches,
aligners in particular, do see the difference: a batch can be filled entirely with reads from
one region, such as a long satellite array. Shuffle before aligning if that matters to you.
one region, such as a long satellite array. Measured on a 21x simulation of T2T-CHM13 (#870),
bwa-mem2 peaked at 166 GB on the reads as written and at 90 GB on the same reads shuffled,
with the same alignments either way. Shuffle before aligning a large genome.
#870 tracks shuffling in `eidolon` itself.

With `seqkit`, shuffle each mate separately **with the same seed**. seqkit has no paired mode,
Expand Down
Loading