Mercurial > repos > iuc > bmtagger
diff test-data/gen.sh @ 0:49a1cbbe5767 draft default tip
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/main/tools/bmtagger commit e3d492d96b0ffe79370ca090b3f749b0869e8b60
| author | iuc |
|---|---|
| date | Wed, 12 Nov 2025 12:03:46 +0000 |
| parents | |
| children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/gen.sh Wed Nov 12 12:03:46 2025 +0000 @@ -0,0 +1,18 @@ +echo "> host" > host.fa +head -c 10000000 /dev/urandom | tr -dc 'ACGT' | head -c 100000 >> host.fa +echo "> contaminant" > contaminant.fa +head -c 10000000 /dev/urandom | tr -dc 'ACGT' | head -c 100000 >> contaminant.fa + +art_illumina -ss HS25 -i host.fa -p -l 75 -f 1 -m 200 -s 10 -o host.fq +art_illumina -ss HS25 -i contaminant.fa -p -l 75 -f 1 -m 200 -s 10 -o contaminant.fq + +cat host.fq1.fq contaminant.fq1.fq > host_and_contaminant.fq1.fq +cat host.fq2.fq contaminant.fq2.fq > host_and_contaminant.fq2.fq + +gzip -c host_and_contaminant.fq1.fq > host_and_contaminant.fq1.fq.gz +gzip -c host_and_contaminant.fq2.fq > host_and_contaminant.fq2.fq.gz + +# use 10-mers to reduce size +bmtool -d host.fa -o host.bitmask -w 10 +srprism mkindex -i host.fa -o host.srprism +makeblastdb -in host.fa -dbtype nucl \ No newline at end of file