I've just aligned some MiSeq paired-end reads to a reference assembly generated in Velvet. The assembly has ~12,000 contigs in it, and when I view my alignment sam files, the header section includes information about each contig. However, when I convert my sam files to bam using samtools view, the log information for each conversion states that there are only 2500 sequences read in the SAM header. Anyone know what's happening?
-
-
Which version of samtools are you using? -
I'm using 0.1.18.Comment
-
1) Check whether you really have 2500 or 12000 reference header lines
2) Post the samtools commands you used and the full log report.Code:grep -E "^@SG" /path/to/your/SamFile.sam | wc -l
3) Try the latest version of samtools (v1.0)Comment
-
1b) Convert your bam file (the one which should only have 2500 seqs in the header) back to sam
and count againCode:samtools view -h yourBamFile.bam > newSamFile.sam
Comment
-
NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification,...10-30-2026, 06:46 AM -
Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?
Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating...10-14-2026, 09:55 AM
Comment