Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,846,019 | G→A | 13.5% | intergenic (+35/+56) | yfhL → / ← mutM | 4Fe‑4S cluster‑containing protein/formamidopyrimidine‑DNA glycosylase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,846,019 | 0 | G | A | 13.5% | 84.2 / 6.0 | 37 | intergenic (+35/+56) | yfhL/mutM | 4Fe‑4S cluster‑containing protein/formamidopyrimidine‑DNA glycosylase |
Reads supporting (aligned to +/- strand): ref base G (17/15); new base A (3/2); total (20/17) | |||||||||||
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.10e-01 |
GCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCC > NC_002947/5845876‑5846163 | gCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCCCGGTCCCTGTGAGAGCGGGCATg > 1:494236/1‑150 (MQ=255) tGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGGGAGCAGGCATGCCCGCTCTCACAggga < 1:732771/150‑4 (MQ=255) tGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGttt < 1:956255/150‑1 (MQ=255) tGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGttt < 1:956249/150‑1 (MQ=255) gACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGc > 1:541651/1‑150 (MQ=255) cATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCCCGGCCCCGGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCaa < 1:29128/150‑1 (MQ=255) gCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGtt < 1:91594/150‑1 (MQ=255) aTGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAg > 1:845531/1‑150 (MQ=255) tGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGc < 2:541651/150‑1 (MQ=255) gcgcATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAAc > 1:425070/1‑150 (MQ=255) cATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGtt < 1:1017632/150‑1 (MQ=255) gacgaaCAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCa < 2:285473/147‑1 (MQ=255) aagaaCAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAg > 2:940905/1‑150 (MQ=255) tGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGCGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAg < 2:740188/112‑1 (MQ=255) tGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAg > 1:740188/1‑112 (MQ=255) gATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGAAGt > 1:138266/1‑150 (MQ=255) gCCAAGTATCGCCGGATCACTGGCAATGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCg < 2:587259/98‑1 (MQ=255) gCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCg > 1:587259/1‑98 (MQ=255) ccAAGTATCGCCGGATCACTGGCAAGGGCTGCTCATTGCCCTTCCCTGCACGGTCCCTGTGGGAGCAGGCATGCCCGCTCTcaca > 2:490614/1‑85 (MQ=37) ccAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGGGAGCAGGCATGCCCGCTCTcaca < 1:490614/85‑1 (MQ=38) ccAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACg < 2:248694/150‑1 (MQ=255) aaGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCt > 2:806510/1‑150 (MQ=255) aCTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAg < 2:406210/150‑1 (MQ=255) agagCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACt > 1:79724/4‑150 (MQ=255) cccTTCCCTGCGCGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGt > 2:775628/1‑150 (MQ=255) ccTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAgcgc > 2:287548/1‑94 (MQ=255) ccTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAgcgc < 1:287548/94‑1 (MQ=255) ccTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCggg < 2:439758/97‑1 (MQ=255) ccTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCggg > 1:439758/1‑97 (MQ=255) tGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATa > 1:754776/1‑150 (MQ=255) ccGGTCGCTGCGAGAGCGGGAATGCCTGGTCCCACGGATTAGGTTGCGCTGCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCt > 2:853040/2‑150 (MQ=255) cGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCtt < 2:845531/150‑1 (MQ=255) ctaaacTGTGGGAGCAGGCATGCCCGCTCTCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTg > 2:481009/6‑150 (MQ=255) gTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGc > 1:826964/1‑150 (MQ=255) ccTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACg > 1:275871/1‑150 (MQ=255) cTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACgg < 2:425070/150‑1 (MQ=255) ggagCAGGCATGCCCGCTCTCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGAGCGGGCAGGACACGCTCGCACTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGcc > 2:533566/2‑150 (MQ=255) | GCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATCACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCC > NC_002947/5845876‑5846163 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |