search

lh3 / seqtk

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

Rename the read ID #99

Open jiangzy26 opened 6 years ago

jiangzy26 commented 6 years ago

I got a problem when i markduplicates using picardtools , it shows Values was put into PairInfoMap more than oce, i think it's the same reads ID problem , so i want to know how to handle this problem with using the command rename of seqtk, Thank you very much

tseemann commented 6 years ago

You can rename your reads as follows:

seqtk rename old.fq.gz prefix_ | gzip > new.fq.gz

This wil name all the reads prefix_1, prefix_2, .... , prefix_346341 etc