Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,023,591 | 0 | T | G | 51.6% | 32.3 / 45.2 | 31 | V191G (GTG→GGG) | yedA | drug/metabolite exporter YedA |
Reads supporting (aligned to +/- strand): ref base T (4/11); new base G (16/0); total (20/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.61e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTTGGCTCAGTTTATGGCTCGCGCATTACCTTACCTGTAGGGATGATGGCGGGTGCGATTGAGATGCTGGCGGCAGGCGTGGTGTTAATGATCGCGTCGATGATTGCGGGTGAAAAACTGACGGCGCTCCCTTCCCTTTCAGGCTTCCTTGCGGTCGGCTATCTGGCGC > NZ_CP009273/2023509‑2023677 | tttGGCTCAGTTTATGGCTCGCGCATTACCTTACCTGTAGGGATGATGGCGGGGGCGGTTGAGGTGCTGGGGGGGGGCGGGGGGGTAATg > 2:792230/1‑90 (MQ=255) tGGCTCAGTTTATGGCTCGCGCATTACCTTACCTGTAGGGATGATGGCGGGGGGGGTTGAGGTGGTGGGGGGAGGGGGGGGGTTTATGAt > 2:535421/1‑90 (MQ=255) tCAGTTTATGGCTCGCGCATTACCTTACCTGTAGGGATGATGGCGGGTGCGATTGAGATGCTGGCGGCAGGCGTGGTGTTAATGATcgcg < 1:776012/90‑1 (MQ=255) tCGCGCATTACCTTCCCTGTAGGGATGATGGCGGGTGCGATTGAGATGCTGGCGGCAGGCGTGGTGTTAATGATCGCGTCGATGATTGCg < 2:12740/90‑1 (MQ=255) tCGCGCATTACCTTACCTGTAGGGATGATGGCGGGTGCGATTGAGATGCTGGCGGCAGGCGTGGTGTTAATGATCGCGTCGATGATTGCg < 2:705989/90‑1 (MQ=255) cgcgCATTACCTTACCTGTAGGGATGATGGCGGGGGGGGTTTTGAAGCTGGGGGGGGGGGGGGGGGTAATGAGCGCGGCGGGGAtttggg > 2:1888040/1‑86 (MQ=255) cgcgCATTACCTTACCTGTAGGGATGATGGCGGGGGCGGTTTAGATGCTGGGGGGAGGGGGGGGGGTAAATGTTGCGGCGATGAtttggg > 2:659951/1‑86 (MQ=255) cgcgCATTACCTTACCTGTAGGGATGATGGCGGGGGCGATTGAGATGCTGGCGGCAGGGGGGGGGGTAAAGATCGCGGCGGTGGTTGCgg > 2:120002/1‑90 (MQ=255) taccttacctGTAGGGATGATGGCGGGGGGGGGTGAGATGCTGGGGGCCGGCGGGGGGGTAAAGATCGCGGCGGTGGTTGGGGGGGgaaa > 1:2023405/1‑90 (MQ=255) ttacctGTAGGGATGATGGCGGGTGCGATTGAGATGCTGGCGGCAGGCGTGGGGTTAATGATCGCGTCGATGATTGCGGGTGAAAAACTg > 2:1371930/1‑90 (MQ=255) tacctGTAGGGATGATGGCGGGTGCGATTGAGATGCTGGCGGCAGGCGTGGTGTTAATGATCGCGTCGATGATTGCGGGTGAAAAACTGa < 2:1695677/90‑1 (MQ=255) ctGTAGGGATGATGGCGGGTGCGATTGAGATGCTGGCGGCAGGCGTGGTGTTAATGATCGCGTCGATGATTGCGGGTGAAAAACTGACgg < 1:1103984/90‑1 (MQ=255) tGTAGGGATGATGGCGGGGGCGATTGAGATGCTGGCGGCAGGCGTGGGGTTAATGATCGCGTCGGTGATTGCGGGGGAAAAACTGACGgc > 1:632506/1‑90 (MQ=255) tAGGGATGATGGCGGGTGCGATTGAGATGCTGGCGGCAGGCGTGGTGTTAATGATCGCGTCGATGATTGCGGGTGAAAAACTGACGgcgc > 2:466164/1‑90 (MQ=255) tAGGGATGATGGCGGGTGCGATTGAGATGCTGGCGGCAGGCGGGGGGTTAATGATCGCGGCGGTGATTTCGGGGGGAAAACTGACGgcgc > 1:1894046/1‑90 (MQ=255) tAGGGATGATGGCGGGGGCGATTGAGATGCTGGCGGCAGGCGTGGGGTTAATGATAGCGTCGGTGATTGCGGGGGAAAAACTGACGgggc > 1:1549416/1‑90 (MQ=255) tAGGGATGATGGCGGGGGCGATTGAGATGCTGGCGGCAGGCGTGGGGGTAATGATCGCGTCGGTGATTGCGGGGGGAAAACTGACGgggc > 1:7297/1‑90 (MQ=255) aGGGATGATGGCGGGTGCGATTGAGATGCTGGCGGCAGGCGTGGGGTTAATGATCGCGTCGATGATTGCGGGGGAAAAACTGACGGCGCt > 2:1017738/1‑90 (MQ=255) atgatgGCGGGGGCGATTGAGATGCTGGCGGCAGGCGGGGGGTTAATGATCGCGTCGGTGATTGGGGGGGAAAAACTGACGGCGCtccct > 1:2070375/1‑90 (MQ=255) tgatgGCGGGTGCGATTGAGATGCTGGCGGCAGGCGTGGTGTTAATGATCGCGTCGATGATTGCGGGTGAAAAACTGACGGCGCtccctt < 2:1300831/90‑1 (MQ=255) tgatgGCGGGTGCGATTGAGATGCTGGCGGCAGGCGTGGTGTTAATGATCGCGTCGATGATTGCGGGTGAAAAACTGACGGCGCtccctt > 2:805446/1‑90 (MQ=255) tgatgGCGGGTGCGATTGAGATGCTGGCGGCAGGCGGGGGGGTAATGATCGCGTCGATGATTGCGGGGGAAAAAATGACGGCGCtccctt > 1:456874/1‑90 (MQ=255) tgatgGCGGGGGCGATTGAGATGCTGGCGGCAGGCGTGGGGTTAATGATCGCGGCGGTGATTGGGGGGGAAAAAATGACGGGGCtccctt > 1:1734157/1‑90 (MQ=255) tgGCGGGTGCGATTGAGATGCTGGCGGCAGGCGTGGGGTTAATGATCGCGGCGATGATTGGGGGGGAAAAACTGACGGCGCtcccttccc > 2:1192092/1‑90 (MQ=255) gCGATTGAGATGCTGGCGGCAGGCGTGGTGTTAATGATCGCGTCGATGATTGCGGGTGAAAAACTGACGGCGCTCCCTTCCCTTTCAGGc < 2:1894046/90‑1 (MQ=255) gCGATTGAGATGCTGGCGGCAGGCGTGGTGTTAATGATCGCGTCGATGATTGCGGGTGAAAAACTGACGGCGCTCCCTTCCCTTTCAGGc < 1:1302626/90‑1 (MQ=255) tGCTGGCGGCAGGCGTGGTGTTAATGATCGCGTCGATGATTGCGGGTGAAAAACTGACGGCGCTCCCTTCCCTTTCAGGCTTCCTTGCgg < 2:632506/90‑1 (MQ=255) cTGGCGGCAGGCGTGGTGTTAATGATCGCGTCGATGATTGCGGGTGAAAAACTGACGGCGCTCCCTTCCCTTTCAGGCTTCCTTGCGGTc < 2:1757143/90‑1 (MQ=255) tGGCGGCAGGCGTGGTGTTAATGATCGCGTCGATGATTGCGGGTGAAAAACTGACGGCGCTCCCTTCCCTTTCAGGCTTCCTTGCGGTCg > 2:1818572/1‑90 (MQ=255) gcggcAGGCGTGGTGTTAATGATCGCGTCGATGATTGCGGGTGAAAAACTGACGGCGCTCCCTTCCCTTTCAGGCTTCCTTGCGGTCGGc < 2:7297/90‑1 (MQ=255) gcAGGCGTGGTGTTAATGATCGCGTCGATGATTGCGGGTGAAAAACTGACGGCGCTCCCTTCCCTTTCAGGCTTCCTTGCGGTCGGCTAt > 2:72909/1‑90 (MQ=255) tggtgTTAATGATCGCGTCGATGATTGCGGGTGAAAAACTGACGGCGCTCCCTTCCCTTTCAGGCTTCCTTGCGGTCGGCTATCTGgcgc < 1:1469468/90‑1 (MQ=255) | TTTGGCTCAGTTTATGGCTCGCGCATTACCTTACCTGTAGGGATGATGGCGGGTGCGATTGAGATGCTGGCGGCAGGCGTGGTGTTAATGATCGCGTCGATGATTGCGGGTGAAAAACTGACGGCGCTCCCTTCCCTTTCAGGCTTCCTTGCGGTCGGCTATCTGGCGC > NZ_CP009273/2023509‑2023677 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |