Comments (5)
Hello,
Looking more closely to this issue, it appears that the bam file produced by longshot also contains supplementary and secondary alignments, even though these are filtered out before SNV discovery. Filtering these out, I get 29240 unphased reads, but I should be getting 26476, so there are ~ 3 thousand additional reads. Wat else could I be missing?
Cheers,
Juan D. Montenegro
from longshot.
Longshot does output all alignments therefore the total number of reads in the output should be exactly the same as in the input bam. The statistics are only for the filtered reads. Can you confirm if the extra reads are duplicates?
from longshot.
Hello,
Sorry for the late reply. The input and output do have the same number of records. So I am probably missing some filter to make the numbers match. Besides A30 and secondary/supplementary alignments, what else is being ignored by Longshot? I amasking this, beacuse I am usually, reassemblying the phased and unphased reads. Usually what I see is that one end of the original contig cannot be phased, while the other end is split into two homologous haplotypes, but in some cases I do get a very complex partition after assemblying, so I guess some of the unphased reads actually help complete gaps between haplotypes but were too short and did not contain enough polymorphisms to be phased.
Any help selecting the appropriate set of unphased reads would be very helpful.
Cheers,
Juan D.
from longshot.
Longshot filters out reads with low mapping quality in addition to secondary/supp. alignments. I have copied the list of filters from the code below:
record.is_quality_check_failed()
|| record.is_duplicate()
|| record.is_secondary()
|| record.is_unmapped()
|| record.mapq() < min_mapq
|| record.is_supplementary()
All these reads will be output as 'unphased'.
from longshot.
from longshot.
Related Issues (20)
- Question about recommendation HOT 1
- Problem when running Longshot with too few reads/too small coverage using -A HOT 2
- Error on --version
- Using BED file as regions HOT 1
- Number of True Variants in HG001 VCF is too low HOT 1
- can we use longshot to test haploid SNP, even organelles HOT 4
- incorrectly calling homozygous genotype with low alternate alleles HOT 3
- Parameter for multi-threading? HOT 2
- Availability through conda? HOT 6
- Advice for phasing reads of a triploid plant species?
- outputting haplotagged bam fails if input bam doesn't contain HP tags HOT 3
- Homozygous sites HOT 1
- call variants with bed file HOT 1
- Potential vcf file format HOT 1
- Does it support multithreading? HOT 1
- joint variant calling HOT 1
- VCF output file format (ALT allele frequency)
- failed to create pileup HOT 1
- Memory usage on PromethION genomes HOT 1
- thread 'main' panicked at 'index out of bounds
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from longshot.