2085 days ago
Samtools commands for bioinformatician !
## count mapped reads samtools view -c -F 260 mapping_file.bam ### converting sam file into fasta samtools fasta reads_mapped.sam > reads.fasta ### converting sam...1613 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=$(b...1055 days ago