<h3>Q1</h3>
<p>Probably not - this is probably just displaying the contents of the input (I presume a <code>bam</code> file) exactly as they are encoded.</p>
<h3>Q2</h3>
<p>The splicing looks a little weird, but <strong>the exon-mapping pattern in the 5' part of the gene makes it seem like RNA rather than DNA</strong>. I think that the large intron is what makes you think that it is DNA, correct?</p>
<ul>
<li>Explanations for weird splicing
<ul>
<li>The weird read is a rare splice variant but a true mRNA.</li>
<li>The weird read is not completely matured and has an intron that has not yet been spliced. The IsoSeq library prep protocol picked this up somehow and it ended up being sequenced.</li>
</ul>
</li>
</ul>
<h2>Q3</h2>
<ul>
<li>Explanations for why <strong>Seq length: 6751</strong> and <strong>Seq length on ref: 5267</strong>
<ul>
<li>There is probably a large amount of clipping on the 5' or 3' end of the read. This could be polyA or some other transcribed-but-not-yet-completely-matured seq. <strong>Or</strong>, the extra sequence could be sequencing error, PacBio adapters, or a chimeric molecule.</li>
</ul>
</li>
</ul>
<h3>things to try</h3>
<p>How to check what is going on:</p>
<p>I'm guessing that you had to upload a <code>bam</code> file to JBrowse to make the IsoSeq track. Use the sequence name that you got from Jbrowse and run the following command to extract the alignment info from the bam file:</p>
<pre><code>samtools view /path/to/your/isoseq/file.bam | grep 'weird_read_id'
</code></pre>
<p>Look at the <a href="http://bioinformatics.cvr.ac.uk/blog/tag/cigar-string/" rel="noreferrer">cigar string</a> to see what it says about 5' or 3' clipping and corroborate with the sequence data.</p>

CC BY-SA 4.0
2018-12-13T16:24:54.133Z
you are right the first 3702 bases are soft clipped.
user1238097
CC BY-SA 4.0
2018-12-13T16:53:53.780Z
I'm curious to know what the clipped sequence is once you figure it out!
conchoecia
CC BY-SA 4.0
2018-12-17T22:47:32.280Z
I am still a mystified why it is clipping.  I am starting to notice that many of the reads that display this property look chimeric where the start site is seen twice in the middle of the sequence. I also note that this was the 4kb and up fraction.  So it might be assembly error you mentioned earlier I am feeling.
user1238097
CC BY-SA 4.0
2018-12-17T23:45:35.687Z
Thanks for following up. Please feel free to post another question about the assembly error if you'd like!
conchoecia
CC BY-SA 4.0
2018-12-13T00:06:42.297Z
2018-12-13T00:17:09.517Z
2020-06-18T08:30:41.743Z
conchoecia
<p>I recently had some Iso-Seq sequencing done on my organism catfish on the new Sequel platform and got weird alignments for a size selected 4 Kilobase and up fraction after running the isoseq3 pipeline.  When I investigated the exon intron junctions using data from (Pac Bio Iso-Seq, Illumina RNA-seq) against our private reference genome, I get PacBio Iso-Seq reads where the sequence length is greater than the seq length on the reference for my bam alignments:</p>

<p><strong>Seq length:</strong> 6751</p>

<p><strong>Seq length on ref:</strong> 5267</p>

<p>When blasting the sequence in NCBI it shows 5% similarity to the publicly-available genome.  I'll attach a JBrowse image with the same RNA sequenced with Illumina paired RNA-seq.  I also included a image and circled the pacbio contig and Illumina is above.</p>

<p><a href="https://i.sstatic.net/vtf1j.jpg" rel="noreferrer"><img src="https://i.sstatic.net/vtf1j.jpg" alt="Pac Bio Iso-seq and Illumina RNA-seq"></a></p>

<h2>Question 1:</h2>

<p>Is this a problem with JBrowse?</p>

<h2>Question 2:</h2>

<p>Does this mean that some DNA was sequenced with the catfish cDNA?</p>

<h2>Question 3:</h2>

<p>Why is this aligned sequence's length longer than the length of the reference genome region to which it is aligned?</p>

<p>Thanks</p>

<rna-seq><phylogenetics><pacbio>
Why is a PacBio read length larger than the aligned reference region?
# Why is a PacBio read length larger than the aligned reference region?

(Asked by: username_0 on 2018-12-12)

I recently had some Iso-Seq sequencing done on my organism catfish on the new Sequel platform and got weird alignments for a size selected 4 Kilobase and up fraction after running the isoseq3 pipeline. When I investigated the exon intron junctions using data from (Pac Bio Iso-Seq, Illumina RNA-seq) against our private reference genome, I get PacBio Iso-Seq reads where the sequence length is greater than the seq length on the reference for my bam alignments:

**Seq length:** 6751

**Seq length on ref:** 5267

When blasting the sequence in NCBI it shows 5% similarity to the publicly-available genome. I'll attach a JBrowse image with the same RNA sequenced with Illumina paired RNA-seq. I also included a image and circled the pacbio contig and Illumina is above.

[![Pac Bio Iso-seq and Illumina RNA-seq](https://i.sstatic.net/vtf1j.jpg)](https://i.sstatic.net/vtf1j.jpg)

## Question 1:

Is this a problem with JBrowse?

## Question 2:

Does this mean that some DNA was sequenced with the catfish cDNA?

## Question 3:

Why is this aligned sequence's length longer than the length of the reference genome region to which it is aligned?

Thanks

**Tags:** rna-seq, phylogenetics, pacbio

---

## 1 answers

### Answer by username_1 on 2018-12-13. Score: 5 | Accepted answer

Q1
Probably not - this is probably just displaying the contents of the input (I presume a bam file) exactly as they are encoded.
Q2
The splicing looks a little weird, but the exon-mapping pattern in the 5' part of the gene makes it seem like RNA rather than DNA. I think that the large intron is what makes you think that it is DNA, correct?

Explanations for weird splicing

The weird read is a rare splice variant but a true mRNA.
The weird read is not completely matured and has an intron that has not yet been spliced. The IsoSeq library prep protocol picked this up somehow and it ended up being sequenced.



Q3

Explanations for why Seq length: 6751 and Seq length on ref: 5267

There is probably a large amount of clipping on the 5' or 3' end of the read. This could be polyA or some other transcribed-but-not-yet-completely-matured seq. Or, the extra sequence could be sequencing error, PacBio adapters, or a chimeric molecule.



things to try
How to check what is going on:
I'm guessing that you had to upload a bam file to JBrowse to make the IsoSeq track. Use the sequence name that you got from Jbrowse and run the following command to extract the alignment info from the bam file:
samtools view /path/to/your/isoseq/file.bam | grep 'weird_read_id'

Look at the cigar string to see what it says about 5' or 3' clipping and corroborate with the sequence data.

> **Comments:**
> **username_0**: you are right the first 3702 bases are soft clipped.  
> **username_1**: I'm curious to know what the clipped sequence is once you figure it out!  
> **username_0**: I am still a mystified why it is clipping.  I am starting to notice that many of the reads that display this property look chimeric where the start site is seen twice in the middle of the sequence. I also note that this was the 4kb and up fraction.  So it might be assembly error you mentioned earlier I am feeling.  
> **username_1**: Thanks for following up. Please feel free to post another question about the assembly error if you'd like!  

---

