Skip to content

Instantly share code, notes, and snippets.

perl -ne 'if (/^(\S+\/1)\:([ATCGN]+)\:(\S+)$/){print "@",$1,"\n",$2,"\n","+\n",$3,"\n";}' a.txt > a.fastq
#!/usr/bin/perl -sw
# sample inputs:
#ILLUMINA-A93428 66 1 1 7848 1267 CGATGT 1 ATGTTGGCTGGCGGTGAAATAAATCTCAAACGTACCTGTTACAAAATTGTGTTAATCCTTTCAGATTCGCAG ggggggggggggggcfdffdgggfefffffcfcffcccacBeddddded_cffefgggggBBBBBBBBBBBB chrX.fa 122287709 R 40C21GAC7 269 Y
# output for sample input:
#@ILLUMINA-A93428:66:1:1:7848:1267:#CGATGT/1
#ATGTTGGCTGGCGGTGAAATAAATCTCAAACGTACCTGTTACAAAATTGTGTTAATCCTTTCAGATTCGCAG
#+ILLUMINA-A93428:66:1:1:7848:1267:#CGATGT/1
#ggggggggggggggcfdffdgggfefffffcfcffcccacBeddddded_cffefgggggBBBBBBBBBBBB
# set up flags for Numpy C extentions compiling
export CFLAGS="-arch i386 -arch x86_64"
export FFLAGS="-m32 -m64"
export LDFLAGS="-Wall -undefined dynamic_lookup -bundle -arch i386 -arch x86_64"
export CC=gcc-4.2
export CXX="g++ -arch i386 -arch x86_64"
pip install numpy
# success!