Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,096,090 | 0 | A | G | 46.3% | 40.4 / 11.0 | 41 | I83V (ATT→GTT) | hisF | imidazole glycerol phosphate synthase, catalytic subunit with HisH |
Reads supporting (aligned to +/- strand): ref base A (6/16); new base G (19/0); total (25/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.33e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGG > NC_000913/2095942‑2096239 | gCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGtt > 1:867297/1‑151 (MQ=255) gCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGtt > 1:888777/1‑151 (MQ=255) gCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGAGCGACATTCCGTTTTGTGTGGCGGGTGGGAtt < 2:1047027/151‑1 (MQ=255) cGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGtt > 1:310652/1‑149 (MQ=255) aaaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAA‑GTCTCTgg > 1:95133/1‑151 (MQ=255) aaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTGA‑GTatcttg > 2:1053710/1‑148 (MQ=255) acgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGc < 1:526183/151‑1 (MQ=255) tGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTTA‑GGTTCTGGAAGATGCcgc > 2:806671/1‑150 (MQ=255) agaagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGaaa < 1:901473/151‑1 (MQ=255) gcTGACGAACTGGTGTTCTACGATATCCCCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAA‑GTCTCTGGGAGCAGTCGCGAAAATTCttt > 2:790273/1‑151 (MQ=255) cGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCttt < 2:172591/151‑1 (MQ=255) tACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGTGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATa < 2:146816/148‑1 (MQ=255) cGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTGA‑GTCCCTGGGAGATGCCGCGAAAATTTTTTCCTTTGGCGCGGGTaaaaat > 2:120450/1‑149 (MQ=255) gATATCACCGCTTCGAGCGATGGACGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGACGTGATCGACATTCCGTTTTGTGGGGCGGGGGGGAGTAAGGTAGGTGG‑AGAAGACGCGGAATATCTTTTCCTTGGGGCGGAAAAAAt > 2:179144/1‑149 (MQ=255) atCACCACTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCCCGCGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAA‑TAAACCGGCAGATGCCCCGTAAAATCTTTTCTTTGTGGCGGagacaaatt > 2:152562/1‑140 (MQ=255) ccGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCaa < 2:231333/151‑1 (MQ=255) cTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAA‑GTCTCTGGGAGGAGCCGCGAAAAATTTTTCCTTTGGGGCGGATAAAATTTTCATCaa > 2:27837/1‑148 (MQ=255) ccAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTcc < 1:1000333/151‑1 (MQ=255) gATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTTCTTTGGCGCGGATAAAATTTTCATCAACTCTcc > 1:200356/1‑146 (MQ=255) cGTGTGGGAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGTGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTACCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGc < 2:931210/149‑1 (MQ=255) tgtgGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATTAA‑GGCTCTGGAAGATGCCGCGAAAATTTTTTTCTTTGGCGCGGGGAAAAATTTCATCAACCCTCCTGCGCTGGCaaa > 2:55671/1‑148 (MQ=255) gtgGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTTGGATTAA‑GTGTCTGGAAGATGCCCCGGAAATTTTTTCCTTTTGCGCGGGTAAAATTTCCCTCAACTCTCCTGCGGTGGccaacc > 2:379403/1‑146 (MQ=255) gtgGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAA‑GGCTCTGGGAGAAGCCCCGAAAATTCTTTTCTTTGGGGCGGATAAAATTTCCCTCAACTCTCCTGCGCTGGGAGAcc > 2:135556/1‑151 (MQ=255) tgGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTTA‑GTCTCTGGGAGAAGCCGCGAAAATTCTTTTCTTTTGCGCGGATAAAATTTCCCTCAACCCTCCTGCGCTGGCCGAccc > 2:692838/1‑151 (MQ=255) aTAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTTA‑GGCCCTGGAAAAAGCCGCGGAAATTGTTTCCTTTTGGGCGGATAAAATTTCCATCAACTCTCCCGCGCCTGCAGACCCaa > 2:120164/1‑147 (MQ=255) aTAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACatt < 1:696540/151‑1 (MQ=255) aaaaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACatta < 1:338209/150‑1 (MQ=255) gggTATCTCGCTTGTCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTGCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACt < 2:798329/149‑1 (MQ=255) tctcGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTTGCGCGGATAAAATTTCCCTCAACTCTCCTGCGCTGGCAGACCCAACATTAATTAATCGCCTgg > 1:1054823/1‑151 (MQ=255) cgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGAt < 2:888175/151‑1 (MQ=255) cGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCtt < 1:630962/151‑1 (MQ=255) cGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCtt < 2:937948/151‑1 (MQ=255) ggAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGATAA‑GACTCTGGAAGATGCCCCCAAAATTTTTTCCTTTGGCGGGGAGAAAAATTTCCTCAACTCTCCTGCGGTGGGAGACCCCAAATTAATTACTTGCCTTGCCGATCGCttt > 2:683110/1‑151 (MQ=255) gAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGATAA‑GTCTCTGGAAGATGGCCCGAAAATTCTTTTCTTTGGCGCGGTGAAAATTTCCATCAAATCTCCCTCGCTGGCAGACCCAACATTTATTACTCGCCTTGGCGAACGCctt > 2:935195/1‑150 (MQ=255) gAAGTGATCGAAATTCCGTTTTGTGTGGCGGGCGGGGTTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTTGCGCGGGTAAAATTTCCATCAAATCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTCGCCGATCGc > 2:97883/1‑147 (MQ=255) cATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAgtgt < 2:928614/151‑1 (MQ=255) aTTCCGTTTTGTGTGGCGGGGGGGGTTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTTCATCAACTCTCCTGCGCCGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAg > 2:882687/1‑147 (MQ=255) tgtgtgGCGGGTGGGATAAA‑ATCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTTGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTAATCCCCTGGCCGCTCGCTTTTGCGTGCAGTGTATTGTGGTc > 2:84532/1‑150 (MQ=255) gCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTg > 2:263897/1‑151 (MQ=255) gggATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCt > 2:274806/1‑151 (MQ=255) gATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTgg < 1:183175/151‑1 (MQ=255) | GCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGG > NC_000913/2095942‑2096239 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |