Bash script to simulate a genome !
# Reference https://github.com/chhylp123/hifiasm/issues/33...e0.fasta.gz # make reference for randomreads.sh # randomreads.sh part of BBTools/BBMap http...max 12000bp /genetics/elbers/bbmap-38.86/randomreads.sh build=1 \ ow=t seed=1 \...912 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.sam812 days ago
Bash command to count reads in fastq file !
#ref fastq file #ref1_1.fq) echo $(( $(wc -l786 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 >...m #Paired_End_Layout samtools view --threads $PROCESSORS -b -f 2 in.bam >...584 days ago