Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,846,239 | T→C | 15.5% | T217A (ACC→GCC) | mutM ← | formamidopyrimidine‑DNA glycosylase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,846,239 | 0 | T | C | 15.5% | 59.0 / 4.5 | 33 | T217A (ACC→GCC) | mutM | formamidopyrimidine‑DNA glycosylase |
Reads supporting (aligned to +/- strand): ref base T (16/11); new base C (3/2); total (20/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-01 |
GGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGATGTTGCCC > NC_002947/5846091‑5846384 | gggCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGggtg < 2:79724/150‑1 (MQ=255) tACACGCTCGCACTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGCTAAAGTCGGGCAGGGTGGTACCg > 1:277869/1‑150 (MQ=255) gCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAACAAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTg > 2:875368/1‑150 (MQ=255) tCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGcc > 2:377358/1‑149 (MQ=255) ccTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGAt < 2:979433/150‑1 (MQ=255) tGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGgcc < 2:754776/150‑1 (MQ=255) ccAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGgccgccg > 1:547402/1‑150 (MQ=255) cTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGCGGTACCGCCCTGCCCGATGGCCGCCGACAGCAccc > 2:149986/1‑150 (MQ=255) ttCGCGCAACTCAGTGCCGCATACCTTGCCCGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGGCGCGCAGGGAGGTACCGCCCTGCCCGCTGGCCGCCGGCCGCAgccc > 2:333029/1‑149 (MQ=255) ttCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGCGGTACCACCCTGc > 1:528753/1‑127 (MQ=255) ttCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGCGGTACCACCCTGc < 2:528753/127‑1 (MQ=255) cgcgCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGtt < 2:826964/150‑1 (MQ=255) ccTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCa < 1:80589/144‑1 (MQ=255) ccTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCa > 2:80589/1‑144 (MQ=255) ccTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCAGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCa < 1:853040/150‑1 (MQ=255) cACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAg > 2:601480/1‑150 (MQ=255) gggCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAgcgcgc < 2:304518/148‑1 (MQ=255) gCCCGCCACGCCCGTACACAAACAGTTCGTGCTGGAGGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGccc > 1:753383/1‑150 (MQ=255) ccGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCgg < 2:547402/150‑1 (MQ=255) ccACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGtt > 1:910225/1‑112 (MQ=255) ccACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGtt < 2:910225/112‑1 (MQ=255) cGTACACACACCGTCCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGAGCAAGCCGCTCAGGGTGGTCCCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATcccgc < 2:277869/150‑1 (MQ=255) acaAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCcgcccgcc > 1:768063/1‑150 (MQ=255) tCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCCGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGt < 1:875368/150‑1 (MQ=255) tgctgGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCTCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCTAGTTTCAGGTAGCGAGCCCGGGAAATCCCGGCCGCCTCTCGGCGTggg > 2:293510/1‑150 (MQ=255) ccccGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACcggc > 2:380508/1‑150 (MQ=255) ccGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTGGATACcggccg > 2:424234/1‑148 (MQ=255) cccGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCGCCAGCGCCAGTTTCAGGTAGCGCGCCCCGGAAAGCCCGCCCGCGTCTCGGCGTCGGCCGATACCGGCGGGAaacaa > 2:917333/1‑150 (MQ=255) ccGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAaacaac > 1:650318/1‑150 (MQ=255) cACCGATAAAGTCGCGCAGGGTGGGACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCgg > 2:943348/1‑150 (MQ=255) gcgcAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGATGtt < 2:650318/150‑1 (MQ=255) gcAGGGCGGTACCACCCTGCGCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACcggcgg < 1:164002/120‑1 (MQ=255) gcAGGGCGGTACCACCCTGCGCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACcggcgg > 2:164002/1‑120 (MQ=255) aGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGATGTTGccc > 2:650101/1‑150 (MQ=255) | GGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACACAAACAGTTCCTGCTGGAAGTACCCCGGCTGCCCGTCGCCACCGATAAAGTCGCGCAGGGTGGTACCGCCCTGCTCGATGGCCGCCGCCAGCACCCGTTTGATCTCGATCGCCAGTTTCAGGTAGCGCGCCCGGGAAATCCCGCCCGCCTCTCGGCGTGGGTCGATACCGGCGGCAAACAACGCCTCGGTGGCGTAGATGTTGCCC > NC_002947/5846091‑5846384 |
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: ‑ |