Results for "Samtools"

Wire posts

  • Extract list of read ID from bam file: samtools view mbi.sorted.bam |perl -ne '/ID/ && print' > extractedID.sam #Extract #Ids #SAM #BAM #NGS

    3625 days ago

  • Count the contigs length in index fasta file with samtools. $ awk '{s+=$2}END{print s}' scaffolds.fasta.fai #Fasta #Index #Length #Count #Size #Samtools

    2134 days ago

  • #BAM file to fasta #conversion with #Samtool: samtools view filename.bam | \ awk '{OFS="\t"; print ">"$1"\n"$10}' - > filename.fasta

    2955 days ago

  • #BAM file to #fastaq #conversion with #Samtool: samtools view filename.bam | awk '{OFS="\t"; print "@"$1"\n"$10"\n+\n"$11}' >| filename.fastq

    2955 days ago

  • #Extract all #Reads from #BAM file for a #region Chr10:18000-45500 using #samtools: samtools view input.bam "Chr10:18000-45500" > output.bam

    2862 days ago

  • #Extract #tabbed info from #BAM file with samtools: samtools view test.bam|awk '{print $1 "\t" $2 "\t" $3 "\t" $4 "\t" $4+length($10)-1}' > file.xls

    2863 days ago

  • #Extract #Chromosome #Number And Start Position #Reads: samtools view bamfile.bam|awk '{print $3 "\t" $4 "\t" $4+length($10)-1}' > newfile.tab

    2857 days ago

  • Check if BAM file is sorted or not use: for F in `find . -name "*.bam"`; do echo $F; samtools index ${F} ; done #Sort #BAM #Trick

    2859 days ago

  • Extract reference nucleotide from BAM file: samtools view your.bam | awk '{print substr( $10, 100, 1)}' #SAM #BAM #Samtools #Tricks #NGS

    2856 days ago

  • Extract a fasta sequence from multifasta file: samtools faidx multitest.fa and then samtools faidx multitest.fa fasta_id > out.fa #Samtools #Fasta #Extract #Multifasta

    2738 days ago