| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 129,708 | 1 | . | T | 30.8% | 18.3 / 12.4 | 13 | intergenic (‑85/+300) | spa/sarH1 | immunoglobulin G binding protein A/staphylococcal accessory regulator H1 |
| Reads supporting (aligned to +/- strand): ref base . (6/3); new base T (2/2); total (8/5) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCTTTTTCAAATTAATACCCCCTGTATGTATTTGTAAAGTCATCATAATATAACGAATTATGTATTGCAATACTAAAATCTATATTTATAATTAAATTTAAAGGTAAGTTTTACAACTTATAAAATAAATA‑TTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGAGATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAA > CP000730/129578‑129841 | tctttttcAAATTAATACCCCCTGTATGTATTTGTAAAGTCATCATAATATAACGAATTATGTATTGCAATACTAAAATCTATATTTATAATTAAATTTAAAGGTAAGTTTTACAACTTATAAAATAAATATTTTTGCTGAAAGATTTAtt > 2:266570/1‑151 (MQ=255)tctttttcAAATTAATACCCCCTGTATGTATTTGTAAAGTCATCATAATATAACGAATTATGTATTGCAATACTAAAATCTATATTTATAATTAAATTTAAAGGTAAGTTTTACAACTTATAAAATAAATATTTTTGCTGAAAGATTTAtt > 2:266569/1‑151 (MQ=255) tACCCCCTGTATGTATTTGTAAAGTCATCATAATATAACGAATTATGTATTGCAATACTAAAATCTATATTTATAATTAAATTTAAAGGTAAGTTTTACAACTTATAAAATAAATA‑TTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCt < 1:582921/151‑1 (MQ=255) gtatTTGTAAAGTCATCATAATATAACGAATTATGTATTGCAATACTAAAATCTATATTTATAATTAAATTTAAAGGTAAGTTTTACAACTTATAAAATAAATA‑TTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCaaa > 2:343203/1‑151 (MQ=255) catcatAATATAACGAATTATGTATTGCAATACTAAAATCTATATTTATAATTAAATTTAAAGGTAAGTTTTACAACTTATAAAATAAATA‑TTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTc > 2:502139/1‑151 (MQ=255) aTTATGTATTGCAATACTAAAATCTATATTTATAATTAAATTTAAAGGTAAGTTTTACAACTTATAAAATAAATA‑TTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTGTTTCGATTAATAAAAACAGa < 2:218932/151‑1 (MQ=255) aTTGCAATACTAAAATCTATATTTATAATTAAATTTAAAGGTAAGTTTTACAACTTATAAAATAAATATTTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAat < 1:266570/151‑1 (MQ=255) aTTGCAATACTAAAATCTATATTTATAATTAAATTTAAAGGTAAGTTTTACAACTTATAAAATAAATATTTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAat < 1:266569/151‑1 (MQ=255) atttattttatAAGTTGTAAAACTTATAAAATAAATA‑TTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGAGATATTTTTTATAATGAAAGCGGTTTa > 1:651449/11‑151 (MQ=255) ggTAAGTTTTACAACTTATAAAATAAATA‑TTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTc < 1:402993/89‑1 (MQ=255) ggTAAGTTTTACAACTTATAAAATAAATA‑TTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTc > 2:402993/1‑89 (MQ=255) aCAACTTATAAAATAAATA‑TTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGAGATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTaaaaa > 1:223725/1‑151 (MQ=255) cAACTTATAAAATAAATA‑TTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGAGATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTaaaaaa > 2:1018754/1‑151 (MQ=255) | TCTTTTTCAAATTAATACCCCCTGTATGTATTTGTAAAGTCATCATAATATAACGAATTATGTATTGCAATACTAAAATCTATATTTATAATTAAATTTAAAGGTAAGTTTTACAACTTATAAAATAAATA‑TTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGAGATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAA > CP000730/129578‑129841 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |