Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,921,173 | 0 | T | C | 42.1% | 50.1 / 11.2 | 38 | I310V (ATT→GTT) | gudP | putative D‑glucarate transporter |
Reads supporting (aligned to +/- strand): ref base T (17/5); new base C (0/16); total (17/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.04e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATC‑AGCCAG‑TCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGAC > NC_000913/2921024‑2921322 | ggCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATC‑AGCCAG‑TCGGAAATAATc < 2:708388/151‑1 (MQ=255) tAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATC‑AGCCAG‑TCGGAAATAATCCCACCCAGcac > 2:55163/1‑151 (MQ=255) tAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAAAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATC‑AGCCAG‑TCGGAAATAATCCCACCCAGcac > 2:247716/1‑151 (MQ=255) aaaGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATC‑AGCCAG‑TCGGAAATAATCCCACCCAGcac > 2:530706/1‑150 (MQ=255) aaGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATC‑AGCCAG‑TCGGAAATAATCCCACCCAGCacac < 2:824607/151‑1 (MQ=255) aaGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATC‑AGCCAG‑TCGGAAATAATCCCACCCAGCacac < 2:776361/151‑1 (MQ=255) aGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATC‑AGCCAG‑TCGGAAATAATCCCACCCAGCACAcc > 1:196750/1‑151 (MQ=255) atcCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATC‑AGCCAG‑TCGGAAATAATCCCACCCAGCACACCGCCGATAaaac > 2:130058/1‑151 (MQ=255) acccACCCAACGTTGCCGTGGTTGCAGACCCCCACCACCATCGCCAGCACCATGCCCATTACGATCGGTGTTTTCCGCGCAATGTTCAGCGATCCGGTGCGGCGCATA‑AGCCAG‑TCGGAAATAACCCCCCCCAGCACACCGCCGATaaaa < 2:461296/148‑1 (MQ=255) ctaaaCGTTGACGTGGTTGCAGACCACCACCACCATCGACGGCAACATCCCCATTACGATCGGTGTTTTCCGCGCAATGTCAAGCGATCCCGGGCGGCGAAAC‑CCCCAG‑TCGGAAATAACCCCCCCCACCACACCGCCGATAaaaccgc < 1:365269/146‑1 (MQ=255) aCGTGGACGTAGTTGCAGAACCCCATCCCCATCGCCAGCAACATGCCAATTACGATCGGTTTTTTCCGCGCATTTTTCAGCGATCCCGTGCGGCGCACC‑AGCCAG‑TCGAAAATAACCCCCCCCAGCACACCGCCGATAaaaccgcaaaccg < 2:720789/151‑1 (MQ=255) cGTTGACGTAGTGCCAGAACACCATCCCCACCGACACCAACATGCCCATTCCGCTCGTTGTTTTCCGCCCAATGTTCAGCGATCCCGTGCGGCGCACC‑ACCCAG‑TCGAAAATAACCCCCCCCAGCACACCGCCGATAaaaccgcaaaccg < 1:149144/150‑1 (MQ=255) ttGCAAACCCCCACCCCCATCCACACCAACAGCCCCATTACGATCGGTGTTTTCCGCCCAAGGTACACCGACCCCGGGGGGCCCACC‑AGCCAG‑TGGGAAAAAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacgga < 1:806729/149‑1 (MQ=255) acacCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATC‑AGCCAG‑TCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACaa > 2:392021/1‑150 (MQ=255) atcaccatcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATC‑AGCCAG‑TCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCg > 2:816521/1‑151 (MQ=255) acaccatcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATC‑AGCCAG‑TCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGc > 1:686571/2‑151 (MQ=255) atcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATC‑AGCCAG‑TCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc > 2:75812/1‑150 (MQ=255) aCAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATC‑AGCCAG‑TCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAAt > 1:169194/1‑151 (MQ=255) aaCATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATC‑AGCCAG‑TCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGAc > 1:810311/1‑150 (MQ=255) cccATTCCGTTGGTTGTTTTCCGCGAAATGTCCAGCGCCCCCGGGCGGCGAACC‑ACCCAG‑CGGAAAATACCCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCa < 1:192072/136‑1 (MQ=255) ttACGATCGGTGATTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATC‑AGCCAG‑TCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGc > 1:28432/1‑151 (MQ=255) acagtCGGTTTTTTCGGCGCAATTTCCAGCGACCCCGGGCGGCGCATA‑GCCCAG‑TCGAAAATACTCCCCCCCGCCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCt < 1:606780/147‑1 (MQ=255) gATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATC‑AGCCAG‑TCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGc > 1:32392/1‑151 (MQ=255) cGGTGTTTTGCGCGTAATGCTCAGCGCTGCCGCGCGTCGCATC‑TGTCAG‑CATGTGATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCAACAAAACCAGCTTTGAGAATCGACATACCGCGTGCCTGc > 1:323318/1‑148 (MQ=255) gtttttttACGCGAATTTTCCAGCGACCCCGTGCGCCGCACC‑A‑CCAGTTCGAAAAAAACCCCCCCCCCCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCAc < 1:247675/145‑1 (MQ=255) ttccgcgCAATTTTCAGCGACCCCGTGCGGCGCACC‑ACCCAG‑TCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATaa < 2:757498/147‑1 (MQ=255) ttccgcgAAATGTCCAGCGACCCGGGGGGGCCCTTC‑GCCCAG‑GGGAAAAAAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATaa < 1:568510/147‑1 (MQ=255) gcgaattGTTCGGCGACCCGGGGCGGGCAATA‑AGCCGG‑CCGAAAATAATCCCCCCCACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAAcc < 2:142905/143‑1 (MQ=255) aatttccAGCGATCCCGTGCGGCG‑ATCAAGCCGG‑TCGAAAATAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAAc < 2:1412/145‑1 (MQ=255) ttCACCGACCCGGGGCGGCGCATC‑CCCCAG‑CCGAAAAAACCCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAg < 1:661002/150‑1 (MQ=255) tccAGCGCCCCCGGGCGCCGCACC‑ACCCAG‑CCGAAAATACCCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAg < 1:399329/148‑1 (MQ=255) tCCCGTGCGGCGCATC‑AGCCAG‑TCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaataa > 1:638620/1‑149 (MQ=255) cGTGCGGCGCATC‑AGCCAG‑TCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCAcc > 2:136251/1‑121 (MQ=255) ggccCATC‑ACCCAG‑TCGAAAATAACCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagtaa < 2:678142/149‑1 (MQ=255) cccATC‑AGCCAG‑TCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCa < 2:226034/149‑1 (MQ=255) tC‑AGCCAG‑TCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGc > 2:532098/1‑151 (MQ=255) tCGAAAATACCCCCCCCCACCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGAt < 2:169218/149‑1 (MQ=255) aaaTAACCCCCCCCACCAACCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATa < 1:135913/146‑1 (MQ=255) aaaatCCCCCCCACCACGCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACt < 2:353607/149‑1 (MQ=255) aatCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTg > 2:222666/1‑150 (MQ=255) atCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGAc > 1:212195/1‑151 (MQ=255) | GGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATC‑AGCCAG‑TCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGAC > NC_000913/2921024‑2921322 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |