Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,319,929 | 0 | T | G | 52.4% | 17.0 / 31.7 | 21 | G114G (GGT→GGG) | atoB | acetyl‑CoA acetyltransferase |
Reads supporting (aligned to +/- strand): ref base T (0/10); new base G (11/0); total (11/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.84e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.82e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGT > NZ_CP009273/2319842‑2320015 | tAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTAt < 2:149935/90‑1 (MQ=255) tATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGgaa > 2:279884/1‑90 (MQ=255) tATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGgaa > 2:226284/1‑90 (MQ=255) tgtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGgaaaa > 2:45174/1‑89 (MQ=255) tgtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAat > 2:71073/1‑90 (MQ=255) ggTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATg < 1:572768/90‑1 (MQ=255) gTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGa > 2:245121/1‑90 (MQ=255) gggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGttt < 1:584783/90‑1 (MQ=255) tgtgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAGCCCCCtactt > 2:11195/1‑90 (MQ=255) gtgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATgg < 1:429610/64‑1 (MQ=255) gtgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATgg > 2:429610/1‑64 (MQ=255) tGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGc > 2:295835/1‑90 (MQ=255) cccAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAg < 1:679046/90‑1 (MQ=255) aGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCAc < 1:434581/90‑1 (MQ=255) cATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGctc < 1:45174/90‑1 (MQ=255) gcaggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTAt < 1:295835/90‑1 (MQ=255) gCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGAc < 1:71073/90‑1 (MQ=255) aTTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCCACTTAATCGATGCAAAAACACGCCCTGGTTATTGTCTTGGGGACGGACAg > 2:436788/1‑90 (MQ=255) ggggggGATGGGAAATGTGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGGTATCGTCTTGGAGACGGGCAGGGTTATGa > 2:155924/1‑90 (MQ=255) ggggggGATGGGAAAAATGAGTTTAGCCCCCCACTTACTCGATGCAAAAACACGCCCTGGTTATCGTCTTGGGGACGGACAGGTTTATGa > 1:283908/1‑90 (MQ=255) gggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGt < 2:402727/90‑1 (MQ=255) | TAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGT > NZ_CP009273/2319842‑2320015 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |