Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | KT2440 | 6,063,995 | 0 | A | C | 35.4% | 76.7 / 13.4 | 48 | V162G (GTG→GGG) | cheB_5 | Protein‑glutamatemethylesterase/protein‑glutamine glutaminase |
Reads supporting (aligned to +/- strand): ref base A (4/27); new base C (17/0); total (21/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.41e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AACGCCTTGATCACACATGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCAGGGCTTTCTTT > KT2440/6063863‑6064123 | accGCCTTGATCACCCATGGCTCCATGCCACTGGCCCGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCtt < 4:9235/135‑1 (MQ=255) aaCGCCTTGATCACACATGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCtt < 8:267277/137‑1 (MQ=255) aaCGCCTTGATCACACATGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCtt < 3:24990/137‑1 (MQ=255) gCCTTGATCACACATGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCGTGAc > 5:159258/1‑137 (MQ=255) gCCTTGATCACACATGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGGGGGCGCCGCCCCCATgcc > 5:281313/1‑135 (MQ=255) ccTTGATCACACATGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGGGGCGGCGCCGCCCCCATGAcc > 1:35176/1‑137 (MQ=255) aTCACACATGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTc < 8:228506/137‑1 (MQ=255) aTCACACATGGCTGCAGGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGTAGAATCCACCTTGACCGCTTc < 5:88486/137‑1 (MQ=255) caTGGCTGGATGCCACTGTCCAGGCGCAGCTGGCCGTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCTGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCtt < 8:252737/137‑1 (MQ=255) aTGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCTTGACCGCTTCGGGCTTg > 6:376766/1‑137 (MQ=255) tGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCACCGCCACCTTGACCGCTTCGGGCTTGc > 4:235284/1‑137 (MQ=255) gCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCATGCCCCCTTCGGGCGTGccc > 7:82425/1‑137 (MQ=255) aTGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTGGGCGCCGCGCCCCTGACCGCTTCGCGCTTgcgcgccgg > 3:42963/1‑134 (MQ=255) gCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGGGGGCGCCGCCCCCGAGACCGCTTCGGGCGTGGCCGcccgcc > 1:372980/1‑132 (MQ=255) tgcccAGGCGCAGCTGGCCCTGGCCACGTACTCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCCAAGGCTTGGCCGGTGCCGGTGTCGCCGCCGCCCCCCTGACCCCGACCGGCGTGACCGCAGTCAgcgcc > 6:325408/4‑136 (MQ=255) gCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCCCAGTCAGcgcggt > 6:172905/1‑135 (MQ=255) aGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCaa < 4:247140/137‑1 (MQ=255) gcAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAg < 4:354595/137‑1 (MQ=255) cAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGGCGGCGCCGCCCCCGTGCCCGCGTCGGGCGTGCCCGCCGGCCGCGCGGTGAACGAgg > 8:29833/1‑136 (MQ=255) cAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGGGTCGGCGCCGCCCCCGTGCCCGCGTCGGGCGTGCCCGCAGTCAGCGCGTGCAACGAgg > 4:254413/1‑136 (MQ=255) ggccctggccACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCTAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTg < 5:325408/137‑1 (MQ=255) ggccctggccACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTg < 3:254413/137‑1 (MQ=255) gcccgggccACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTgg < 7:136676/137‑1 (MQ=255) cACGGCCGCTCGGCGCGTGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGcc < 7:128930/137‑1 (MQ=255) cACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCGTGACCGCGTCGGGCTTGGCCGCAGGGCGCGGGGTCCACGAGGGGGTGGGAAAGGcc > 2:92676/1‑137 (MQ=255) aCGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCCCGTCGGGCTTGGCCGCAGTCAGCGCGTTGACAGAGTGGTTGGGAAAGGCCg > 8:328767/1‑137 (MQ=255) aCGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGGCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCg < 8:343939/137‑1 (MQ=255) cGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTGGGCGCCGCCCCCATGACCGCGTCGGGCTTGGCCGCCGGGCGCGCGGTGAACGCGGCGTTGGGGAAGGCCGGGGTg > 4:76330/1‑137 (MQ=255) cGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTg < 6:31017/137‑1 (MQ=255) cGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGGGTGGGCGGCGCCCCCCTGACCGGTTCGGGCGTGCCCCCCGTCAGCGCGTTGGACGGGGCGTTGGGAAAGGCCggggtt > 4:380354/1‑134 (MQ=255) ggCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGTGCCGCCACCTTGACCGCATCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCtac > 5:367236/1‑135 (MQ=255) ggCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCGTGACCGCTTCGGGCGTGCCCGCAGGCCGCGCGGTGCACGAGTCGGTGGCAAAGGCCGGGGTGGCTcg > 5:90784/1‑137 (MQ=255) ggCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGGGTCGGCGCCGCCCCCGTGACCGCGTCGGGCTTGCCCGCAGTCAGCGCGGTCAACGAGGCGTTGGCGAAGGCCGGGTTGGGTcg > 7:130230/1‑137 (MQ=255) gcggcAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGcc < 6:326650/137‑1 (MQ=255) gcggcAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGcc < 8:297045/137‑1 (MQ=255) gcggcAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGcc < 4:366516/137‑1 (MQ=255) cGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATc < 2:207927/137‑1 (MQ=255) cgccgcGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGACCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGc < 7:331554/137‑1 (MQ=255) gcgACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTc < 1:222826/137‑1 (MQ=255) gACCAAAGGCTTGGCCGGTGCCGGGGTCGGCGCCGCCCCCTTGGCCGCGTCGGGCTTGCCCGCAGGGCGCGCGGTGAACGAGTCGTTGGCAAAGGGCGGGGTGGGGCGCGCCGGGGCGCTGGCCAGGAGGCGGtggg > 2:319979/1‑134 (MQ=255) aCCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCCCCCTGGCCGCGTCGGGCTTGCCCGCCGTCAGCGCGGTCAACGCGTCGTTGGGAAAGGGCGGGTTGGGGGGCGCCGGGGCGCTGGGCCTCAGGCGGGCGAg > 2:322148/1‑137 (MQ=255) cAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCt < 1:322148/137‑1 (MQ=255) gccggtgccggtGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCtt < 7:379218/137‑1 (MQ=255) gccggtgTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCGAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCAg < 6:281313/137‑1 (MQ=255) gccggtgTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCAg < 3:76330/137‑1 (MQ=255) gccggtgTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCAg < 7:293694/137‑1 (MQ=255) gccggtgTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCAg < 6:159258/137‑1 (MQ=255) gcgcCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCAGGGctttcttt < 2:195370/137‑1 (MQ=255) | AACGCCTTGATCACACATGGCTGCATGCCACTGGCCAGGCGCAGCTGGCCCTGGCCACGGCCGCTCGGCGCGGGGGCGGCAGGTGCGGGCGCCGCGACCAAAGGCTTGGCCGGTGCCGGTGTCGGCGCCGCCACCTTGACCGCTTCGGGCTTGCCCGCAGTCAGCGCGTTCAACGAGTCGTTGGCAAAGGCCGGGTTGGCTCGCGCCGGGCCGATGGCCATCAGGCTGTCGAGCTTACCGATCTTGGTCAGGGCTTTCTTT > KT2440/6063863‑6064123 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |