File: generate_sequence.sh

package info (click to toggle)
jellyfish 2.3.1-5
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 3,276 kB
  • sloc: cpp: 35,703; sh: 995; ruby: 578; makefile: 397; python: 165; perl: 36
file content (23 lines) | stat: -rwxr-xr-x 620 bytes parent folder | download | duplicates (2)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
#! /bin/sh

mkdir -p tests-data; cd tests-data
. ../compat.sh

${DIR}/generate_sequence -v -o seq10m -m 10 -m 22 -s 3141592653 10000000
${DIR}/generate_sequence -v -o seq1m -s 1040104553 1000000 1000000 1000000 1000000 1000000

for i in 0 1 2 3 4; do
    gzip -c seq1m_$i.fa > seq1m_$i.fa.gz
done

${DIR}/generate_sequence -v -q -o seq10m -s 1473540700 10000000
ln -sf seq10m.fq seq10m.fastq
${DIR}/fastq2sam seq10m.fastq
if [ -n "$SAMTOOLS" ]; then
    $SAMTOOLS view -b -o seq10m.bam seq10m.sam
    $SAMTOOLS view -C -o seq10m.cram seq10m.sam
fi

if [ -n "$UNIX2DOS" ]; then
    $UNIX2DOS -n seq10m.fa seq10mDOS.fa
fi