Hello
Does anyone know tools to see the contents in PacBio movie files. I get movie files in three different formats: fasta, fastq and bam. But these files are different from the Illumina files and it is not possible to read the contents using the cat or head command.
Is the raw data from sequel instrument comes as bam files which is converted as movie files (fasta, fastq and bam).
Thanks
•
written
by
priya.bmg •
1 answer
Hello,
fasta & fastq files (unless they are gzipped) are regular text files and can be read with the cat or head command. In the case they are gzipped, you can still read them with the "zcat" command.
BAM files are binary compressed version of SAM files and can be read with the program samtools: https://github.com/samtools/samtools
There is a good tutorial on manipulating bam files here: http://quinlanlab.org/tutorials/samtools/samtools.html
The raw data comes off the instrument as a bam file and is subsequently converted to both .fasta and .fastq. Alot of people only need the fasta or fastq text representation of the sequencing data, but service providers will often provide all three to users because there is additional data in the *.bam file that certain analyses require.
Log in to answer this question.
-
Genome assembly reference .fasta + Illumina reads
written by A_heath •
-
pbmm2 alignment of long reads
written by priya.bmg •
-
Need an original PacBio SMRT Sequel fastq file from 2020
written by Michael Dondrup
-
Merging bowtie output
written by Sbrillo •
-
Analyzing PacBio sequel data and SRA
written by Smandape
-
Demultiplex PacBio CCS reads on Mac
written by manthony5955 •
-
How to split BAM files by samples from PacBio
written by misterie
-
PacBio methylation information in bam file
written by gabri
-
PacBio Quality checking tool
written by vahapel
-
SRA File Content
written by me •
PacBio data files should be no different than Illumina files except that each read may be insanely long compared to short illumina reads. They are still going to be following formats for files you mention.
There is one exception to make a note of.
Calculated Q scores are no longer provided for PacBio data (see this link):
This is only true for subreads data. Nearly all PacBio data is HiFi these days, which comes with PHRED scores.
Thanks for the clarification. I have not worked with PacBio data of late but remembered the page from a past discussion.
If you are with PacBio then you may want to consider making this information prominently available on company site. There is no definitive page to be found in a simple google search that lists this from official PacBio site e.g. "pacbio q score" via a google search.