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
-
The immune system’s power comes from its genetic diversity, allowing myriad threats to be neutralized through first recognizing foreign antigens. That diversity is also what makes the immune system so difficult to study. Recent advances in sequencing technology and computational biology, however, are giving researchers new tools to understand immune responses and immune-related diseases in greater detail.
This convergence of genetics, immunology, and computation...09-01-2026, 05:41 AM
Comment