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
-
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 -
Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.
We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing...
Comment