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 | 95.0% | 152.6 / ‑4.2 | 40 | I966S (ATC→AGC) | rpoB | DNA‑directed RNA polymerase subunit beta |
Reads supporting (aligned to +/- strand): ref base T (1/1); new base G (22/16); total (23/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 = 9.97e-01 |
CAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCG > REL606/4163681‑4163823 | cAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTc > 2:60888/1‑76 (MQ=255) aaaCGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCg > 1:1470045/1‑76 (MQ=255) aaaCGTGCGATGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCg > 2:1317909/1‑76 (MQ=255) aCGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGaa < 1:522095/76‑1 (MQ=255) cGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAg > 1:818749/1‑76 (MQ=255) gTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGc > 1:1359763/1‑76 (MQ=255) tGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCg > 1:306853/1‑76 (MQ=255) cTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTc < 2:709109/76‑1 (MQ=255) cTGAAAATCAAAGAAAGCCAACTCAAACAGGCGAAGAAAGACCTTTCTGAAGAACTGCAGATCCTCGAAGCGGGTc < 2:709093/76‑1 (MQ=255) tGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCt > 2:741502/1‑76 (MQ=255) tGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCt < 2:1099332/76‑1 (MQ=255) aaTCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTc < 1:240844/76‑1 (MQ=255) aTCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCa < 1:1060089/76‑1 (MQ=255) aTCGAAGAAAGTCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCa < 2:109817/76‑1 (MQ=255) tCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAg < 2:932903/76‑1 (MQ=255) tCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAg < 1:1228453/76‑1 (MQ=255) gaagaaATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTATGTTCAGcc > 1:51724/1‑76 (MQ=255) aaaTGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTAt < 1:1249115/76‑1 (MQ=255) aTGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATcc < 1:604985/76‑1 (MQ=255) aTGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATcc > 2:504135/1‑76 (MQ=255) cTCAAACAGGCGAAGAAAGACCTCTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGCATCCgagctg > 1:577310/1‑76 (MQ=255) aaaCAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtgc > 2:97656/1‑76 (MQ=255) agcAGGAGAAGAAAGACCGATCTGAAGAACCGCAGAGCCGCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtgct < 2:324008/74‑1 (MQ=255) aCAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGtgctgtgctg < 2:914802/76‑1 (MQ=255) cAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGTCGTATCCGTGCTGTGCTgg < 2:330805/76‑1 (MQ=255) cAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTgg > 2:618496/1‑76 (MQ=255) aGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGt < 2:997865/76‑1 (MQ=255) aGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGt < 2:1128095/76‑1 (MQ=255) ggCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTa < 1:1237307/76‑1 (MQ=255) cGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGc < 1:849556/76‑1 (MQ=255) acaaaagtCCTGTCAGAACACCGCCAGACCCCCGACGCGGCTCTGTCCAGCCGTATCCGTGCTGTACTGGTAGCCg < 2:579274/73‑1 (MQ=255) aaagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTg > 2:953966/1‑76 (MQ=255) aaagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTg < 1:1050252/76‑1 (MQ=255) aagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTgg > 2:1019980/1‑76 (MQ=255) aagaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTCCTGGTAGCCGGTgg > 2:237407/1‑76 (MQ=255) gaCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCg > 2:812255/1‑76 (MQ=255) ccTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGtt > 2:159266/1‑76 (MQ=255) cTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTg > 2:446779/1‑76 (MQ=255) gTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGaa > 2:1106885/1‑76 (MQ=255) aggaaCAGAAGAGCCTCGAACGGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTga < 2:162518/74‑1 (MQ=255) aagaaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTga > 2:307247/1‑76 (MQ=255) aagaaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTga > 1:814744/1‑76 (MQ=255) gaaCTCCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTgaga < 2:517175/76‑1 (MQ=255) tcAGAGCCTCAAAGCGGGTGTTTCCAGCCGTATCCGTGCTGTACTGGTAGCCGGTGGCGTTGAAGCCGAGAAGCTc < 2:1313202/75‑1 (MQ=255) cAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCg > 1:267955/1‑76 (MQ=255) cAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCg > 2:264930/1‑76 (MQ=255) | CAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCG > REL606/4163681‑4163823 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |