1438 days ago
1211 days ago
Corona variant calling steps !
...em trimmed_fastq_small/$filename > results/sam/$basename.aligned.sam samtools view -S -b results/sam/$basename.aligned.sam > results/bam/$basename.aligned.bam samtools sort -o results/bam/$basename...1041 days ago
992 days ago
Commandline to Extract a list of specific read IDs from a bam file
#Save all the IDs in IDs.txt file #Run it on BAM file samtools view file.bam | fgrep -w -f IDs.txt992 days ago
Bash script to simulate a genome !
...nlm.nih.gov/genomes/all/GCA/003/401/745/GCA_003401745.1_ASM340174v1/GCA_003401745.1_ASM340174v1_genomic.fna.gz # Use BCFtools 1.10.2 and SAMtools 1.10 conda activate bcftools...914 days ago
863 days ago
Extract the mapped and unmapped reads !
PROCESSORS=20 #Single_End_Layout: samtools view --threads $PROCESSORS -b -F 4 in.bam > mapped.bam samtools view --threads $PROCESSORS -b -f 4 in.bam > unmapped.bam #Paired_End_Layout samtools view --threads $PROCESSORS -b...586 days ago