复制行我有这样类似项的文件:与SED和剪切命令
@HWI-ST750:151:C1C6AACXX:8:1101:10000:14114/1
AAACATACCATTGCATATTTCGTATAATTCGAATTTAAGCAAGTTGTATTGCCAACAGCTCAGATCAATTAAGAATCTTTTCAAATTTATCAATTTA
+/1
DDHBFHHIAHHJEGJFGICHG;CGHGIJIIAIHHIJIGAGIE9FDFHGHG:[email protected]=;@FHI:DACHE>>[email protected];[email protected];)6(>[email protected]>[email protected]
@HWI-ST750:151:C1C6AACXX:8:1101:10000:37802/1
GTTCCTTTACTTGCTCCCATGTCCCAATCTAACACACATGAGTCGTCGCACGCCAGTAAATCTAAAGATTTCAGACACTCAGACTCGTACAGTGAAA
+/1
[email protected]DEEEEDDDDDDDDDCDDCCCDDDABBDD>CC>A
我希望削减每2号线,从1日至45字。所以我用下面的命令:
sed "0~2p" BT8_A.fastq | cut -c1-45 | head -n 12
它工作正常,但削减1日45个字同一条线路后,得到追加到下一line.Like本 -
@HWI-ST750:151:C1C6AACXX:8:1101:10000:14114/1
AAACATACCATTGCATATTTCGTATAATTCGAATTTAAGCAAGTT
AAACATACCATTGCATATTTCGTATAATTCGAATTTAAGCAAGTT
+/1
DDHBFHHIAHHJEGJFGICHG;CGHGIJIIAIHHIJIGAGIE9FD
DDHBFHHIAHHJEGJFGICHG;CGHGIJIIAIHHIJIGAGIE9FD
@HWI-ST750:151:C1C6AACXX:8:1101:10000:37802/1
GTTCCTTTACTTGCTCCCATGTCCCAATCTAACACACATGAGTCG
GTTCCTTTACTTGCTCCCATGTCCCAATCTAACACACATGAGTCG
+/1
[email protected]
[email protected]
我不能够解决这个。我在哪里错了。帮我。
谢谢。
我所需的输出是:
@HWI-ST750:151:C1C6AACXX:8:1101:10000:14114/1
AAACATACCATTGCATATTTCGTATAATTCGAATTTAAGCAAGTT
+/1
DDHBFHHIAHHJEGJFGICHG;CGHGIJIIAIHHIJIGAGIE9FD
@HWI-ST750:151:C1C6AACXX:8:1101:10000:37802/1
GTTCCTTTACTTGCTCCCATGTCCCAATCTAACACACATGAGTCG
+/1
[email protected]
BINGO!谢谢你..确实是我想要的.. :) 但是上面的代码出了什么问题?! – PoisonAlien 2013-02-20 10:17:44
@poisonAlien:通常情况下,'sed'在执行命令后总是打印模式空间的内容(除非用'-n'告诉它)。所以'0〜2p'告诉它打印偶数行a * second *时间,除了打印它将要做的行。另外,你将所有'sed'的输出传递给'cut',所以奇数行也会被减少到前45个字符。 – OpenSauce 2013-02-20 13:22:41