Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_002947 | 5,241,585 | A→C | V96G (GTG→GGG) | hmgC ← | maleylacetoacetate isomerase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,241,585 | 0 | A | C | 96.4% | 92.6 / ‑3.3 | 28 | coding (287/633 nt) | hmgC | maleylacetoacetate isomerase |
Reads supporting (aligned to +/- strand): ref base A (0/0); major base C (15/12); minor base . (0/1); total (15/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.64e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01 |
TCCACTGCCGCGAGCCCTTGGCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGC > NC_002947/5241447‑5241728 | tCCACTGCCGCGAGCCCTTGGCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGaa > 4:156837/1‑148 (MQ=255) cgcgAGCCCTTGGCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGcc > 4:82993/1‑148 (MQ=255) gCCCTTGGCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGGGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGc < 7:137626/148‑1 (MQ=255) cccTTGGCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGcc > 4:170106/1‑148 (MQ=255) ccTTGGCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCt > 7:122235/1‑148 (MQ=255) ttGGCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCCGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTcg > 3:171438/1‑148 (MQ=255) gCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGGTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTcgg < 1:62314/146‑1 (MQ=255) cTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTcggcgg > 1:20024/1‑148 (MQ=255) ccACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCa > 8:11772/1‑148 (MQ=255) aCTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGACCCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAgg > 8:131086/1‑148 (MQ=255) cccTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAAc < 4:178508/146‑1 (MQ=255) cccTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACct > 4:152948/1‑148 (MQ=255) cTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGCCCCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCCCCGCAGGCTGTGGATAAACctct > 1:174354/1‑148 (MQ=255) tCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCtt > 1:16282/1‑148 (MQ=255) tCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTcc < 8:153269/147‑1 (MQ=255) tCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCACCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCATTCTGTGGATAAACCTCTTc < 6:16617/146‑1 (MQ=255) gCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCACCCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCga > 7:133113/1‑148 (MQ=255) aGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCGTCCAGGTACTCGATGATCGCCGGCGCCTGcc > 6:104299/1‑147 (MQ=255) gCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCAc < 5:104299/148‑1 (MQ=255) gCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCC‑CCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCAc < 6:47532/147‑1 (MQ=255) cacaCTGCCGTTGGGCAGCGGATTAATGGCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCAcc < 7:15935/148‑1 (MQ=255) acacTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCa > 3:76504/1‑148 (MQ=255) cacTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAg < 8:122235/148‑1 (MQ=255) acTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGc < 3:156837/148‑1 (MQ=255) cGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCAcc < 2:16282/146‑1 (MQ=255) aGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGt < 4:76504/148‑1 (MQ=255) cgccCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAGGGCCGCGTAC‑CGGCCTTGc < 3:82993/148‑1 (MQ=255) cgccCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCGGCCAGGTCCTCGATGAGCGCCGGCGCCTGCCCCAGCACCTCCCCACCGTCGGTGCGCAAGGCCG‑GTCCGCGGCCCTGc > 1:58113/1‑148 (MQ=255) | TCCACTGCCGCGAGCCCTTGGCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCG‑GTACGCGGCCCTGC > NC_002947/5241447‑5241728 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |