Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,251,032 | C→A | 20.0% | G26V (GGT→GTT) | cidA ← | murein hydrolases holin regulator |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,251,032 | 0 | C | A | 20.0% | 37.4 / 4.9 | 20 | G26V (GGT→GTT) | cidA | murein hydrolases holin regulator |
Reads supporting (aligned to +/- strand): ref base C (6/10); new base A (2/2); total (8/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01 |
AGCCATCAGCCAGCCGGCCCCCAATTGCAGCATGGCCGGTTTGAGCACACCCGAAGCGAACAGTATCAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCGGGGTTCCTCGACAGGCCCTCATTTTAAGATTTGGCTGCCTAGGCCAGAAGCGAATTGTTAGACTGCTC > NC_002947/5250899‑5251180 | agccatcagccaGCCGGCCCCCAATTGCAGCATGGCCGGTTTGAGCACACCCGAAGCGAACAGTATCAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAgcgc > 2:53796/1‑150 (MQ=255) catcagccaGCCGGCCCCCAATTGCAGCATGGCCGGTTTGAGCACACCCGAAGCGAACAGTATCAACAGCAGTGCCAGGCCCATCACACAACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGcaca > 1:83008/1‑150 (MQ=255) atcagccaGCCGGCCCCCAATTGCTGCATGGCCGGTTTGAGCACACCCGAAGCGAACAGTATCAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACag > 1:176200/1‑150 (MQ=255) gccGGCCCCCAATTGCAGCATGGCCGGTTTGAGCACACCCGAAGCGAACAGTATGAACAGCAGTGCCACGCTCATCACTGCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGc < 2:147928/150‑1 (MQ=255) tGCAGCATGGCCGGTTTGAGCACACCCGAAGCGAACAGTATCAACAGCAGTGCCAGGCCCATCACACAACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACaaac < 2:83008/150‑3 (MQ=255) aGCACACCCGAAGCGAACAGTATCAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCAACAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTc > 1:295509/1‑150 (MQ=255) aGCACACCCGAAGCGAACAGTATCAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCAACAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTc > 1:295473/1‑150 (MQ=255) gAAGCGAACAGTATCAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCa < 1:129213/150‑1 (MQ=255) gCGAACAGTATCAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGt > 2:57388/1‑150 (MQ=255) tCAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCGATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCgg < 2:176200/150‑1 (MQ=255) cAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCAACAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCggg < 2:295473/150‑1 (MQ=255) cAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCAACAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCggg < 2:295509/150‑1 (MQ=255) aaCAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACGTGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCgggg < 1:442247/150‑1 (MQ=255) cAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCGGGGTTCCTCGACAg > 1:583246/1‑150 (MQ=255) cAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCGGGATTCCTCGACAg < 2:104286/150‑1 (MQ=255) cacaCCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCGGGGTTCCTCGACAGGCCCTCAtt < 2:557840/150‑1 (MQ=255) gccAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCGGGGTTCCTCGACAGGCCCTCATTTTAAGATTTGGCTGCCTAGGCCAGaa < 1:57388/150‑1 (MQ=255) cAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCGGGGTTCCTCGACAGGCCCTCATTTTAAGATTTGGCTGCCTAGGc < 2:211673/143‑1 (MQ=255) cAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCGGGGTTCCTCGACAGGCCCTCATTTTAAGATTTGGCTGCCTAGGc > 1:211673/1‑143 (MQ=255) tGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATACTGCAGGTTTCATGGCGGGGTTCCTCGACAGGCCCTCATTTTAAGATTTGGCTGCCTAGGCCAGAAGCGAAt < 1:53796/150‑1 (MQ=255) ccaacaAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTCTTTTCAATCATGCAGGTTTCATGGCGGGGTTCCTCGACAGGCCCTCATTTTAAGATTTGGCTGCCTAGGCCAGAAGCGAATTGTTAGACTGc < 2:201416/150‑1 (MQ=255) accaAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCGGGGTTCCTCGACAGGCCCTCATTTTAAGATTTGGCTGCCTAGGCCAGAAGCGAATTGTTAGACTGCTc < 2:56142/150‑1 (MQ=255) accaAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGAtttttt < 2:27270/60‑1 (MQ=255) accaAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGAtttttt > 1:27270/1‑60 (MQ=255) | AGCCATCAGCCAGCCGGCCCCCAATTGCAGCATGGCCGGTTTGAGCACACCCGAAGCGAACAGTATCAACAGCAGTGCCAGGCCCATCACACCACCGGGGATCGGCCAGCCCAGCCAACTGGCAATCTGGCCACCAAGCAGGAACAGCGCACAGAGAACCGCCAGCTCGACAAGCAGGCGCAGGGCTTTTTTCAATAATGCAGGTTTCATGGCGGGGTTCCTCGACAGGCCCTCATTTTAAGATTTGGCTGCCTAGGCCAGAAGCGAATTGTTAGACTGCTC > NC_002947/5250899‑5251180 |
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: ‑ |