Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,089,571 | 0 | A | G | 53.3% | 9.7 / 31.1 | 30 | I83V (ATT→GTT) | hisF | imidazole glycerol phosphate synthase subunit HisF |
Reads supporting (aligned to +/- strand): ref base A (2/12); new base G (16/0); total (18/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.05e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.37e-05 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGTCTCTGGAAGATGCCGC‑G‑AAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAA > NZ_CP009273/2089486‑2089655 | ttCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTaa < 2:1383145/90‑1 (MQ=255) tCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAg < 2:747425/90‑1 (MQ=255) ccAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGGTTa > 1:1240247/1‑87 (MQ=255) ccAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTa < 2:1240247/87‑1 (MQ=255) ccAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGt < 2:868308/90‑1 (MQ=255) cAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGtc < 1:256461/90‑1 (MQ=255) cGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGTTAAATctct > 2:803117/1‑90 (MQ=255) tGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGATTAATGCTCTggg > 1:848957/1‑89 (MQ=255) tgtgGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGTTAAGTCTCTTGGAGAAGc > 1:1434232/1‑90 (MQ=255) tgtgGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGGGGC‑GGGGGGGGTTAAGGGTCTGGGAggggc > 1:907134/1‑86 (MQ=255) gTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGGTGAGTCTCTTGGAGATGCcgc‑g > 1:1193568/1‑90 (MQ=255) gTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGGTAAGGCTCTGGAAGAAGCcgc‑g > 1:633311/1‑90 (MQ=255) aaaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGTTAAATTTCTGGGAGGAGCCCC‑G‑AAAAttt > 1:971619/1‑89 (MQ=255) aaaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGGTAAGTCTCTGGGAGATGGCGC‑G‑GAAAttt > 1:587178/1‑89 (MQ=255) aGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGGTTAAGGCTCTTGAAGGTGCCGC‑G‑AAAAATCtt > 2:1470111/1‑90 (MQ=255) aGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGGTTAAATCTCTGGTAAACGCCGC‑C‑AAAAtttct > 2:423215/1‑87 (MQ=255) gggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGTTAAGAGTCTGGGAGAATCCGCGG‑AAAATTTTTTct > 1:773165/1‑89 (MQ=255) ggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTTAGGTTTTTGGAGGTGGCCCG‑G‑AAAATTTTTTctt > 1:116183/1‑88 (MQ=255) ggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGTTTAGGTTCTGGGAGAAGCCGC‑GCAAAATTTTTTctt > 1:89640/1‑88 (MQ=255) gcgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGTTAAATCTCTGGAAGAAGCCGC‑G‑AAAAATTTTTTCTTTTGGGCgg > 2:303762/1‑90 (MQ=255) tGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGTTAAGTCTCTGGGAGATGCCCC‑G‑AAAAATTTTTTCTTTTGCGCGGGTa > 2:677336/1‑90 (MQ=255) ggCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGTCTCTGGAAGATGCCGC‑G‑AAAATTCTTTCCTTTGGCGCGGATaa < 2:971619/90‑1 (MQ=255) aCATTCCGTTTTGTGTGGC‑GGGGGGGGTTAAGTCTCTGGAAGATGCCGC‑G‑AAAATTTTTTTCTTTTGGGCGGGTAAAAATTTCATCAAct > 2:470714/1‑90 (MQ=255) ttCCGTTTTGTGTGGC‑GGGTGGGATTAAGTCTCTGGAAGATACCGC‑G‑AAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACtctc < 1:803117/90‑1 (MQ=255) gTTTTGTGTGGC‑GGGTGGGATTAAGTCTCTGGAAGATGCCGC‑G‑AAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTgc < 1:470714/90‑1 (MQ=255) tttGTGTGGC‑GGGTGGGATTAAGTCTCTGGAAGATGCCGC‑G‑AAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTgcgc < 1:1470111/90‑1 (MQ=255) gtgtgGC‑GGGTGGGATTTAATCTCTGGGAGAAGCCGC‑G‑AAAATTCTTTTCTTTTGCGCGGGAAAAATTTCCATCAACCCTCCTGCGCTgg > 2:810566/1‑90 (MQ=255) gC‑GGGTGGGATTAAGTCTCTGGAAGATGCCGC‑G‑AAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGAc < 2:1117441/90‑1 (MQ=255) ggTGGGATTAAGTCTCTGGAAGATGCCGC‑G‑AAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCa < 1:1092006/90‑1 (MQ=255) gTGGGATTAAGTCTCTGGAAGATGCCGC‑G‑AAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCaa < 2:388046/90‑1 (MQ=255) | TTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGTCTCTGGAAGATGCCGC‑G‑AAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAA > NZ_CP009273/2089486‑2089655 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |