| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,039,701 | 0 | A | C | 40.0% | 8.8 / 13.9 | 15 | N75T (AAC→ACC) | msrP | periplasmic protein‑L‑methionine sulfoxide reductase catalytic subunit |
| Reads supporting (aligned to +/- strand): ref base A (9/0); new base C (0/6); total (9/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.00e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGCGCTGGAGTTCAGCAAG‑CCT‑GCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAACTTCTATGAATTCGGGCTGGAT > NC_000913/2039663‑2039765 | ggcgcTGGATTTCAGCAAG‑CCT‑GCCGCCTGGAAAAATACCCTGCCACTGACGCCAGCAGATAAAg < 1:3353088/64‑1 (MQ=255)cgcgcTGGATTTCACAAAC‑CCT‑GCCGCCGGGAAAAAAACCCTGCCACTGACGCCAGCAGATAAAg < 1:3381373/64‑1 (MQ=255) gcTGGAGTTCAGCAAG‑CCT‑GCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGttt > 1:540163/1‑65 (MQ=255) ggAGTTC‑GCAAGCCCT‑CCCTCGTGGGAAACAACCCTGCCACTGACGCCAGCAGATAAAGTTTCCg < 1:1676863/65‑1 (MQ=255) aGTTCAGCAAG‑CCT‑GCCGCCTGGCAAAATACCCTGCCACTGACGCCAGCAGATAAAGTTTCCGGt < 1:3659416/65‑1 (MQ=255) aGTTCAGCAAG‑CCT‑GCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGt > 1:957612/1‑65 (MQ=255) tCAGCAGG‑CCG‑GCCGCCTGGCAAAATACCCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTtat < 1:822315/65‑1 (MQ=255) tCAGCAAG‑CCT‑GCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTtat > 1:1111098/1‑65 (MQ=255) tCAGCAAG‑CCG‑GCCGCCTGGCAAAATACCCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTtat < 1:514459/65‑1 (MQ=255) cgCAAG‑CCCGGCCGCCTGAAAAAATACCCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTtata < 1:295081/64‑1 (MQ=255) aaG‑CCT‑GCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTAtaataa > 1:2281437/1‑65 (MQ=255) gg‑CCT‑GCCGCCGGAAAAAATACCCTGCCGCTGACGCCAGCAGATAAAGTTTCCGGTTATAATAAc < 1:1014893/64‑1 (MQ=255) g‑CCT‑GCCGCCTGGCAAAATACCCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAACt < 1:98196/65‑1 (MQ=255) g‑CCT‑GCCGCCGGGCAAAATACCCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAACt < 1:813945/65‑1 (MQ=255) gggt‑GCCGCCTGGCAAAATACCCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAACt < 1:3838238/62‑1 (MQ=255) ccT‑GCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAACtt > 1:2392081/1‑65 (MQ=255) t‑GCCGCCTGGCAAAATACCCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAACTTCt < 1:2448329/65‑1 (MQ=255) ggccgccTGGAAAAATACCCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAACTTCt < 1:3492837/64‑1 (MQ=255) ccTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAACTTCTATGaa > 1:1086547/1‑65 (MQ=255) cAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAACTTCTATGAATTCgg > 1:2469581/1‑65 (MQ=255) cAAAAAACCCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAACTTCTATGAATTCgg < 1:1330167/65‑1 (MQ=255) aTAACCTGCCACTGACGCCAGCAGCTAAAGTTTCCGGTTACAATAACTTCTATGAATTCGGGCTg > 1:2120667/1‑65 (MQ=255) tAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAACTTCTATGAATTCGGGCTgg > 1:2260103/1‑65 (MQ=255) aCCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAACTTCTATGAATTCGGGCTGGAt > 1:2576322/1‑65 (MQ=255) | AGCGCTGGAGTTCAGCAAG‑CCT‑GCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAACTTCTATGAATTCGGGCTGGAT > NC_000913/2039663‑2039765 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |