You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
I want to convert my bam file mapped on genomic assembly (bam by bowtie2) to my predicted genes (gff by prokka). Next I want to use this bam as salmon input.
But when i use the mudskipper bulk command line no reads are mapped on the transcriptome.
I thinks, they have a bug when the gff is loading because This information are logged :
2024-04-18T16:04:35Z INFO mudskipper::bam] Number of reference sequences: 0
I also test our gff in tests folder and i have the same results.
mudskipper bulk -t 1 -a mudskipper/tests/NC_002333.2.sam -g mudskipper/tests/NC_002333.2.gff -o demo.bam
[2024-04-18T16:04:35Z INFO mudskipper] Mudskipper started...
[2024-04-18T16:04:35Z INFO mudskipper::annotation] reading the gff file and building the tree.
█████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████ 102/102
[2024-04-18T16:04:35Z INFO mudskipper::annotation] building the tree
[2024-04-18T16:04:35Z INFO mudskipper::annotation] Time to build the tree: 6.540013ms
[2024-04-18T16:04:35Z INFO mudskipper::bam] Number of reference sequences: 0
[2024-04-18T16:04:35Z INFO mudskipper::bam] thread count: 1
[2024-04-18T16:04:35Z INFO mudskipper] Mudskipper finished.
When i use the gtf test file reference sequence are detected and sequence are mapped:
mudskipper bulk -t 1 -a mudskipper/tests/NC_002333.2.sam -g mudskipper/tests/NC_002333.2.gtf -o demo.bam
[2024-04-18T16:07:56Z INFO mudskipper] Mudskipper started...
[2024-04-18T16:07:56Z INFO mudskipper::annotation] reading the gtf file and building the tree.
█████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████ 100/100
[2024-04-18T16:07:56Z INFO mudskipper::annotation] building the tree
[2024-04-18T16:07:56Z INFO mudskipper::annotation] Time to build the tree: 9.174297ms
[2024-04-18T16:07:56Z INFO mudskipper::bam] Number of reference sequences: 24
[2024-04-18T16:07:56Z INFO mudskipper::bam] thread count: 1
[2024-04-18T16:07:56Z INFO mudskipper] Mudskipper finished.
Regards,
Steven
The text was updated successfully, but these errors were encountered:
steven-bioinfo
changed the title
Number of reference sequences: 0
No mapped read in output bam
Apr 18, 2024
steven-bioinfo
changed the title
No mapped read in output bam
No mapped read in output bam with gff file
Apr 18, 2024
Hello,
I want to convert my bam file mapped on genomic assembly (bam by bowtie2) to my predicted genes (gff by prokka). Next I want to use this bam as salmon input.
But when i use the mudskipper bulk command line no reads are mapped on the transcriptome.
I thinks, they have a bug when the gff is loading because This information are logged :
2024-04-18T16:04:35Z INFO mudskipper::bam] Number of reference sequences: 0
I also test our gff in tests folder and i have the same results.
When i use the gtf test file reference sequence are detected and sequence are mapped:
Regards,
Steven
The text was updated successfully, but these errors were encountered: