Question: seqtk subseq command not pulling full fastq reads
23 months ago by
kratnasiri • 0
kratnasiri • 0 wrote:
When I run the seqtk subseq to extract sequences, it completely disregards the quality scores and the sequence.
>>>./seqtk subseq file.fastq IDs.txt @NS500126:798:HWTLTAFXX:1:11101:10771:1083:1-1 1:N:0:AATGGTCG+NGACCATT G + A @NS500126:798:HWTLTAFXX:1:11101:2195:1084:1-1 1:N:0:AATGGTCG+NGACCATT G + A @NS500126:798:HWTLTAFXX:1:11101:2581:1084:1-1 1:N:0:AATGGTCG+NGACCATT T + A
It's pulling the first character in the quality score and read; however, it's not pulling the full line.
Is there a way to fix that?
ADD COMMENT • link •