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.bam1277 days ago
Corona variant calling steps !
...me "$my_file" .fasta) echo $basename mkdir -p results/sam results/bam results/bcf results/vcf bwa...ed.sam samtools view -S -b results/sam/$basename.aligned.sam > results/bam/$basename.aligned.bam samto...1039 days ago
990 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.txt990 days ago
861 days ago
Extract the mapped and unmapped reads !
...#Single_End_Layout: samtools view --threads $PROCESSORS -b -F 4 in.bam > mapped.bam samtools view --threads $PRO...#Paired_End_Layout samtools view --threads $PROCESSORS -b -f 2 in.bam > mapped.bam samtools view -...584 days ago