Question: Problems with STAR mapping output
0
gravatar for manekineko
4.2 years ago by
manekineko130
Bulgaria
manekineko130 wrote:

Hi,

I have problems seeing the output from STAR mapper. I have done the mapping, then made an sorted bam for IGV but I see very few reads mapped, and I have this log:

 Number of input reads |       116876253
                      Average input read length |       20
                                    UNIQUE READS:
                   Uniquely mapped reads number |       93912191
                        Uniquely mapped reads % |       80.35%
                          Average mapped length |       19.85
                       Number of splices: Total |       129488
            Number of splices: Annotated (sjdb) |       0
                       Number of splices: GT/AG |       129455
                       Number of splices: GC/AG |       33
                       Number of splices: AT/AC |       0
               Number of splices: Non-canonical |       0
                      Mismatch rate per base, % |       0.00%
                         Deletion rate per base |       0.00%
                        Deletion average length |       1.00
                        Insertion rate per base |       0.00%
                       Insertion average length |       1.00
                             MULTI-MAPPING READS:
        Number of reads mapped to multiple loci |       20734865
             % of reads mapped to multiple loci |       17.74%
        Number of reads mapped to too many loci |       2072050
             % of reads mapped to too many loci |       1.77%
                                  UNMAPPED READS:
       % of reads unmapped: too many mismatches |       0.00%
                 % of reads unmapped: too short |       0.00%
                     % of reads unmapped: other |       0.13%

 

I'm particullary interested to see the reads that are splitted in exon-intron sites (Number of splices: GT/AG |       129455)

I have noticed that when converting from SAM to BAM the size decreased from 16GB to only 60MB which is strange?!

rna-seq star • 3.0k views
ADD COMMENTlink modified 4.1 years ago by mark.ziemann1.2k • written 4.2 years ago by manekineko130

Well 16 GB to 60 Mb is strange. Have you checked if your BAM file is correct? Try running "samtools view Input.bam" and see if it throws "EOF marker is absent" error.

ADD REPLYlink written 4.2 years ago by Ashutosh Pandey11k

In addition to Ashutosh Pandey's recommendation, can you do a 

samtools view -c input.bam

 to count the rows in the bam and 

samtools view -c -q10  input.bam

to count the uniquely mapped reads. If you have fewer than 93912191 reads then the sam->bam conversion failed somehow.

ADD REPLYlink written 4.1 years ago by mark.ziemann1.2k

it returned error:

 

samtools view -c -q10 Aligned.out.bam
[bam_header_read] EOF marker is absent. The input is probably truncated.
[main_samview] truncated file.

but I made the conversion as I always did (of course I tried STAR for the first time, but it should output a standart SAM?)

I did it with - Samtools view -b -S SAM > BAM

ADD REPLYlink written 4.1 years ago by manekineko130
1

Yes, almost all the aligner these days will output standard SAM. I think there was a problem while converting your sam file to bam file. You should check the log. It could be insufficient space on your system, or a malformatted line in the sam file. The error "EOF marker is absent" suggests that the bam conversion process terminated before the complete conversion. 

ADD REPLYlink modified 4.1 years ago • written 4.1 years ago by Ashutosh Pandey11k
0
gravatar for mark.ziemann
4.1 years ago by
mark.ziemann1.2k
Australia/Mebourne/Geelong/Deakin
mark.ziemann1.2k wrote:

STAR 2.4 can output alignments directly in binary BAM format, so I'd recommend an update.

https://github.com/alexdobin/STAR/blob/master/doc/STARmanual.pdf

ADD COMMENTlink written 4.1 years ago by mark.ziemann1.2k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 725 users visited in the last hour