Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 819,244:1 | +C | 11.4% | coding (395/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,244 | 1 | . | C | 11.4% | 150.6 / 15.6 | 44 | coding (395/627 nt) | PP_0705 | DNA‑3‑methyladenine glycosylase |
Reads supporting (aligned to +/- strand): ref base . (18/21); new base C (3/2); total (21/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.58e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01 |
CAGTTGCTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGGAGACG > NC_002947/819099‑819389 | cAGTTGCTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATc < 2:761357/150‑1 (MQ=255) 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 < 1:352428/106‑6 (MQ=255) cTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTCGATC‑AACGCCtcgtca > 2:352428/1‑101 (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:629637/150‑1 (MQ=255) tCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTgg > 2:718480/1‑150 (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‑GAt > 2:809372/1‑50 (MQ=255) cTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GAt < 1:809372/50‑1 (MQ=255) cTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTggg > 2:442795/1‑119 (MQ=255) cTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTggg < 1:442795/119‑1 (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) | CAGTTGCTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTT‑GATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGGAGACG > NC_002947/819099‑819389 |
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: ‑ |