Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 4,390,131:1 | +T | 39.1% | coding (27/309 nt) | hfq → | global sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 4,390,130 | 1 | . | T | 39.1% | 7.3 / 38.3 | 23 | coding (26/309 nt) | hfq | global sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication |
Reads supporting (aligned to +/- strand): ref base . (9/5); new base T (7/2); total (16/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.57e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01 |
AGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGA‑TCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAA > CP009273/4389986‑4390271 | aGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGA‑TCCg > 1:317092/1‑149 (MQ=255) gttgttGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGA‑TCCGTTCCTGAACg > 1:316123/1‑149 (MQ=255) gttgGTGCTATCTCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGA‑TCCGTTCCTGAACGCAc < 1:346697/149‑1 (MQ=255) tgGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATTCCGTTCCTGAAc < 2:27819/144‑1 (MQ=255) gTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGA‑TCGGTTCCTGCACGCa > 2:262994/1‑144 (MQ=255) gTGCTATCGCAGGCTGAATGTGTACAATTGAGACGCATAGAGCGCAATATTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAAACTTGACAAGG‑TCCGTTCAGGAACGCAATGCg > 1:36468/1‑149 (MQ=255) tGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGA‑TCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTc < 1:234360/149‑1 (MQ=255) gAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATACGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGA‑TCCGTTCCTGAACGCACTGCGTCGGGAACCGGTTcc > 1:154976/1‑149 (MQ=255) gtgtACAATTGAGACGTATCGTGCGCAGTTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATACGGAAAAGAGAGAATGGCTAAGGGGCAATCTTCACAAGA‑TCCGTTCCTGAAAGCACTGCGTCGGGAACGGGTTc > 2:182816/1‑144 (MQ=255) gtgtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATATTTACAAGATTCCGTTCCTGAACGCACTGCGTCGGGAACGGGtt > 2:166548/1‑144 (MQ=255) gtACAATTGAGACATATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATTACGATCCTGAACGCACTGCGTCGGGCACGTGTTCCAGgtt > 1:281397/1‑149 (MQ=255) tACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTc > 1:63589/1‑149 (MQ=255) cgcaTCGTGCGCAATTTTTACAGAACCGAAAGGTTCAAAGTACAAACAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGA‑TCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGt < 1:236969/146‑1 (MQ=255) ttttttCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAAGGGCTAAGAGGCAATCTTTACAACATTCCGTTCCTGAACGCACAGCGTCGGGAAAGTGGTCCAGTTTCTATTTAGTTGGTGGATGGTa > 2:231360/1‑144 (MQ=255) tACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATTCCGTTCCTGCACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTGGGTGAATGGTATTAAGCTGCAAGGGCCAATCGAGTCTTTTGa > 1:330548/1‑149 (MQ=255) tACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTgt > 1:366612/1‑148 (MQ=255) taAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGCTTCTATTTATTTGGTGAATGGGATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGt > 2:38439/1‑144 (MQ=255) ggAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGA‑TCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGtt > 1:360486/1‑149 (MQ=255) agagaATGGCTAAGGGGCAATCTTTACAAGA‑TCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAac > 1:140094/1‑149 (MQ=255) aaTGGCTAAGGGGCAATCTTTACAAGA‑TCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGaaaaa > 2:313025/1‑144 (MQ=255) aaTCTTTACAAGA‑TCCGTTCCTGACCGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGtt < 1:154531/149‑1 (MQ=255) tACAAGATTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGTTCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACaa < 1:260325/149‑1 (MQ=255) cAAGA‑TCCGTTCCGGAACGCACGGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTAc < 2:203464/144‑1 (MQ=255) | AGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGA‑TCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAA > CP009273/4389986‑4390271 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |