File tree Expand file tree Collapse file tree
Expand file tree Collapse file tree Original file line number Diff line number Diff line change @@ -12,7 +12,7 @@ REF=ecoliMG1655.fa
1212
1313
1414# r.namesorted.fastq.md5sum FORCE
15- r : r.sorted.bam.md5sum r.unsorted.bam.md5sum r.fastq.md5sum r.namesorted.bam.md5sum r.single.sorted.bam.md5sum r.single.unsorted.bam.md5sum r.single.fastq.md5sum r.single.namesorted.bam.md5sum r.repeat.unsorted.bam.md5sum r.repeat.sorted.bam.md5sum r.repeat.mixed.bam.md5sum r.repeat.fastq.md5sum r.fasta.md5sum r.single.noqual.bam.md5sum FORCE
15+ r : r.sorted.bam.md5sum r.unsorted.bam.md5sum r.fastq.md5sum r.namesorted.bam.md5sum r.single.sorted.bam.md5sum r.single.unsorted.bam.md5sum r.single.fastq.md5sum r.single.namesorted.bam.md5sum r.repeat.unsorted.bam.md5sum r.repeat.sorted.bam.md5sum r.repeat.mixed.bam.md5sum r.repeat.fastq.md5sum r.fasta.md5sum r.single.noqual.bam.md5sum r.single.noqual.fastq.md5sum FORCE
1616 for x in r* .md5sum; do echo $$ x $$ (cat $$ x); done
1717
1818# retrieve ref from internet if not already available
@@ -24,7 +24,7 @@ ${REF}.bwt: ${REF}
2424 bwa index ${REF}
2525
2626r1.fastq : ${REF}.bwt
27- ${WGSIM} -N 10 ${REF} r1.fastq r2.fastq
27+ ${WGSIM} -N 100000 ${REF} r1.fastq r2.fastq
2828
2929
3030r2.fastq : r1.fastq
@@ -104,6 +104,10 @@ r.fastq.md5sum: r1.fastq r2.fastq FORCE
104104r.single.fastq.md5sum : r1.fastq FORCE
105105 ${FASTQBIN} -P r1.fastq > r.single.fastq.md5sum
106106
107+ r.single.noqual.fastq.md5sum : r1.fastq FORCE
108+ ${FASTQBIN} -P -Q r1.fastq > r.single.noqual.fastq.md5sum
109+
110+
107111r.repeat.fastq.md5sum : r1.fastq r2.fastq FORCE
108112 ${FASTQBIN} r1.fastq r2.fastq r1.fastq r2.fastq > r.repeat.fastq.md5sum
109113
You can’t perform that action at this time.
0 commit comments