Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,319,929 | 0 | T | G | 55.0% | 1.5 / 22.9 | 20 | G114G (GGT→GGG) | atoB | acetyl‑CoA acetyltransferase |
Reads supporting (aligned to +/- strand): ref base T (1/8); new base G (11/0); total (12/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.14e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGA > NZ_CP009273/2319851‑2320012 | tgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAatat > 1:24729/1‑90 (MQ=255) tCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATTAgg > 2:292056/1‑89 (MQ=255) cGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATGAggt > 1:310660/1‑88 (MQ=255) cGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAggt > 1:173691/1‑88 (MQ=255) gggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGttt < 2:706686/90‑1 (MQ=255) ttAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGccc < 2:682771/90‑1 (MQ=255) ttAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGccc < 2:24729/90‑1 (MQ=255) gtgtGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTTGCCCCCcact > 1:525729/1‑90 (MQ=255) gtgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTTACCCCCtactta > 1:324375/1‑90 (MQ=255) gcgcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTACGGAAAATATGAGTTTAGCCCCCTACTTACTc < 1:532323/90‑1 (MQ=255) gcgcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGGTTTGCCCCCCACCTACTc > 2:145877/1‑90 (MQ=255) ttCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAATATGAGGTTAGCCCCCTCCTTACTCGATGCAAAAGCACGCTCTg > 2:158143/1‑90 (MQ=255) gcaggTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTAGCCCCCCACTTACTCGATGCAAAAACACGCCCTGGTTAt > 2:214108/1‑90 (MQ=255) gTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTc > 1:179398/1‑90 (MQ=255) gCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGAc < 1:621144/90‑1 (MQ=255) aTTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAg < 1:292056/90‑1 (MQ=255) gCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTAt < 2:475070/90‑1 (MQ=255) ggggggGATGGGAAATATGAGTTTAGCCCCCCACTTACTCGATGCAAAAGCACGCCCTGGTTATCGTCTTGGAGGCGGGCAGGGTTATGa > 1:475070/1‑90 (MQ=255) gTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCaa < 2:391417/42‑1 (MQ=255) gTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCaa > 1:391417/1‑42 (MQ=255) | TGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGA > NZ_CP009273/2319851‑2320012 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |