| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,089,571 | 0 | A | G | 46.0% | 66.2 / 16.1 | 50 | I83V (ATT→GTT) | hisF | imidazole glycerol phosphate synthase subunit HisF |
| Reads supporting (aligned to +/- strand): ref base A (1/26); new base G (23/0); total (24/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.22e-13 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGTCTC‑TGGAAGAT‑GCCGCGAAAATTCTTTCCTTTG‑GCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACAT > NZ_CP009273/2089485‑2089658 | cTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTa < 1:1554074/90‑1 (MQ=255)cTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTa < 1:1565469/90‑1 (MQ=255) tCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAg < 2:2184484/90‑1 (MQ=255) cAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGtc < 1:157887/90‑1 (MQ=255) cAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGtc < 1:2024649/90‑1 (MQ=255) gCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGtctc < 2:329809/90‑1 (MQ=255) cGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGTTAAATctc‑t > 1:474722/1‑90 (MQ=255) aTGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGGTTAAGTCTC‑Tgg > 2:1869906/1‑90 (MQ=255) cGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGGTTAGGCTC‑TTGAAGAt > 2:1528638/1‑90 (MQ=255) gtgtGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGGTAAGGGTC‑TGGGAGAT‑g > 2:967916/1‑90 (MQ=255) tgtgGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGTTAAAGCTC‑TTGGAGAT‑Gc > 2:444717/1‑90 (MQ=255) aGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGTCTC‑TGGAAGAT‑GCCGCGaa < 1:1406116/90‑1 (MQ=255) aGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGTCTC‑TGGAAGAT‑GCCGCGaa < 2:1359752/90‑1 (MQ=255) aGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGTTAAGTCTC‑TTGTAGAT‑GCCGCGaa > 1:413217/1‑90 (MQ=255) aGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGTTAAGTCTC‑TGGGAGAT‑GCCGCGaa > 1:1128974/1‑90 (MQ=255) gATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGTTAAATTTC‑TGGGAGAA‑GCCCCGaaa > 1:941976/1‑90 (MQ=255) tAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGTTAAGGCTC‑TGGAAGAG‑GCCGCCaaaaa > 1:1451101/1‑89 (MQ=255) aaaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGTCTC‑TGGAAGAT‑GCCGCGAAAATTc < 1:2151320/90‑1 (MQ=255) aaaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGTTAAGTCTC‑TGGAAGAT‑GCCGCGAAAAttt > 1:794169/1‑89 (MQ=255) aaaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGTTAAGTCTC‑TGGAAGAT‑GCCCCGaaaaatt > 1:2255048/1‑87 (MQ=255) aGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGTCTC‑TGGAAGAT‑GCCGCGAAAATTCtt < 1:564654/90‑1 (MQ=255) gCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGTTAAGGCTC‑TGGAAGAT‑GGCGCGGAAAtttttt > 2:1530595/1‑90 (MQ=255) gggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGTTAAATCTC‑TGGGAGATGGCCCGGAAAATTTTTTct > 2:696982/1‑89 (MQ=255) tctcGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGTTAAATCTC‑TGGAAGAT‑GCCGCGAAAATTCTTTTCTTTT‑gc > 1:1418644/1‑90 (MQ=255) tcGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGTTAAGTCTC‑TGGGAGAA‑GCCCCGGAAAATTTTTTCTTTT‑gcgc > 2:608684/1‑90 (MQ=255) cgcgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGGTTAAGtc > 1:1113742/1‑53 (MQ=255) cgcgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGtc < 2:1113742/53‑1 (MQ=255) cgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGGGTTAGTGTC‑TGGGAGGT‑GGCCCGGAAATTCTTTTCTTTT‑GGGCggg > 2:1034092/1‑89 (MQ=255) cgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGTTAAGTTTC‑TGGGAGAT‑GCCGCGAAAATTTTTTTCTTTT‑GGGCggg > 1:954922/1‑89 (MQ=255) gATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTTAAGTTTT‑GGGAGAAG‑GCCCGGAAAATTTTTTCTTTTGCGCGGGGA‑AAAATTTcat > 2:173420/1‑88 (MQ=255) gATCGACATTCCGTTTTGTGTGGC‑GGGGGGGGTTAAATGTCTTGGAAAAG‑GCC‑CGAAAAATTTTTTCTTTT‑TCGCGGGTAAAAATTTCAt > 2:1607652/1‑90 (MQ=255) aTCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGTCTC‑TGGAAGAT‑GCCGCGAAAATTCTTTCCTTTG‑GCGCGGATAAAATTTCCATc < 1:696982/90‑1 (MQ=255) gACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGTCTC‑TGGAAGAT‑GCCGCGAAAATTCTTTCCTTTG‑GCGCGGATAAAATTTCCATCAAc < 2:1451101/90‑1 (MQ=255) gACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGTCTC‑TGGAAGAT‑GCCGCGAAAATTCTTTCCTTTG‑GCGCGGATAAAATTTCCATCAAc < 1:925713/90‑1 (MQ=255) gACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGTCTC‑TGGAAGAT‑GCCGCGAAAATTCTTTCCTTTG‑GCGCGGATAAAATTTCCATCAAc < 1:222051/90‑1 (MQ=255) gACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGTCTC‑TGGAAGAT‑GCCGCGAAAATTCTTTCCTTTG‑GCGCGGATAAAATTTCCATCAAc < 1:323817/90‑1 (MQ=255) gACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGTCTC‑TGGAAGAT‑GCCGCGAAAATTCTTTCCTTTG‑GCGCGGATAAAATTTCCATCAAc < 2:2273137/90‑1 (MQ=255) cATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGTCTC‑TGGAAGAT‑GCCGCGAAAATTCTTTCCTTTG‑GCGCGGATAAAATTTCCATCAActc < 1:1389949/90‑1 (MQ=255) ttCCGTTTTGTGTGGC‑GGGGGGGGGTAAGTCTC‑TGGGAGAT‑GCCGCGAAAAATTTTTTCTTTT‑GCGCGGGTAAAATTTTCATCAACtctc > 2:1888145/1‑90 (MQ=255) ttGTGTGGC‑GGGGGGGGTTAAGTCTC‑TGGAAGAA‑GCCGCGAAAAATTTTTTCTTTT‑GCGGGGGTAAAAATTTCATCAAATCTCCTGCGCt > 2:1152300/1‑90 (MQ=255) gtgtgGC‑GGGGGGGGTTAAATCTC‑TTGGAGAT‑GCCGCGAAAATTTTTTTCTTTG‑GGGCGGGTAAAATTTTCCTCAACTCTCCTTCGCTgg > 2:2126787/1‑90 (MQ=255) tgtgGC‑GGGTGGGATTAAGTCTC‑TGGAAGAT‑GCCGCGAAAATTCTTTCCTTTG‑GCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGc < 2:474722/90‑1 (MQ=255) tgtgGC‑GGGTGGGATTAAGTCTC‑TGGAAGAT‑GCCGCGAAAATTCTTTCCTTTG‑GCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGc < 2:941976/90‑1 (MQ=255) gtgGC‑GGGTGGGATTAAGTCTC‑TGGAAGAT‑GCCGCGAAAATTCTTTCCTTTG‑GCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCa < 2:1229757/90‑1 (MQ=255) gtgGC‑GGGTGGGATTAAGTCTC‑TGGAAGAT‑GCCGCGAAAATTCTTTCCTTTG‑GCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCa < 2:957125/90‑1 (MQ=255) gggTGGGATTAAGTCTC‑TGGAAGAT‑GCCGCGAAAATTCTTTCCTTTG‑GCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGAccc < 1:2166962/90‑1 (MQ=255) gggTGGGATTAAGTCTC‑TGGAAGAT‑GCCGCGAAAATTCTTTCCTTTG‑GCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGAccc < 1:1659303/90‑1 (MQ=255) ggTGGGATTAAGTCTC‑TGGTAGAT‑GCCGCGAAAATTCTTTCCTTTG‑GCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCa < 2:413217/90‑1 (MQ=255) gTGGGATTAAGTCTC‑TGGAAGAT‑GCCGCGAAAATTCTTTCCTTTG‑GCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCaa < 1:1006363/90‑1 (MQ=255) ggATTAAGTCTC‑TGGAAGAT‑GCCGCGAAAATTCTTTCCTTTG‑GCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACat > 2:924725/1‑90 (MQ=255) | CTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGC‑GGGTGGGATTAAGTCTC‑TGGAAGAT‑GCCGCGAAAATTCTTTCCTTTG‑GCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACAT > NZ_CP009273/2089485‑2089658 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |