| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 3,288,116 | G→C | 14.3% | R14G (CGC→GGC) | xthA ← | exodeoxyribonuclease III/apurinic/apyrimidinic endodeoxyribonuclease VI |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 3,288,116 | 0 | G | C | 14.3% | 74.2 / 6.1 | 28 | R14G (CGC→GGC) | xthA | exodeoxyribonuclease III/apurinic/apyrimidinic endodeoxyribonuclease VI |
| Reads supporting (aligned to +/- strand): ref base G (11/13); new base C (2/2); total (13/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01 | |||||||||||
GGCCTTTCTGGCCATGGTAGTGCACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTGGGGGTGAACGAAG > NC_002947/3287972‑3288263 | ggCCTTTCTGGCCATGGTAGTGCACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGccc > 2:24889/1‑150 (MQ=255) tCTGGCCATGGTAGTGCACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGc > 2:569396/1‑150 (MQ=255) cATGGTAGTGCACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCTGGCGGGCCCGCAGGCCGTTGa < 1:361058/150‑1 (MQ=255) gTAGTGCACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGtt < 2:631968/150‑1 (MQ=255) cACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACa < 1:24889/150‑1 (MQ=255) aGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGAAATTGACCTCCCCTACGGCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATg < 1:828997/150‑1 (MQ=255) ggCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCg < 1:448655/150‑1 (MQ=255) cGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTaaaa < 2:667051/150‑1 (MQ=255) gtcgCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCtg > 1:318082/1‑149 (MQ=255) gtcgCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTa > 2:399312/1‑150 (MQ=255) ccTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACAt > 1:463073/1‑150 (MQ=255) cTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCCCGCCCGCatca < 1:54960/80‑4 (MQ=255) cTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCCCGCCCGCatca > 2:54960/1‑77 (MQ=255) gTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGAcc < 1:399312/150‑1 (MQ=255) cAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCAc < 2:114080/150‑1 (MQ=255) cGATCACGTCCGGCTGGTGCTTTTCGATCGGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACATGCATTACCAGCa > 1:380785/1‑150 (MQ=255) gATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTa > 1:312957/1‑150 (MQ=255) tCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAAc < 1:569396/150‑1 (MQ=255) ccgccAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTg > 1:207040/1‑150 (MQ=255) cAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTgg > 1:583824/1‑150 (MQ=255) cAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTgg < 2:312957/150‑1 (MQ=255) aGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTg > 1:505293/1‑148 (MQ=255) aGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTg < 2:505293/148‑1 (MQ=255) gCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGCTGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTgggg < 2:207040/150‑1 (MQ=255) ggcctGCGGGCCCGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTggggg > 1:755116/5‑150 (MQ=255) gATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTgg > 2:197781/1‑145 (MQ=255) gATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTgg < 1:197781/145‑1 (MQ=255) cctGCGGGCCCGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTGGGGGTg < 2:755116/148‑1 (MQ=255) gcgggcCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTGGGGGTGAACGAAg < 2:463073/150‑1 (MQ=255) | GGCCTTTCTGGCCATGGTAGTGCACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTGGGGGTGAACGAAG > NC_002947/3287972‑3288263 |
| 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: ‑ |