Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 4,164,097 | C→T | P1081L (CCG→CTG) | 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,164,097 | 0 | C | T | 97.4% | 133.8 / ‑3.0 | 38 | P1081L (CCG→CTG) | rpoB | DNA‑directed RNA polymerase subunit beta |
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (17/20); total (17/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01 |
AAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCG > REL606/4164024‑4164167 | aaaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGa < 2:797954/76‑1 (MQ=255) aCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATc < 1:45018/76‑1 (MQ=255) aTCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGAGCAATTTCTAAGATCAACCAGATCGAAGata > 1:1276951/1‑76 (MQ=255) tCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGatat < 1:1112154/76‑1 (MQ=255) ccAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATg > 1:132165/1‑76 (MQ=255) ccAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATg > 2:1196233/1‑76 (MQ=255) catccTTGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGc < 2:130103/73‑1 (MQ=255) tcTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCtt < 2:677890/75‑1 (MQ=255) gcAAGGTGGCAGGTGGTCGCGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGa < 1:1415273/75‑1 (MQ=255) aCAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGa > 2:179895/1‑76 (MQ=255) aGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGaaaa < 1:661877/76‑1 (MQ=255) tGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACgg > 2:1313180/1‑76 (MQ=255) tGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACgg < 2:843379/76‑1 (MQ=255) ggCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGt > 1:1505502/1‑76 (MQ=255) gCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTa < 1:164251/76‑1 (MQ=255) cAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTAc < 2:760161/76‑1 (MQ=255) cAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTAc > 2:912061/1‑76 (MQ=255) cAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTAc > 1:48708/1‑76 (MQ=255) aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACg < 1:620565/76‑1 (MQ=255) aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACg < 2:29378/76‑1 (MQ=255) cgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCgg > 1:379060/1‑76 (MQ=255) gtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGt < 1:390565/76‑1 (MQ=255) cGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGAc < 1:1615766/76‑1 (MQ=255) tAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATc < 1:1017685/76‑1 (MQ=255) aaCAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCg > 1:1351305/1‑76 (MQ=255) cAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTa > 2:1222920/1‑76 (MQ=255) gggTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTg < 1:455445/76‑1 (MQ=255) gggTGTAATTTCTAAGATCAACCTGATCGAAGAGATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTg < 2:1437926/76‑1 (MQ=255) gtAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAAcc > 1:1160060/1‑76 (MQ=255) gtAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAAcc < 2:1396511/76‑1 (MQ=255) gtAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAAcc > 1:1456262/1‑76 (MQ=255) aTTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGc < 1:686212/76‑1 (MQ=255) tttCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCt > 2:1055235/1‑76 (MQ=255) tCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTgg > 1:5498/1‑76 (MQ=255) cTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTggg < 1:1036338/76‑1 (MQ=255) tAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGc > 1:1600889/1‑76 (MQ=255) tAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGc > 2:1573876/1‑76 (MQ=255) aGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGt < 2:1005173/76‑1 (MQ=255) gATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTa < 1:821930/76‑1 (MQ=255) aTCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTAc > 1:77367/1‑76 (MQ=255) cAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCg < 2:366542/76‑1 (MQ=255) | AAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCG > REL606/4164024‑4164167 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |