Tweet
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 field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...04-22-2024, 07:01 AM -
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...04-04-2024, 04:25 PM
Comment