| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,185,043 | T→C | 13.3% | N489D (AAT→GAT) | aceK ← | isocitrate dehydrogenase kinase/phosphatase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 5,185,043 | 0 | T | C | 13.3% | 70.6 / 3.4 | 30 | N489D (AAT→GAT) | aceK | isocitrate dehydrogenase kinase/phosphatase |
| Reads supporting (aligned to +/- strand): ref base T (15/11); new base C (2/2); total (17/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01 | |||||||||||
ATATCGGCAAACAGAAACGGCGGAAACTCTTCGGGGAACACATCATGCGGCCCGATCGAGTACCACGGCTCGCCCGACATTTCGTCCTCCGGGTAACGCGGTGGCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCA‑CACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATG‑TTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGC > NC_002947/5184922‑5185188 | atatCGGCAAACAGAAACGGCGGAAACTCTTCGGGGAACACATCATGCGGCCCGATCGAGTACCACGGCTCGCCCGACATTTCGTCCTCCGGGTAACGCGGTGGCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCa < 2:329836/150‑1 (MQ=255) ggAAACTCTTCGGGGAACACATCATGCGGCCCGATCGAGTACCACGGCTCGCCCGACATTTCGTCCTCCGGGTAACGCGGTGGCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATATTTGTAGAACACGACCCgg > 2:359135/1‑150 (MQ=255) gAAACTCTTCGGGGAACACATCATGCGGCCCGATCGAGTACCACGGCTCGCCCGACATTTCGTCCTCCGGGTAACACGGGGGCGGAATGTGGCGGAAATTCACGTCGGAAAGGAAGTTGATGTCTGCATAGTCGTAGAACACGACCCGGc > 1:354621/1‑150 (MQ=255) aaaCTCTTCGGGGAACACATCATGCGGCCCGATCGAGTACCACGGCTCGCCCGACATTTCGTCCTCCGGGTAACGCGGTGGCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGcc > 1:416080/1‑150 (MQ=255) ttCGGGGAACACATCATGCGGCCCGATCGAGTACCACGGCTCGCCCGACATTTCGTCCTCCGGGTAACGCGGTGGCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCg < 2:662805/118‑1 (MQ=255) ttCGGGGAACACATCATGCGGCCCGATCGAGTACCACGGCTCGCCCGACATTTCGTCCTCCGGGTAACGCGGTGGCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCg > 1:662805/1‑118 (MQ=255) ttCGGGGAACACATCATGCGGCCCGATCGAGTACCACGGCTCGCCCGACATTTCGTCCTCCGGGTAACGCGGTGGCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACg > 2:444915/1‑150 (MQ=255) ggggAACACATCATGCGGCCCGATCGAGTACCACGGCTCGCCCGACATTTCGTCCTCCGGGTAACGCGGTGGCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCt > 1:531634/1‑108 (MQ=255) ggggAACACATCATGCGGCCCGATCGAGTACCACGGCTCGCCCGACATTTCGTCCTCCGGGTAACGCGGTGGCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCt < 2:531634/108‑1 (MQ=255) ggggAACACATCATGCGGCCCGATCGAGTACCACGGCTCGCCCGACATTTCGTCCTCCGGGTAACGCGGTGGCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGt < 1:537010/150‑1 (MQ=255) ggggAACACATCATGCGGCCCGATCGAGTACCACGGCTCGCCCGACATTTCGTCCTCCGGGTAACGCGGTGGCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGt > 1:148488/1‑150 (MQ=255) ggAACACATCATGCGGCCCGATCGAGTACCACGGCTCGCCCGACATTTCGTCCTCCGGGTAACGCGGTGGCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGa < 2:557371/143‑1 (MQ=255) ggAACACATCATGCGGCCCGATCGAGTACCACGGCTCGCCCGACATTTCGTCCTCCGGGTAACGCGGTGGCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGa > 1:557371/1‑143 (MQ=255) atcatGCGGCCCGATCGAGTACCACGGCTCGCCCGACATTTCGTCCTCCGGGTAACGCGGTGGCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCA‑CACCaaa > 2:436859/1‑150 (MQ=255) aGTACCACGGCTCGCCCGACATTTCGTCCTCCGGGTAACGCGGTGGCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCA‑CACCAAAGTTCTTCAGCAACATGt > 2:322659/1‑150 (MQ=255) ccACGGCTCGCCCGACATTTCGTCCTCCGGGTAACGCGGTGGCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCA‑CACCAAAGTTCTTCAGCAACATGTCGcc > 1:15794/1‑150 (MQ=255) aCGGCTCGCCCGACATTTCGTCCTCCGGGTAACGCGGTGGCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCA‑CACCAAAATTCTTCAGCAACACGTCGCCgg < 2:374042/150‑1 (MQ=255) cATTTCGTCCTCCGGGTAACGCGGTGGCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCA‑CACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATG‑TTcgc < 2:148488/150‑1 (MQ=255) cGGGTAACGCGGTGGCGGAATGTGGCGGAAATCCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCAGAACCCGAG‑AGTGCAGCAACATGTCGCCGGGGAAGATG‑TTCGCCGCCGCCAACTg > 2:584104/1‑150 (MQ=255) cGGGTAACGCGGTGGCGGAATGTGGCGGAAATCCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTAA‑CGAGGAGGTTCGTCAGCAACATGTCGCCGGGGAAGATG‑TTCGCCGCCGCCAACTg > 2:584213/1‑150 (MQ=255) gTAACGCGGTGGCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCA‑CACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATT‑TTCGCCGCCGCCAACTGCtt < 1:525796/150‑1 (MQ=255) gTAACGCGGTGGCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCA‑CACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATG‑TTCGCCGCCGCCAACTGCtt < 2:15794/150‑1 (MQ=255) tGGCGGAATGTGGCGGAAATCCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCA‑CACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATG‑CGGGCCGCCGCCAACTGCTTGATGGCCAg < 1:584213/150‑1 (MQ=255) tGGCGGAATGTGGCGGAAATCCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCA‑CACCAAAGTTCTTCAGCAACATGTCGCCGGGGAA‑ATGCGCCGTCGCCGCCAACTGCTTGATGGCCAg < 1:584104/150‑1 (MQ=255) gCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCA‑CACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATG‑TTCGCCGCCGCCAACTGCTTGATGGCCAGGc > 2:96818/1‑150 (MQ=255) ggAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCA‑CACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATG‑TTCGCCGCCGCCAACTGCTTGATg > 1:77607/1‑141 (MQ=255) ggAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCA‑CACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATG‑TTCGCCGCCGCCAACTGCTTGATg < 2:77607/141‑1 (MQ=255) aaTGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCA‑CACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATG‑TTc > 1:587709/1‑118 (MQ=255) aaTGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCA‑CACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATG‑TTc < 2:587709/118‑1 (MQ=255) aaaTTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCA‑CACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATG‑TTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAgc > 1:399901/1‑150 (MQ=255) | ATATCGGCAAACAGAAACGGCGGAAACTCTTCGGGGAACACATCATGCGGCCCGATCGAGTACCACGGCTCGCCCGACATTTCGTCCTCCGGGTAACGCGGTGGCGGAATGTGGCGGAAATTCACCTCGGTCAGGAAGCTGATTTCGTCATAGTCGTAGAACACGACCCGGCCATGACGGGTCA‑CACCAAAGTTCTTCAGCAACATGTCGCCGGGGAAGATG‑TTCGCCGCCGCCAACTGCTTGATGGCCAGGCCGTAGTCCTCCAGC > NC_002947/5184922‑5185188 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |