Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,921,173 | 0 | T | C | 45.0% | 49.8 / 10.3 | 40 | I310V (ATT→GTT) | gudP | putative D‑glucarate transporter |
Reads supporting (aligned to +/- strand): ref base T (19/3); new base C (0/18); total (19/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.19e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGAC > NC_000913/2921032‑2921322 | ccATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATcccaccc > 1:567838/1‑149 (MQ=255) aTAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGca < 1:736220/151‑1 (MQ=255) aTAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGc > 2:560044/1‑150 (MQ=255) ccGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACAccgc < 1:37979/151‑1 (MQ=255) tgatCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGAt > 2:84127/1‑151 (MQ=255) cGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaaccgcaaaccgc > 2:331153/1‑151 (MQ=255) tGACGAAGTTGCAGAACCCCACCCCCATCGACAGAAACATCCCCATTCCGACCGTTTTTTCCCGCCCAATGTCCGGCGACCCCGGGCGGCGAATCAGCCAGCCGGAAATAATCCCCCCCATCACACCGCCGATAAAACCGCAAAccgcc < 1:63332/149‑1 (MQ=255) aCGTAGTTGCAAAACCCCACCACCTCCGACAGCACCATGCCCATTACGACCGGTGTTTTACGCGAAATGTCCAGCGACCCCGGGCGGCGCAACAGCCAGTCGGAAATAACCCCACCCAGAACACCGCCGATAAAACCGCAAACCGCcggaa < 2:636586/151‑1 (MQ=255) cGTAGTGGCAAACCCCCACCCCCACCGACAGCAACTTCCCCTTTACAACCGGTTTTTTCCGCGCAATGTTCAGCGACCCCGTGCGGCGCATCAGCCAGTCGGAAATAACCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCcgga < 2:256725/149‑1 (MQ=255) aGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacgg > 1:333924/1‑150 (MQ=255) tccaaaaCCCCACCCCCATCGACAGCAACATCCCCATTACAACCGGGTTTTTCCGCCCATTGTTCAGCGACCCCGTGCGGCGCACCACCCGGCGGGAAAAAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAg < 2:703656/145‑1 (MQ=255) gCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGcc > 1:538695/1‑151 (MQ=255) cAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGAGGCGCATCAGTCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGAAAACCGCCGGAACGGAAGcc > 1:571563/1‑150 (MQ=255) aaCACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCcaca > 1:516931/1‑150 (MQ=255) acacCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACaa > 1:104697/1‑150 (MQ=255) catcaccaccGACAGCAACATGCCCATTCCGGTGGGTGTTTTCCGCGCAATGTTCAGCGACCCGGTGCGGCGCATCAGCCAGTCGGAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCGCAAAGccc < 2:492050/151‑1 (MQ=255) caccGACAGAAAAATGCCCATTACGAACGGTTTTTTCCGCGAATTGTCCAGCGCCCCGGTGCGGCGCACCGGCCAGTGGAAAAAAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc < 2:758139/148‑1 (MQ=255) cAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATc > 1:231784/1‑151 (MQ=255) aGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCg > 1:381104/1‑151 (MQ=255) aaCATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacg > 1:519019/1‑120 (MQ=255) ccccatacctaccGGTTTTTTACGCGAAAGGTACACCGCTCCCGCGGGGCCCACCACCCGGCCGGAAAAACCCCCCCCCAGCACCCCGCCGATAAAACCGCAACCCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcg < 1:669124/139‑1 (MQ=255) cccATTACGATCGGTTTTTTACGCGAAATGTTCAGCGCCCCCGTGCGGCGCATCAGCCAGTCGAAAATAACCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATccc < 2:122463/149‑1 (MQ=255) ccaaTTACGATCGTTGTTTTCCGCGAATTGTTCAGCGTCCCGGTGCGGCCCATCACCCGTTCGAAAATATCCCCCCCCACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATccc < 2:567799/146‑1 (MQ=255) ttttCCGCGCATTGTTCACCACCCCCGTGCGGCGCACCAGCCATTCGAAAATAACCCCACCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCa < 2:153442/144‑1 (MQ=255) tACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATaaa > 1:65166/1‑150 (MQ=255) cgcgCAATGCTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGACCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCg > 1:486376/1‑151 (MQ=255) acattGTCCAGCGCCCCCGTGCGGCCCACCAGCCAGCGGAAAATACCCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGa < 1:94183/147‑1 (MQ=255) cAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGa > 1:788458/1‑150 (MQ=255) aaTGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAAc > 1:694397/1‑151 (MQ=255) gTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAAcc > 1:3667/1‑149 (MQ=255) cAGCGACCCCGTGCGGCACACCACCCAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAg < 2:132048/148‑1 (MQ=255) cccGGGCGGCGAACCAGCCAGTCGGAAATAACCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAaag < 2:220325/150‑1 (MQ=255) ccGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAaaga > 2:167696/1‑150 (MQ=255) ccGTGCGGCGCATCAGCCAGCCGAAAATAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAaagaa < 2:455482/151‑1 (MQ=255) gcggccCACCACCCAGTCGAAAAAAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagtaa < 1:778784/151‑1 (MQ=255) ccacccAGTGGAAAAAAACCCCCCCCACCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAggg < 2:281647/146‑1 (MQ=255) cAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCg > 2:766909/1‑151 (MQ=255) caagccGGAAAAAACCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCg < 1:513484/143‑1 (MQ=255) tCGAAAAAAACCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGt < 2:244812/145‑1 (MQ=255) aaaaaaaCCCCCCCCGCCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGt < 1:617268/141‑1 (MQ=255) ataatCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACt > 1:232713/1‑151 (MQ=255) aatCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTg > 1:109395/1‑150 (MQ=255) atCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGAc > 1:663827/1‑151 (MQ=255) | CCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGAC > NC_000913/2921032‑2921322 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |