| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 819,249 | Δ1 bp | 11.6% | coding (400/627 nt) | PP_0705 → | DNA‑3‑methyladenine glycosylase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 819,249 | 0 | G | . | 11.6% | 114.5 / 12.4 | 43 | coding (400/627 nt) | PP_0705 | DNA‑3‑methyladenine glycosylase |
| Reads supporting (aligned to +/- strand): ref base G (17/20); new base . (3/2); total (21/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.56e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.35e-01 | |||||||||||
GTTGCTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGGAGACGGCGCCGA > NC_002947/819101‑819396 | gTTGCTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGa > 2:89512/1‑150 (MQ=255) cGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGc < 1:829763/150‑1 (MQ=255) aCCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTggg < 1:761861/150‑1 (MQ=255) cGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAggg < 1:742675/150‑1 (MQ=255) cTGCGGCTTTTCGGCGAGCAAGCTGGGGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTCCTGCCGACGGCGCAGGCGTTATCAATGGTTGACGAGGAGTT‑GCTCGATCGCGTGATGGCTCTTGTTGTGGTAGGGcggtg < 2:225920/150‑1 (MQ=255) gCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacgg < 1:294890/150‑1 (MQ=255) ttttCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacggtg > 1:803816/1‑150 (MQ=255) gCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGGTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGAtgctgc > 1:115268/1‑150 (MQ=255) aGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGa < 2:168011/150‑1 (MQ=255) gCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGAGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGGGGGGATGGTGCGGag > 1:168376/1‑149 (MQ=255) cTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTCGATC‑AACGCCtcgtca > 2:352428/1‑101 (MQ=255) cTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTCGATC‑AACGCCtcgtca < 1:352428/106‑6 (MQ=255) gTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCtat < 2:702426/150‑1 (MQ=255) tGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacggtgga > 2:211062/1‑133 (MQ=255) cgACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCTTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacggtgga < 1:211062/132‑1 (MQ=255) cGGTCCAGGGCAGTGCTAAAGCCACCCGGGAGGGCCTGGGACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGcc < 2:279761/150‑1 (MQ=255) tCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTgg > 2:718480/1‑150 (MQ=255) tCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTgg < 2:629637/150‑1 (MQ=255) cAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTgggg > 2:89088/1‑150 (MQ=255) cTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTggg < 2:360047/139‑1 (MQ=255) cTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTggg > 1:360047/1‑139 (MQ=255) cTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACa > 2:715786/1‑150 (MQ=255) tCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACAt > 1:672420/1‑150 (MQ=255) cAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACAtt < 1:89088/150‑1 (MQ=255) aCCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCg > 1:205306/1‑150 (MQ=255) cccTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCgg < 2:672420/150‑1 (MQ=255) gAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGAc < 1:718480/150‑1 (MQ=255) cTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTggg < 1:442795/119‑1 (MQ=255) cTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTggg > 2:442795/1‑119 (MQ=255) aCCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGt < 1:118707/150‑1 (MQ=255) gACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGTACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTCGACGGTGGGGATGCTGCTGATCTATCGCCTGGGGCGTTCTGACATTCGGTCGGTGGACGATGTTGGCGTgtg < 2:779988/150‑3 (MQ=255) gACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTgcg > 2:326469/1‑150 (MQ=255) aCCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTc > 1:608908/1‑150 (MQ=255) cATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCaa > 2:655157/1‑150 (MQ=255) tGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAgg < 1:241278/150‑1 (MQ=255) cTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGc > 1:339915/1‑150 (MQ=255) caccaGGCGTTCGATC‑AACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTgg < 1:725124/102‑1 (MQ=255) caccaGGCGTTCGATC‑AACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTgg > 2:725124/5‑106 (MQ=255) ccaGGCGTTCGATC‑AACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGga > 1:597332/3‑150 (MQ=255) gCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGCATCGGCGGCTCAACGGGCTGGAGACg > 2:799082/1‑150 (MQ=255) gCGGT‑GATCGGACGCCTGGTGGCGCTGCGTGGGGTAGGGCAGTGGACGGTGGAGATGCTGCTGATCTAGAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATNTTGGCGTGCGCGAGGAGTATCGGCGGCTCAAGGGGCTGGAGACg < 1:799082/150‑1 (MQ=255) taTCAAACCCCTGGTGGCGCGGCGTGGGGTAGGTAGGAGAACGGAGGAGATGCTGCTGATGTATAGCCTGGGGCGTTCTGACATGTTGCCGGGGGACGGGTGTTGTGTGCGCGCGCGGTAATCACCCCTCAGAGAGCTGGAGACGGTGcc > 2:97493/2‑150 (MQ=255) tCGCACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGGAGACGGCGCCGa < 1:89512/150‑1 (MQ=255) | GTTGCTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGGAGACGGCGCCGA > NC_002947/819101‑819396 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |