Question: Generate count table using Samtools
0
gravatar for vishalchanda364
3 months ago by
vishalchanda36420 wrote:

Hello,

I need to generate a count table for my metagenomic data (10 samples), in an article which I referred it says. Count table can be generated using Samtools, can anyone validate me about this, and may I know whether I can find the command for the same in Samtools manual.

Thanks, Vishal

alignment assembly • 236 views
ADD COMMENTlink modified 7 weeks ago by Biostar ♦♦ 20 • written 3 months ago by vishalchanda36420
2

samtools stats or samtools idxstats ?

ADD REPLYlink written 3 months ago by lieven.sterck7.8k
1

Thank you everyone for the help and support

ADD REPLYlink written 3 months ago by vishalchanda36420

Please include a link to the article you are talking about.

ADD REPLYlink written 3 months ago by Arup Ghosh2.4k

Count table of what? Reads mapped to each gene? contig?

ADD REPLYlink written 3 months ago by Asaf7.6k

You have included almost no information in your post so it is impossible to help you. Retry ? What did you align to ? Did you align data ?

ADD REPLYlink written 3 months ago by colindaven2.2k

Well sorry, I couldn't frame the question properly. The article which I am referring is

The Human Gut Virome Is Highly Diverse, Stable, and Individual Specific

it can be cited here, The question is that, I have similar work to the one done in this work, I too have 10 samples of metagenomic sequences. For the first sample I have done all the step (from pre-processing to aligning it with various database) for getting clean and bacterial free Viral Sequences and to reduce the work load in the next samples, those samples are pre-processed and aligned to the 1st clean sequence. Now with all the 10 sample they explained that a count table is generated using samtools v-0.1.19. I just want to ask, how.?

I might sound a bit lame with the concept, but till now I have only this much knowledge about bioinformatics, Sorry

Please if anyone can help now

Thanks, Vishal

ADD REPLYlink modified 3 months ago • written 3 months ago by vishalchanda36420
2

Try a newer version of samtools like 1.7+

The you can do what @lieven.sterck recommended above, maybe idxstats will be most helpful.

Program: samtools (Tools for alignments in the SAM format)
Version: 1.7 (using htslib 1.7)

Usage:   samtools <command> [options]

Commands:

  -- Statistics
     bedcov         read depth per BED region
     depth          compute the depth
     flagstat       simple stats ************
     idxstats       BAM index stats *******************
     phase          phase heterozygotes
     stats          generate stats (former bamcheck) ******************
ADD REPLYlink written 3 months ago by colindaven2.2k

okay Thank You so much @colindaven. I will do it now.

ADD REPLYlink written 3 months ago by vishalchanda36420
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1808 users visited in the last hour