Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,096,090 | 0 | A | G | 42.5% | 52.6 / 11.7 | 40 | I83V (ATT→GTT) | hisF | imidazole glycerol phosphate synthase, catalytic subunit with HisH |
Reads supporting (aligned to +/- strand): ref base A (4/19); new base G (17/0); total (21/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.94e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGT > NC_000913/2095950‑2096240 | cAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTctct > 1:526713/1‑151 (MQ=255) aaaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCCCCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTgg < 2:62669/151‑1 (MQ=255) aaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGa < 1:721105/151‑1 (MQ=255) cgctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTTAGTGTCTTGaa > 2:430622/1‑150 (MQ=255) gctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGAAGAT‑Gcc > 2:616417/1‑150 (MQ=255) tGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTAAGTCTCTTGAAAAT‑Gcccc > 2:93128/1‑148 (MQ=255) aagaagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATTAAGTCTCTGGAAGAT‑GCcgcg > 1:407850/1‑149 (MQ=255) agaagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑aaa < 1:31194/151‑1 (MQ=255) agaagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑aaa < 2:334715/151‑1 (MQ=255) gcTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCttt < 1:717413/151‑1 (MQ=255) gAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTg < 1:48370/151‑1 (MQ=255) gggtTCTACGATATCACCGCTTCCATCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑gg < 2:168053/149‑1 (MQ=255) tCTACGATATCACTGCTTACAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACGTTCCGTTTTGTGTGGCGGGTAGGACTTAGAATCTATGAGAT‑TCCGCG‑AAAAATTTTTCCTTTGGAGC‑ggg > 2:623934/1‑147 (MQ=255) tCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATa < 1:223503/150‑1 (MQ=255) tCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATa < 1:466538/150‑1 (MQ=255) aCGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGGAGAT‑GCCGCG‑GAAATTTTTTCCTTTTGCGC‑GGGTaaaa > 1:207499/1‑150 (MQ=255) atatCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTc < 2:197312/151‑1 (MQ=255) tatCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGAAGAT‑GCCCCG‑GAAATTTTTTTCTTTGGTGC‑GGGTAAAAttttc > 2:420649/1‑149 (MQ=255) aCCGCTTCCAGCGATGGCCGTGCGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGTGGGTTCAGTCCCTCGAAGTTGGCCACC‑AACTTTTTTTCTTTTTG‑GC‑GGGTACAATTTCCATc > 2:130088/1‑150 (MQ=255) aGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGGAGAT‑GCCGCG‑AAAATTCTTTCCTTTTGCGC‑GGATAAAATTTTCATCAACTCTcc > 2:76644/1‑149 (MQ=255) aGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGGAGAT‑GCCGCG‑GAAATTTTTTTCTTTGGCGC‑GGGTAAAATTTCCCTCAACTCTcc > 1:313977/1‑149 (MQ=255) cGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGGGTCTGGGAGAT‑GCCGCG‑AAAATTTTTTTCTTTGGGGC‑GGGTAAAATTTTCATCAACTCTCCTGCGCTgg > 2:98932/1‑148 (MQ=255) tAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGGTTAGTCTATGGGAGAT‑GCCGCG‑GAAATTCTTTTCTTTGGCGCGGGACAAAATTT‑CATCAACACTCCCTCGCTGGCAC‑ACCCa > 2:514691/1‑149 (MQ=255) tAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTTAGGCTCTGGGAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGGTAAAATTTTCATCAACTCTCCTGCGCTGGCAG‑ACCCaa > 2:212412/1‑150 (MQ=255) gATAAAAGCTGGGTATCTCGCGTGTCGGAAGTGATCGACATTCCGTTTTGTGGGGCGGGGGGGGCTAGGCCTCTGGGAGAT‑GCCCCG‑AAAATTTTTTGTTTTGGGGC‑GGGCAAAATTTTCCTCAAACCTCCTTGGCTTGCAG‑Accccaga > 2:68061/1‑146 (MQ=255) aaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACattaatt < 1:748791/151‑1 (MQ=255) gggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTc < 2:747554/150‑1 (MQ=255) gggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCg < 1:461817/151‑1 (MQ=255) tctcGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGGAGAT‑GCCGCG‑AAAATTCTTTTCTTTGGCGC‑GGATAAAATTTTCCTCAACTCTCCTGCGCTGGCCG‑ACCCAACATTAATTACTCGCCCgg > 2:203518/1‑151 (MQ=255) cgcgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCGCCTGGcc < 1:723992/150‑1 (MQ=255) gcgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCGCCTGGCCGa < 1:231543/151‑1 (MQ=255) gTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTAAGACTCTGGGAGAA‑GCCGCGAAAAATTTTTTCTTTTGGGCG‑GGAAAAAATTTTC‑TCCACTCCCCTGCCCTGGGGGAACCCCAC‑TTAATTATTCTCCTCGcccagc > 2:425809/1‑147 (MQ=255) gATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGt < 1:612460/150‑1 (MQ=255) aTCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGGCCCTTGAAGAT‑GGCCCG‑ACAATTTTTTCCTTTGGCGG‑GGATAAAAATTCCCTCAACTCTCCTGCGCTGGCAG‑ACCCCACATTAATTTTCCCCCTGGCCCATCGCTTTg > 2:502857/1‑145 (MQ=255) tCGACATTCCGTTTTGTGTGGCGGGGGGGGTTTAGTATCTGGGAGAT‑GCCGCG‑GAAATTTTTTTCTTTTGCGC‑GGGTAAAATTTCCCTCAACTCTCCTGCGCTGGCAG‑AACCCACATTTATTACTCCCCTGGCCGATCGCTTTTGCGTg > 2:70789/1‑149 (MQ=255) cGACATTCCGTTTTGTTTTGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCGCCTTGCCGATCGCTTTGGCGTGCa < 2:627978/150‑1 (MQ=255) cATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTTTCTGGGAGAT‑GCCGCG‑AAAATTTTTTCCTTTTGGGC‑GGGTAAAATTTTCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTTCTCGCCTGGGCCATCCCTTTGGCGTGCAgtg > 2:661856/1‑150 (MQ=255) cGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATtgt < 1:455448/151‑1 (MQ=255) gtgttGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATtgtg < 2:115901/148‑1 (MQ=255) gTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATAc > 1:341807/1‑151 (MQ=255) aTTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGt > 1:754902/1‑151 (MQ=255) | CAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGAT‑GCCGCG‑AAAATTCTTTCCTTTGGCGC‑GGATAAAATTTCCATCAACTCTCCTGCGCTGGCAG‑ACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGT > NC_000913/2095950‑2096240 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |