lh3 / seqtk

Toolkit for processing sequences in FASTA/Q formats
MIT License
1.35k stars 310 forks source link

seqtk_telomere_reads_output #204

Open Wenfei-Xian opened 8 months ago

Wenfei-Xian commented 8 months ago

In the original seqtk telo, it only outputs the sequence ID and the telomere coordinates. I modified it because I needed to extract both the original sequence and the corresponding telomere-removed sequence (to avoid the influence of the telomere sequence on alignments).