Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,589,941 | G→A | 8.9% | A193T (GCG→ACG) | mccB → | methylcrotonyl‑CoA carboxylase biotin‑containing subunit beta |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,589,941 | 0 | G | A | 8.9% | 125.3 / 3.3 | 45 | A193T (GCG→ACG) | mccB | methylcrotonyl‑CoA carboxylase biotin‑containing subunit beta |
Reads supporting (aligned to +/- strand): ref base G (22/19); new base A (2/2); total (24/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01 |
GCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGC > NC_002947/4589801‑4590081 | gcggcgCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGc > 1:945847/1‑150 (MQ=255) gcCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGcc > 2:603592/1‑150 (MQ=255) cAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGCCGGCACGTACGCGCCgcca < 1:1280505/148‑4 (MQ=255) cAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGCCGGCACGTACGCGCCgcca > 2:1280505/1‑145 (MQ=255) aaCTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATg > 2:352888/1‑150 (MQ=255) cTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGc < 2:767543/150‑1 (MQ=255) gCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCcga > 2:634681/1‑150 (MQ=255) ccGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCcgacga > 1:1118677/1‑150 (MQ=255) ccGCCAGGACGAAGTGTTCCCCGACCGCGAGCATCTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCcgacga > 1:658467/1‑150 (MQ=255) aGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGt > 1:1419205/1‑150 (MQ=255) ccccGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCca > 2:993388/1‑150 (MQ=255) ccGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGc < 2:1248394/150‑1 (MQ=255) gAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCGTGGTGCGCcagcag < 2:381754/144‑1 (MQ=255) gAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCGTGGTGCGCcagcag > 1:381754/1‑144 (MQ=255) aGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCa < 2:35096/150‑1 (MQ=255) ttCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCGTGGTGCGCCAGCAGGCGACCATCTTc < 2:484377/150‑1 (MQ=255) aTTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATc < 2:949502/138‑1 (MQ=255) aTTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATc > 1:949502/1‑138 (MQ=255) ttttCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGcc > 1:111820/1‑150 (MQ=255) cttcAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCccgcc > 1:46921/1‑150 (MQ=255) ttcAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCccgccg > 2:794948/1‑150 (MQ=255) aaCCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTg < 2:945847/150‑1 (MQ=255) aCCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCACGTACGCGCCgcca < 1:362911/86‑4 (MQ=255) aCCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCACGTACGCGCCgcca > 2:362911/1‑83 (MQ=255) ccAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTGGt > 2:644377/1‑150 (MQ=255) aGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGCCACGCACAGCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTGGTGa > 2:1058866/1‑150 (MQ=255) gCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTGGTGaaa < 2:111820/150‑1 (MQ=255) gCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTGGTGaaa < 1:644377/150‑1 (MQ=255) gCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTGGTGaaa > 1:1200810/1‑150 (MQ=255) gCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTGGTGaaa < 1:993388/150‑1 (MQ=255) cAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTGGTGAAAgc > 2:675798/1‑150 (MQ=255) aaCATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCTGCAGTAGCTCAAGCAGTGCTTACAGTGCATCAGCAGGCGACCATCTTCCTCGCCGTCCCGCCGCTCGTGAAAgcc > 1:878559/1‑150 (MQ=255) gcgcgcAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTGGTGAAAGCCGCCACCg < 1:634681/150‑1 (MQ=255) ggCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTGGTGAAAGCCGCCACCGGCGAAgtg < 1:675798/150‑1 (MQ=255) gCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTGGTGAAAGCCGCCACCGGCGAAgtgg < 1:650840/150‑1 (MQ=255) cATCCCGCAGATCGCCGTGGTGTTGGGCTCGTGCACCGCTGGCGGGGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTGGTGAAAGCCGCCACCGGCGAAgtggt < 2:735385/150‑1 (MQ=255) aTCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTGGTGAAAGCCGCCACCGGCGAAgtggtg < 2:1419205/150‑1 (MQ=255) cccGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATCGTGCGCCAGCAGGCGAACATCGGCGTCACCTTCCTGTCTCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAg < 1:603592/150‑1 (MQ=255) gCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAgcgc > 2:497270/1‑150 (MQ=255) gATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGa < 1:516903/150‑1 (MQ=255) tCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAg > 2:1256299/1‑150 (MQ=255) gCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGAc > 1:1369055/1‑150 (MQ=255) tCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGAt < 1:1256299/150‑1 (MQ=255) aCCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCAt > 2:532286/1‑150 (MQ=255) gCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGc < 2:878559/150‑1 (MQ=255) | GCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTGCACCGCTGGCGGCGCGTACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGCCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGC > NC_002947/4589801‑4590081 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |