Comments (4)
After samtools view -F 4 ERR067664.sam | awk '$3=="NC_000962.3"' | wc -l
is executed, the result is 2206287, which is consistent with the number of TB reads in the tsv file !
Thank you very much for your patient answer!
from clockwork.
Not sure right now why the difference of 1065. Clockwork just counts the number of reads matched to H37Rv, and that's what it reports in TB
.
Your second question: this will be because of read pairs. If a read matches H37Rv, then the both reads of the pair are kept regardless of what the read's mate did (it could be unmapped or mapped to a contamination sequence). Either both reads of a pair are kept, or they are both removed. The surviving read pairs are what is written to the output fastq files - the total number of which will be the count in Reads_kept_after_remove_contam
.
from clockwork.
Is my method of using samtools view -F 4 ERR067664.sam | grep NC_000962.3 | wc -l
to obtain the number of reads matched to H37Rv consistent with your method? Is there something wrong with my statistical method?
Regarding the second question, I understand, thank you for your answer!
from clockwork.
grep NC_000962.3
will be looking for NC_000962.3
anywhere in the line, not just in column 3. This means it counts a read that is mapped to anything, if its mate is mapped to NC_000962.3
.
Try awk '$3=="NC_000962.3"'
instead of grep NC_000962.3
from clockwork.
Related Issues (20)
- HMPDACC downloads not working HOT 2
- gnuplot error libQt5Core.so.5: cannot open shared object file: No such file or directory HOT 1
- ERROR : Unknown image format/type HOT 4
- QC pipeline with script only HOT 4
- approach to classify mutations HOT 5
- Missed INH-R caused by deletion of the entire katG gene HOT 6
- maximize cores and threads when running pipeline HOT 7
- the minimum read depth for variant calling HOT 3
- Update nextflow to scripts to use dsl2
- enaDataGet not working: HTTP Error 400 HOT 2
- Samples with a small number of reads causes Clockwork to crash HOT 1
- AssertionError: When I run clockwork remove_contam command HOT 5
- Error running this command: /bioinf-tools/cortex/scripts/calling/run_calls.pl HOT 11
- Why map to the H37Rv reference genome twice using minimap? HOT 6
- Error running this command: fqtools count HOT 7
- Citation HOT 2
- Error when running qc.nf HOT 4
- Error connecting to database HOT 1
- I don't know what's wrong HOT 3
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 clockwork.