Skip to content

Instantly share code, notes, and snippets.

@gatoravi
Created June 30, 2016 20:20
Show Gist options
  • Save gatoravi/40fd41f11ef9053e9777416c32fdfcdf to your computer and use it in GitHub Desktop.
Save gatoravi/40fd41f11ef9053e9777416c32fdfcdf to your computer and use it in GitHub Desktop.
cis_ase_generate_testdata
rm -f simulated_hets.bam && python simulate_bam_reads.py ~/dat/ref/hs37d5/all_sequences.fa 0.5 1:100000 1:100200 1:10040 0 1:100600 1:100800 | samtools view -bh > simulated_hets.bam && samtools index simulated_hets.bam
rm -f simulated_hets.bam && python simulate_bam_reads.py ~/dat/ref/hs37d5/all_sequences.fa 1.0 1:100000 1:100200 1:10040 0 1:100600 1:100800 | samtools view -bh > simulated_homs.bam && samtools index simulated_homs.bam
rm -f cis_ase_tumor_rna.bam && python simulate_bam_reads.py ~/src/regtools/tests/integration-test/data/fa/test_chr22.fa 0 22:1050 22:1550 22:1950 22:2050 22:2550 22:2950 22:3050 22:3550 22:3950 | samtools view -bh > cis_ase_tumor_rna.bam && samtools index cis_ase_tumor_rna.bam
rm -f cis_ase_tumor_dna.bam && python simulate_bam_reads.py ~/src/regtools/tests/integration-test/data/fa/test_chr22.fa 0.5 22:1050 22:1550 22:1950 22:2050 22:2550 22:2950 22:3050 22:3550 22:3950 | samtools view -bh > cis_ase_tumor_dna.bam && samtools index cis_ase_tumor_dna.bam
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment