Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 816,005 | 0 | G | T | 61.5% | 4.9 / 16.9 | 13 | L222L (CTG→CTT) | ybhL | Bax inhibitor‑1/YccA family protein |
Reads supporting (aligned to +/- strand): ref base G (5/0); new base T (0/8); total (5/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.77e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTGGCGCGTTAACCTTGTATCTGGACTTCATCAACCTGTTCCTGATGTTGTTGCGGATCTTCGGCAACCGCCGTTAATTCTTTTTCACCATCACTATCCCGGCCATGTGTCGGGATTTTTTTCTTACTAAAAACGACACTTAACAAGTACCAGGTATAAAA > NZ_CP009273/815912‑816128 | tatGGGTGAGTATGTTCACCCCCGCGACCCGTCAACCCTCCAAAAATTCCCCATCTTTGGCTCTTTACCCTTTTATCTGGTCTCCACCAACCTTTTCTTGTTGTTGTGGCGGTTCTCCGCAAACCGCCGTTAATTCttt < 2:406857/139‑1 (MQ=255) tGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTGGCGCGTTAACCTTGTATCTGGACTTCATCAACCTGTTCCTGATGTTGTTGCGGATCTTCGGCAACCGCCGTTAATTCttttt > 1:197854/1‑139 (MQ=255) gggggCAGATTGATCCCCGCGACACTTCGACCCTCCGAAAATTTCCCTTTGTTGGCGCTTTACCCTTGTATCTGGACTTCATCACCTTTTTCTTGTTTTTGTGGCGGTTCTTCGCCAACCGCCGTTAATTCTTTTTCAc < 2:197854/135‑1 (MQ=255) gTGAGTAGTTTGATCCCCGCGACTCGTCGACCCTGCAAAAATTTCCCTTTCTGGGCGCTTTACCTTTGTATCTGGACTCCATCATCTTTTTCTTGTTTTTGTTGCGGTTCTTCGGAAACCGCCGTTAATTCTTTTTCAc < 2:156506/139‑1 (MQ=255) gTGAGTAGTTTGACCCCCGCGACACGTCGACCCTGCAAAATTATTCCTTTCTTGCCGTTTTACCTTTTTATCTGCCCTTCATCACCTTTTTCTTTATGTTGTTGCGTATCTGCGGCAACCGCCGTTAATTCTTTTTCAc < 2:548100/139‑1 (MQ=255) gTGAGCAGTTTGATCCCCGCGACACGTCGACCCTGCGAAAATATTCCATTCTTGGCGCTTAATCCTTTTTTCTGGTCTTCACCACCCTTTTCTTGTTGTTTTTGCGGATCTTCGGAAACCGCCGTTAATTCTTTTTCAc < 1:585667/139‑1 (MQ=255) aGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTGGCGCGTTAACCTTGTATCTGGACTTCATCAACCTGTTCCTGATGTTGTTGCGGATCTTCGGCAACCGCCGTTAATTCTTTTTCACCAt > 1:53597/1‑139 (MQ=255) tataTTGACCCCCGCCACTCTTCCAACCTGCACAAATATTCCTTTCTGGGCTCTTTACCCTTTTATCTTCACTCCACCACCTTTTTCTTGTTGTTTTTGCGACCCTCCGGAAACCGCCGTTAATTCTTTTTCACCATCa < 2:23714/136‑1 (MQ=255) gcgACACGTCGAACCTGCGCAAATATTCCATTCTTGGCGCGTTAACCTTGTATCTGGACTTCATCAACCTGTTCCTGATGTTGTTGCGGATCTTCGGCAACCGCCGTTAATTCTTTTTCACCATCACTATCCCGGCCAt > 2:494552/1‑139 (MQ=255) acctttGTTTCTGCACTTCATCACCTTTTTCTTTTTTTTGTTGCGGTTCTGCGGCAACCGCCGTTAATTCTTTTTCACCATCACTATCCCGGCCATGTGTCGGGATTTTTTTCTTACTAAAAACGACACTTAACAAGt < 1:354423/134‑1 (MQ=255) cctttttATCTGAACTTCATCACCTTTTTCCTGTTGTTGTGGCGCATCTCCGGCAACCGCCGTTAATTCTTTTTCACCATCACTATCCCGGCCATGTGTCGGGATTTTTTTCTTACTAAAAACGACACTTAACAAGTAc < 2:53597/133‑1 (MQ=255) ggACTTCATCAACCTGTTCCTGATGTTGTTGCGGATCTTCGGCAACCGCCGTTAATTCTTTTTCACCATCACTATCCCGGCCATGTGTCGGGATTTTTTTCTTACTAAAAACGACACTTAACAAGTACCAGGTATaaaa > 1:592041/1‑139 (MQ=255) gACTTCATCAACCTGTTCCTGATGTTGTTGCGGATCTTCGGCAACCGCCGTTAATTCTTTTTCACCATCACTATCCCGGCCATGTGTCGGGATTTTTTTCTTACTAAAAACGACACTTAACAAGTACCAGGTATaaaa > 1:136340/1‑138 (MQ=255) | TATGGGTGAGCAGATTGATACCCGCGACACGTCGAACCTGCGCAAATATTCCATTCTTGGCGCGTTAACCTTGTATCTGGACTTCATCAACCTGTTCCTGATGTTGTTGCGGATCTTCGGCAACCGCCGTTAATTCTTTTTCACCATCACTATCCCGGCCATGTGTCGGGATTTTTTTCTTACTAAAAACGACACTTAACAAGTACCAGGTATAAAA > NZ_CP009273/815912‑816128 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |