Bash script to alignment of short reads against reference genome !
bwa mem -t 40 -R '@RG\tID:K12\tSM:K12' \ E.coli_K12_MG1655.fa SRR1770413_1.fastq.gz SRR1770413_2.fastq.gz \ | samtools view -b - >SRR1770413.raw.bam sambamba s...1569 days ago
To convert just one specific read group to fastq
# Stop script on error. set -uex # The SRR BioProject number for the sequencing data. PROJECT=PRJNA257197 # The number of datasets to subselect from the project....1552 days ago
Print the screenshot of read coverage remotely !
#Chrom_2:9700000-9760000 ASCIIGenome -ni --fasta A_vaga.NDPD.fasta -r Chrom_2:9700000-9760000 -x "save >> test.pdf" A_vaga.NDPD.fasta.out.sam.bam.sorted.bam1293 days ago
Corona variant calling steps !
#!/bin/bash FILE_PATH="data/trimmed_fastq_small/" find "$FILE_PATH" -name "*.fasta" | while IFS= read -r my_file do filename=$(basename "$my_file") basename=...1054 days ago
876 days ago
BBmap the reads with all alignments !
bbmap.sh in=../reference/reference.numbered.fa ambig=all vslow perfectmode maxsites=100000 out=fetch_Ids_for_barcode.sam827 days ago