Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Relevant for me was lh3's comment in another thread:: Use the -A option with mpileup in order to include non-properly-paired reads.
-
I guess, you need to sort your bam file. i had similar problem and I got away with it by simply sorting it.
Leave a comment:
-
Too bad.
You could try to trawl the Samtools mailing list archives (e g http://sourceforge.net/mailarchive/f...=samtools-help) for further clues.
Leave a comment:
-
Thanks a lot for both suggestions. However, when I tried downgrading, all I got was a memory error, and no joy.
Not piping with -ugf gives the following output, but no output file.
[mpileup] 1 samples in 1 input files
?BCa??BCF#chromosom1chromosome2chromosome3 s7_2.bam%##samtoolsVersion=0.1.13 (r926:134)
<mpileup> Set max per-sample depth to 8000
I checked to see whether there was some problem with my SAM to BAM conversion by uploading my test file to Galaxy and running mpileup there, and everything seemed to work. So, could there be some weirdness in the way samtools are compiled on my system? I run Mac OS, but received no errors during compilation.
Leave a comment:
-
I also tried to run the command like that, and I got the same thing. It didn't write anywhere.
So I don't pipe. I run mpileup with -ugf to make a file, (and that step can take a very long time) and then I feed that file to bcftools to make the bcf once it's done.
The second step is so much faster than the first that you don't save much time by piping.
Leave a comment:
-
No output from samtool's mpileup
I am trying to run samtools for SNP detection, but neither mpileup or the deprecated pileup produces any output.
My data are in SAM format, which I convert to BAM using samtools thusly:
./samtools view -bSt FAIFILE.fai INFILE.sam > OUTFILE.bam
No errors, seem to get a functional BAM file.
Now, following the manual,
./samtools mpileup -ugf REFSEQ.fasta INFILE.bam | bcftools view -bvcg - > var.raw.bcf
This produces the following output to sdterr, and no output to screen
[mpileup] 1 samples in 1 input files
<mpileup> Set max per-sample depth to 8000
I have aslo tried a similar approach with pileup:
./samtools pileup -vcf REFSEQ.fasta INFILE.bam
This computes for LONG time, but still produces no output.
I have no idea what is going wrong... Please help!
Sasha
Latest Articles
Collapse
-
by seqadmin
The first FDA-approved CRISPR-based therapy marked the transition of therapeutic gene editing from a dream to reality1. CRISPR technologies have streamlined gene editing, and CRISPR screens have become an important approach for identifying genes involved in disease processes2. This technique introduces targeted mutations across numerous genes, enabling large-scale identification of gene functions, interactions, and pathways3. Identifying the full range...-
Channel: Articles
08-27-2024, 04:44 AM -
-
by seqadmin
Sequencing mRNA provides a snapshot of cellular activity, allowing researchers to study the dynamics of cellular processes, compare gene expression across different tissue types, and gain insights into the mechanisms of complex diseases. “mRNA’s central role in the dogma of molecular biology makes it a logical and relevant focus for transcriptomic studies,” stated Sebastian Aguilar Pierlé, Ph.D., Application Development Lead at Inorevia. “One of the major hurdles for...-
Channel: Articles
08-07-2024, 12:11 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 08-27-2024, 04:40 AM
|
0 responses
16 views
0 likes
|
Last Post
by seqadmin
08-27-2024, 04:40 AM
|
||
New Single-Molecule Sequencing Platform Introduces Advanced Features for High-Throughput Genomics
by seqadmin
Started by seqadmin, 08-22-2024, 05:00 AM
|
0 responses
293 views
0 likes
|
Last Post
by seqadmin
08-22-2024, 05:00 AM
|
||
Started by seqadmin, 08-21-2024, 10:49 AM
|
0 responses
135 views
0 likes
|
Last Post
by seqadmin
08-21-2024, 10:49 AM
|
||
Started by seqadmin, 08-19-2024, 05:12 AM
|
0 responses
124 views
0 likes
|
Last Post
by seqadmin
08-19-2024, 05:12 AM
|
Leave a comment: