| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | REL606 | 4,163,752 | T→G | I966S (ATC→AGC) | rpoB → | DNA‑directed RNA polymerase subunit beta |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 4,163,752 | 0 | T | G | 92.1% | 137.4 / 1.9 | 38 | I966S (ATC→AGC) | rpoB | DNA‑directed RNA polymerase subunit beta |
| Reads supporting (aligned to +/- strand): ref base T (3/0); new base G (27/8); total (30/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01 | |||||||||||
GACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACA > REL606/4163679‑4163826 | gACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGcc > 2:880503/1‑76 (MQ=255) aCAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCt < 2:477782/76‑1 (MQ=255) cAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTc > 1:1323731/1‑76 (MQ=255) aaaCGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCg > 2:58905/1‑76 (MQ=255) aCGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGaa < 1:1692931/76‑1 (MQ=255) cTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTc > 2:1656091/1‑76 (MQ=255) aaTCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTc < 2:528605/76‑1 (MQ=255) tCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGACTGTTCAg > 2:116421/1‑76 (MQ=255) cGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGc > 1:1395492/1‑76 (MQ=255) gaagaaATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGcc < 1:1321974/76‑1 (MQ=255) aaaTGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTAt > 1:537432/1‑76 (MQ=255) aaTGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATc > 2:1572380/1‑76 (MQ=255) cTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctg > 1:714373/1‑76 (MQ=255) aaaCAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtgc < 1:430991/76‑1 (MQ=255) aCAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGtgctgtgctg > 2:1743711/1‑76 (MQ=255) aCAGGCGAACAAATACCTTTCCGAAAAAGGGGAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGtgctgtgctg < 2:226972/76‑1 (MQ=255) cGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGc > 1:1593892/1‑76 (MQ=255) cGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGc < 2:1421450/76‑1 (MQ=255) gaagaaAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGcc < 2:1457575/76‑1 (MQ=255) aagaaagaCCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCg > 2:267743/1‑76 (MQ=255) aaaaaaggCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCg < 2:377097/73‑1 (MQ=255) aCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGt > 1:87077/1‑76 (MQ=255) aCCTGTCGGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGt < 2:1565750/76‑1 (MQ=255) ccTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGtt > 1:363395/1‑76 (MQ=255) cTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTg > 1:935744/1‑76 (MQ=255) tGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGa > 1:1421498/1‑76 (MQ=255) gTCTGAAGAACTGCAGAGCCTCGAAGCGGGGCTGTTCAGCCGTATCCGCGCCGTGCTGGTAGCCGGAGGCGTTGaa > 2:472745/1‑76 (MQ=255) tCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAg > 2:1408358/1‑76 (MQ=255) tCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAg > 2:774449/1‑76 (MQ=255) tCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGACGGTGGCGTTGAAg > 2:1720595/1‑76 (MQ=255) cTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGc > 1:1522605/1‑76 (MQ=255) tGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCt > 2:1127706/1‑76 (MQ=255) tGAAGAACTGCAGAGCCGCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCt < 2:136227/76‑1 (MQ=255) gaagaaCTCCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTg < 2:582676/76‑1 (MQ=255) aagaaCTGCAGCGCCCCGAAGCGGGTCTGTCCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTga < 2:1177862/76‑1 (MQ=255) aagaaCTGCAGAGCCTCGAAGGGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTga < 2:1234225/76‑1 (MQ=255) aagaaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTga > 2:301115/1‑76 (MQ=255) gaaCTGCAGGGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTgaga < 1:942322/76‑1 (MQ=255) gaaCTGCAGAGCCTCGAAGCGGGTTTGTTCATCCGTAGCCGTGCTGTGCTGGTAGCCGGTGGCGTGGAGGCGgaga > 2:98786/1‑76 (MQ=255) aaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGaa > 2:51444/1‑76 (MQ=255) aaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGaa > 2:1606482/1‑76 (MQ=255) aCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAg > 2:519541/1‑76 (MQ=255) tGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCt > 2:841435/1‑76 (MQ=255) ggCAGAGCCGCGAAGGGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCt < 2:147623/75‑1 (MQ=255) cAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCg > 2:541661/1‑76 (MQ=255) agagCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGa > 2:558801/1‑76 (MQ=255) aTCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGAGCTGGTAGCCCGGGGGGGTGACGGTGAGAAGCTCcaca > 1:928604/1‑76 (MQ=255) | GACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACA > REL606/4163679‑4163826 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |