Confused about the number of reads #1924
Unanswered
sodiumnitrate
asked this question in
Q&A
Replies: 2 comments
-
Hi Irem, Mapping long reads (>300b) requires STARlong executable. |
Beta Was this translation helpful? Give feedback.
0 replies
-
Thank you for the info. I compiled and used
Here's the command:
The stdout is as follows:
|
Beta Was this translation helpful? Give feedback.
0 replies
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
-
Hi, I'm trying to map sequences from a
.fasta
file. The file contains about 2,800,000 sequences. However, the output alignment file has 259 reads. Is STAR somehow merging similar sequences? TheLog.final.out
file says only 30 are unmapped (see below).Further, the sequences in the
.fasta
file have names in the following formatname/XXXX/ccs
whereXXXX
is a number. However, the outputSAM
file has all the read names listed asname
. I tried adding-readNameSeparator "_"
but that didn't solve the problem.Here's the command I ran:
What could I be missing? There were no error messages in the Log files or otherwise. This is
Log.final.out
:Beta Was this translation helpful? Give feedback.
All reactions