Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 4,166,437 | C→T | P493S (CCG→TCG) | rpoC → | DNA‑directed RNA polymerase subunit beta' |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,166,437 | 0 | C | T | 97.1% | 128.3 / NA | 35 | P493S (CCG→TCG) | rpoC | DNA‑directed RNA polymerase subunit beta' |
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (18/16); total (18/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.86e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGAC > REL606/4166364‑4166508 | tCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCTcg < 1:464720/76‑1 (MQ=255) cACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCTcgg > 2:514130/1‑76 (MQ=255) cACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCTcgg < 1:268315/76‑1 (MQ=255) cGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCcggcg < 1:78355/76‑1 (MQ=255) gTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTAGCAACAACATCCTGTCCTcggcga < 1:660134/76‑1 (MQ=255) tACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCTcggcgaa > 1:526169/1‑76 (MQ=255) cGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCTcggcgaacgg > 2:549976/1‑76 (MQ=255) gACGCTGGAAGCCCAGCTAGANGCNCGTGCGCTGATGATGTCTACCNACAACGTNCNGTCCTcggcgaacggcgaa > 2:176691/1‑76 (MQ=255) gCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCTCGGCGAACGGCGAACCa > 1:425662/1‑76 (MQ=255) ggAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCTCGGCGAACGGCGAACCAatc < 2:1169805/76‑1 (MQ=255) gAAGCTGAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCTCGGCGAACGGCGAACCAatca < 2:56527/76‑1 (MQ=255) aGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCTCGGCGAACGGCGAACCAATCATCGTTcc > 2:948059/1‑76 (MQ=255) gCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCTCGGCGAACGGCGAACCAATCATCGTTCCg < 2:517366/76‑1 (MQ=255) cTGGAAGCGCATGCGCTCATGATGTCTACCAACAACATCCTGTCCTCGGCGAACGGCGAACGAATCATCGTTCCGt > 1:97424/1‑76 (MQ=255) gAAGCGCGTGCGCTTATGATGTCTACCAACAACATCCTGTCCTCGGCGAACGGCGAACGAATCATAGTTCCGtctc < 1:101667/76‑1 (MQ=255) gAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCTCGGCGAACGGCGAACCAATCATCGTTCCGtctc < 1:850114/76‑1 (MQ=255) gAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCTCGGCGAACGGCGAACCAATCATCGTTCCGtctc > 2:442198/1‑76 (MQ=255) aaGCGCGTGCNCGAATGATGTNTACCAACAACATGCTGTCCTCGGCGAACGGCGAACCAATCATCGTTCCGTCTCa > 2:1189386/1‑76 (MQ=255) aaGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCTCGGCGAACGGCGAACCAATCATCGTTCCGTCTCa < 2:120102/76‑1 (MQ=255) aaGCGCGTGCGCTGATGATGTCTACCAACAACAGCCTGTCCTCGGCGAACGGCGAACCAATCATCGTTCCGTCTCa < 2:307918/76‑1 (MQ=255) aGCGCGTGCGCTGATGATGTCTACCGACGACATCCTGTCCTCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAg < 1:190130/76‑1 (MQ=255) aGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCTCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAg < 2:1039296/76‑1 (MQ=255) cgTGCGCTGATGATGTCTACCAACAACATCCTGTCCTCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACg > 2:323346/1‑76 (MQ=255) cgTGCGCTGATGATGTCTACCAACAACATCCTGTCCTCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACg > 2:428310/1‑76 (MQ=255) gTGCGCTGATGATGTCTACCAACAACATCCTGTCCTCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGt > 2:442808/1‑76 (MQ=255) cgcTGATGATGTCTACCAACAACATCCTGTCCTCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGt < 1:729496/76‑1 (MQ=255) cTGATAATGTCTACCAACAACATCCTGTCCTCGGCGAACGCCGAAGCAATCATCGTTCCCTCTCAGGACGGCGTAc > 1:76908/1‑76 (MQ=255) tgatgaTGTCTACCAACAACATCCTGTCCTCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACt < 2:1112884/76‑1 (MQ=255) tgatgaTGTCTACCAACAACATCCTGTCCTCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACt > 1:746031/1‑76 (MQ=255) atgatgTCTACCAACAACATACTGTCCTCGGCGAATGGCTAACCAATCATCGTTCCGTCTCAGGACGTTGTACTgg < 2:965252/76‑1 (MQ=255) gTCTACCAACAACATCCTGTCCTCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTg < 2:526169/76‑1 (MQ=255) tACCAACAACATCCTGTCCTCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGtac > 2:819521/1‑76 (MQ=255) aCCAACAACATCCTGTCCTCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGtact > 2:278961/1‑76 (MQ=255) aCCAACAACATCCTGTCCTCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGtact > 2:595888/1‑76 (MQ=255) caacaaCATCCTGTCCTCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGtactac < 2:111354/76‑1 (MQ=255) ccTGTCCTCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGt > 1:890740/1‑76 (MQ=255) cTGTCCTCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTg > 2:1092576/1‑76 (MQ=255) gTCCTCGGCGAAAGGCGAACCAATCATCGTTCCATCTCAGGACGTTGTACTGGGTCTGTACNACATGACCCGTGAc > 1:1039431/1‑76 (MQ=255) | TCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGAC > REL606/4166364‑4166508 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |