Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,252,012 | G→A | A250V (GCT→GTT) | nusA ← | transcription elongation factor NusA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,252,012 | 0 | G | A | 97.7% | 160.8 / ‑3.1 | 44 | A250V (GCT→GTT) | nusA | transcription elongation factor NusA |
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (22/21); total (22/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01 |
GACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCG > REL606/3251939‑3252084 | gACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACa < 1:1360110/76‑1 (MQ=255) atatCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTa < 2:221203/76‑1 (MQ=255) atatCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTa < 2:523791/76‑1 (MQ=255) cGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCgg > 1:1756851/1‑76 (MQ=255) aTACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGAt < 2:1776841/76‑1 (MQ=255) tACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATc < 1:1401309/76‑1 (MQ=255) aCGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCg < 1:1160488/76‑1 (MQ=255) aCGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCg < 2:1832063/76‑1 (MQ=255) aCGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCg > 2:1441360/1‑76 (MQ=255) aCGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCg > 2:515734/1‑76 (MQ=255) aCGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCg < 2:613514/76‑1 (MQ=255) aCGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCg > 1:426005/1‑76 (MQ=255) cTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATa > 2:1356376/1‑76 (MQ=255) cGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATCCCTACGCAAACACCTACCGGATCGATCCg > 2:969126/1‑76 (MQ=255) cGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACg < 2:1570988/76‑1 (MQ=255) aCCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTa > 2:899767/1‑76 (MQ=255) aCCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTa > 2:354775/1‑76 (MQ=255) cccAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTAt < 2:987226/76‑1 (MQ=255) ccAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATc < 2:434010/76‑1 (MQ=255) cAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCg > 1:700564/1‑76 (MQ=255) aGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGt > 2:1646868/1‑76 (MQ=255) tCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTgg < 2:998990/76‑1 (MQ=255) tCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTgg > 2:1693729/1‑76 (MQ=255) tCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTgg > 2:1567998/1‑76 (MQ=255) tAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGtttt < 2:281308/76‑1 (MQ=255) aGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTc < 2:904732/76‑1 (MQ=255) gACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCa < 1:937833/76‑1 (MQ=255) ccGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACcgc > 2:885637/1‑76 (MQ=255) cGCCTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACcgcg > 2:819288/1‑76 (MQ=255) cTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGAtt > 2:146847/1‑76 (MQ=255) cTGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGAtt > 1:1985969/1‑76 (MQ=255) tGAACACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGAttt > 2:194518/1‑76 (MQ=255) aaCACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTc < 1:1640524/76‑1 (MQ=255) aaCACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTc > 2:265278/1‑76 (MQ=255) aaCACGCGCGCCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTc < 2:827364/76‑1 (MQ=255) gcCACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGaa < 1:48093/76‑1 (MQ=255) cACGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAAcc > 1:1303223/1‑76 (MQ=255) cGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCg > 2:706872/1‑76 (MQ=255) cGCATACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCg < 2:1442598/76‑1 (MQ=255) aTACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGAt < 2:1195101/76‑1 (MQ=255) aTACCTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGAt < 2:1262384/76‑1 (MQ=255) ccTACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATcgc < 1:885637/76‑1 (MQ=255) tACGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGa > 2:1218250/1‑76 (MQ=255) cGCAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGc > 1:940922/1‑76 (MQ=255) cAAACACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCg < 2:779839/76‑1 (MQ=255) | GACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCG > REL606/3251939‑3252084 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |