| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | REL606 | 4,163,089 | A→G | E745G (GAA→GGA) | 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,089 | 0 | A | G | 94.4% | 127.6 / ‑4.0 | 36 | E745G (GAA→GGA) | rpoB | DNA‑directed RNA polymerase subunit beta |
| Reads supporting (aligned to +/- strand): ref base A (0/2); new base G (13/21); total (13/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.25e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGA > REL606/4163016‑4163157 | ggtggtgTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAg < 2:800106/76‑1 (MQ=255) ggtgTCGTTCAGTACGTGGATGCTTCCNATATCGTNNTCAAAGTTAACGAAGACGAGTTGTATCCTGGGGGAGCAg > 1:826456/1‑76 (MQ=255) gtgtCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAgg < 1:574634/76‑1 (MQ=255) gtCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGTa > 2:771056/1‑76 (MQ=255) gtCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGTa > 1:561486/1‑76 (MQ=255) cGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGTATc < 2:360474/76‑1 (MQ=255) cGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGTATc > 2:191908/1‑76 (MQ=255) tCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGTATCGAc > 2:340190/1‑76 (MQ=255) aGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCAGGGTGAAGCAGGTATCGACAt < 1:421597/76‑1 (MQ=255) gTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGTATCGACATc < 1:442920/76‑1 (MQ=255) tACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGTATCGACATCt > 2:24815/1‑76 (MQ=255) aCGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGTATCGACATCTa < 1:373468/76‑1 (MQ=255) aCGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGTATCGACATCTa > 2:503092/1‑76 (MQ=255) tGGATGCTTCCCGTATCGTAATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGTATCGACATCTACaa > 1:716172/1‑76 (MQ=255) ggATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGTATCGACATCTACAAc < 1:199792/76‑1 (MQ=255) ggATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGTATCGACATCTACAAc < 2:183314/76‑1 (MQ=255) cTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATACGGGTGGAGCAGGTATCGACATCTACAACCTGAc < 1:258134/76‑1 (MQ=255) ttCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGTATCGACATCTACAACCTGAcc < 1:407537/76‑1 (MQ=255) cccGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAACACGTATCGACATCTACAACCTGACCac > 1:129630/1‑75 (MQ=255) ccGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCGGGGTGGAGCAGGTATCGACATCTACAACCTGACCaaa < 1:706586/76‑1 (MQ=255) ccGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGTATCGACATCTACAACCTGACCaaa > 2:391684/1‑76 (MQ=255) gTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAAGTCTCGACATCTACAACCTGACCAAAtc > 1:126207/1‑75 (MQ=255) tATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGTATCGACATCTACAACCTGACCAAATac < 1:159682/76‑1 (MQ=255) tATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGTATCGACATCTACAACCTGACCAAATac < 1:492115/76‑1 (MQ=255) tATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGNATCGACATCTACAACCTGACCAAATac < 2:1426/76‑1 (MQ=255) tATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGANATCTACATCCNGACCAAATac < 2:830025/76‑1 (MQ=255) gTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGTATCGACATCTACAACCTGACCAAATACAccc < 2:259029/76‑1 (MQ=255) ttATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGTATCGACATCTACAACCTGACCAAATACACCCg > 2:289881/1‑76 (MQ=255) tATCAAAGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGt > 1:375889/1‑76 (MQ=255) aaaGTTAACGAAGACGAGATGTATCCGGGTGGAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTa < 2:561486/76‑1 (MQ=255) cGAAGACGAGATGTATCCGGGTGGAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAac > 2:440611/1‑76 (MQ=255) gACGAGATGTATCCGGGTGGAGCAGGTATCGACATCAACTACCTGACCNAATACACCCGTTCTAACCAGAACACCt < 1:739415/76‑1 (MQ=255) aCGAGATGTATCCGGGTGGAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTg > 1:661725/1‑76 (MQ=255) cGAGATGTATCCGGGTGGAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGt < 2:672774/76‑1 (MQ=255) aTGTATCCGGGTGGAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCa < 2:548509/76‑1 (MQ=255) tGTATCCGGGTGGAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCaa < 2:661725/76‑1 (MQ=255) gTATCCGGGTGGAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACTAGAACACCTGTATCAAc < 1:66766/76‑1 (MQ=255) ccGGGTGGAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGa < 2:281577/76‑1 (MQ=255) | GGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGA > REL606/4163016‑4163157 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |