Tuesday, 29 November, 2022
[kevshaw@451]$ history
1 date
2 ls
3 cat hello.c
4 cat hello.c | wc -l
5 cat hello.c | wc
6 cd seq
7 ls
8 cd ..
9 cd dnaseq/
10 ls
11 samtools view yeast.bam | less
12 less yeast.bam
13 samtools view yeast.bam | less
14 samtools view -T yeast.fa -C -o yeast.cram yeast.bam
15 less yeast.cram
16 samtools mpileup -f yeast.fa yeast.cram | less
17 bcftools mpileup -Ou -f yeast.fa yeast.bam | bcftools call -vmO z -o snps.vcf.gz
18 ls
19 gzip -dv snps.vcf.gz
20 ls
21 less snps.vcf
22 grep ^[IVXM] snps.vcf
23 grep ^[IVXM] snps.vcf | less
24 grep ^[IVXM] snps.vcf | wc -l
25 grep ^[IVXM] snps.vcf | less
26 grep ^[IVXM] snps.vcf | grep INDEL
27 grep ^[IVXM] snps.vcf | grep INDEL | wc -l
28 grep ^[^#] snps.vcf
29 grep ^[^#] snps.vcf | less
30 grep ^[IVXM] snps.vcf | less
[kevshaw@451]$