更改 fastq 文件中的“+”。 phred_offset错误qiime

Question

我的一些 fastq 读取有问题：

@1V3F_10526394 M01994:35:000000000-BM49D:1:1106:17684:21227 1:N:0:1 orig_bc=GGAATCTCTATAGCCT new_bc=GGAATCTCTATAGCCT bc_diffs=0
+
CGTACACTCCTGCGGGAGGCAGCAGTGGGGAATATTGCACAATGGGCGCAAGCCTGATGCAGCGACGCCGCGTGCGGGATGACGGCCTTCGGGTTGTAAACCGCTTTTGATCGGGAGCAAGCCTTCGGGTGAGTGTACCTTTCGAATAAGCACCGGCTAACTACGTGCCAGCAGCCGCGGTAATACGTAGGGTGCAAGCGTTATCCGGAATTATTGGGCGTAAAGGGCTCGTAGGCGGTTCGTCGCGTCCGGTGTGAAAGTCCATCGCTTAACGGTGGATCCGCGCCGGGTACGGGCGGGCTTGAGTGCGGTAGGGGAGACTGGAATTCCCGGTGTAACGGTGGAATGTGTAGATATCGGGAAGAACACCAATGGCGAAGGCAGGTCTCTGGGCCGTTACTGACGCTGAGGAGCGAAAGCGTGGGGAGCGAACAGGATTAGATACCCCTGTAGTCCC
+
CCCCCGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGDGGGDGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGFGGGGDGGCFGGGGGGGGGDGGGFGGGGGGDGGGGGGGGGGGGGCFG@CFGFGCFFGGGFGGFDFGGDGGGEFCGGCFGGGFGGGGGGDGGGGGFGGGGGGGGGGGDGGGGGGGFGDFFGGGGGGGGGGGGGGGGDECGGF7EEGGGGGFGGGGGGGGGGGGGFCGGGGEEGGGEEGGGGGGGF@CEGGGGGGGGGGGGGGGGGGGGFBGDGGGGFDGGGGGCGDGGGGGFGGGGGGGGGFGGGGGGGGGGGGGGGGGGGGEGGGGGGGGFCGGGGGGDGGGFGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGEGGFGGGGGGGGGGGGGGGGGGGGGGCCCCC

第一个“+”是个问题。我怎样才能删除它？

N.B。并不是所有的读取都有这个问题所以我不能删除符号“@”之后的每一行，因为我已经尝试过..

Answer 1

也许您可以尝试删除所有“+”行，然后每 3 行添加一个“+”行。这可能比尝试识别“+”是否正确放置更容易。

cat file.fastq | sed '/^+$/d' | awk '{print; if (NR%3==2){print "+"}}' > fixed.fastq

更改 fastq 文件中的“+”。 phred_offset错误qiime

change "+" in fastq file. phred_offset error qiime

fastq