Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 650,464 | 0 | A | T | 44.4% | 48.6 / 15.6 | 36 | R9R (CGT→CGA) | citE | citrate lyase, citryl‑ACP lyase (beta) subunit |
Reads supporting (aligned to +/- strand): ref base A (3/17); new base T (15/1); total (18/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.06e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGC > NC_000913/650330‑650610 | gtgtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGT‑TTTTCGTTGTTTCAGCg > 1:58668/1‑148 (MQ=255) gtgtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGaa < 2:114727/150‑1 (MQ=255) gtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGc < 1:724642/150‑1 (MQ=255) gtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGT‑TTTTCGTTTTTGGAGCGAAGCg > 2:669170/1‑151 (MQ=255) gtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGTGTTTTCGTT‑CTGAAGCGAAg > 2:549635/1‑149 (MQ=255) ttttCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGaa < 2:715581/150‑1 (MQ=255) cACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGAGG‑TTTTTGGTGGTGCAGGGAAGCGGaaagtat > 1:734778/1‑147 (MQ=255) aTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCaa < 1:22813/151‑1 (MQ=255) aCGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGG‑TTTTCGTTGTTTCAGCGAAGCGGGAAAAAATTGCAAtcc > 1:542067/1‑149 (MQ=255) gAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGG‑TTTTCGTTGTTGCAGCGAAGCGGGAAATATTGGCAATCCTccc > 1:437685/1‑150 (MQ=255) aGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGCGG‑TTTTATTTGTTTCAGCGCAGAGGGAATCAATGGCaagacttc > 1:209952/1‑142 (MQ=255) gTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATg < 2:295018/151‑1 (MQ=255) ttCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGAGT‑TTTTCGCTTTTGTAGAGCAACGGGAAACATTTGCCATCCTccc > 2:744281/1‑145 (MQ=255) aaCATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCgg < 1:56090/148‑1 (MQ=255) aCATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGGGGGT‑TTTTCGTTGTTGCAGCGGAACGGGAATCAATGGGAAAcccccccagggcgaggcggg > 1:690397/1‑134 (MQ=255) aGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGGGGGT‑TTTTTGTTTTTGCAGCGCAAAGGAAAACATTTGCAAACCCCCCCTCGCCAGGCCGGGGTAccc > 2:215527/1‑148 (MQ=255) gggTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcag < 2:188182/151‑1 (MQ=255) gATGAAGGTGTTGCTGACCAGCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTTCGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCAGGATACCGCTGGCGCGTGCcagcatg < 2:505622/150‑3 (MQ=255) aTGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGAGGCGGGG‑TTTTAGTTGTTTCTGCGGAAGGGAAAACAATGGCAAACCCCCCCAGGCAGAGCCGGGAAACCCCCCGCGCGTGCcaccccgg > 1:89748/1‑144 (MQ=255) aaGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGT‑TTTTCGTTGTTTCAGCGAAGCGGAAATCATTTGCAATACTCCCCCGGCAGAGCCCGGATACCGCTGGCGCGTGCCCGCAGGGcct > 1:345396/1‑148 (MQ=255) ggAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCa < 2:168745/151‑1 (MQ=255) gACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCAc < 1:46485/151‑1 (MQ=255) tCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGG‑TTTTTGTTTTTGCCGCGAAGCGGGAAACAATTGCAAACCCCCCCTGGCAGAGCCGGGATACCGCCGGCGCGGGCCAGCAGGGGCTCCAGTCGTGCCCGTaa > 2:497897/1‑150 (MQ=255) cgcgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGG‑TTTTCGTTGTTGCAGCGGAGCGGGAATCATTGGCAATCCTCCCCTGGCCGAGCCGGGGAACCGCTGGCGCGGGGCCGCAGGGCTTCCAGGCGTGCACGTaaa > 1:524569/1‑150 (MQ=255) gcgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAAt < 2:78449/151‑1 (MQ=255) cATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGG‑TTTTTGGTGTTGGCACGGCACGGAAAACATTGGGAATCCCTCCCCGGCCGAGGCGCGATTACCCTGGTGGGGGCCCGCAGGGCGTCCAGTCGTGGACgggaa > 2:156288/1‑141 (MQ=255) ccAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTgc < 1:105722/151‑1 (MQ=255) acaaacaaCATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCttt < 1:609556/151‑1 (MQ=255) aacaTGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATc < 1:454492/151‑1 (MQ=255) caTGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGa < 1:435151/151‑1 (MQ=255) ggcggGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGc < 1:324659/151‑1 (MQ=255) cggGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTg > 2:212035/1‑151 (MQ=255) gTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAc > 1:75259/1‑151 (MQ=255) gCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAcg < 1:111324/150‑1 (MQ=255) gT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCcgcg < 1:431351/151‑1 (MQ=255) t‑tttACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCgcgc > 2:285400/1‑151 (MQ=255) | GTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGT‑TTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGC > NC_000913/650330‑650610 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |