Question: Problems with STAR mapping output
0
gravatar for manekineko
5.4 years ago by
manekineko140
manekineko140 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.5k views
ADD COMMENTlink modified 5.4 years ago by mark.ziemann1.3k • written 5.4 years ago by manekineko140

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 modified 13 months ago by _r_am31k • written 5.4 years ago by Ashutosh Pandey12k

In addition to Ashutosh'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 modified 13 months ago by _r_am31k • written 5.4 years ago by mark.ziemann1.3k

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 modified 13 months ago by _r_am31k • written 5.4 years ago by manekineko140
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 5.4 years ago • written 5.4 years ago by Ashutosh Pandey12k
0
gravatar for mark.ziemann
5.4 years ago by
mark.ziemann1.3k
Australia/Mebourne/Geelong/Deakin
mark.ziemann1.3k 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 5.4 years ago by mark.ziemann1.3k
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: 2187 users visited in the last hour