Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
SRA issue m4merg Bioinformatics 2 04-08-2016 10:22 AM
SOAPdenovo issue szimmerman Bioinformatics 1 11-14-2013 02:36 PM
samtools mpileup tview issue shin_enomoto Bioinformatics 0 03-06-2013 07:17 PM
MV emPCR issue ajthomas 454 Pyrosequencing 1 10-23-2012 09:12 PM
SAMTools mpileup issue ashkot Bioinformatics 9 01-12-2012 05:22 PM

Thread Tools
Old 03-20-2017, 08:57 AM   #1
Junior Member
Location: USA

Join Date: Mar 2017
Posts: 2
Default Issue with mpileup

Hey y'all, I am encountering an issue while using mpileup to construct a consensus sequences from bam alignments. I had left the program to execute over the weekend but it was seemingly stuck on the first input file. No error message just:

[mpileup] 1 samples in 1 input files

This is strange because the bam file is not particularly large only a few mb. The command I executed is:
${z}/samtools-1.2/samtools mpileup -A -u -C50 -m3 -f ${bowtie2_path}.fa ${infile}.bam | ${z}/bin/bcftools view -cg - | vcf2fq > ${infile}.fq; {z}/seqtk-master/seqtk seq -A ${infile}.fq > ${infile}.consensus.fas; rm ${infile}.fq

Both samtools and bcftools are versions 1.2. Anyone know what's the cause of this?
quantic is offline   Reply With Quote

bioinformatics, mpileup snps bam multiple, samtools

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 08:44 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO