Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,252,009 | C→A | C251F (TGC→TTC) | nusA ← | transcription elongation factor NusA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,252,009 | 0 | C | A | 96.7% | 107.3 / ‑2.5 | 30 | C251F (TGC→TTC) | nusA | transcription elongation factor NusA |
Reads supporting (aligned to +/- strand): ref base C (1/0); new base A (16/13); total (17/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01 |
AGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGA > REL606/3251937‑3252080 | aGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAg > 2:575920/1‑76 (MQ=255) aCGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGAAAGCAc > 1:299675/1‑76 (MQ=255) ctatCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGAAAGCACCTa < 1:788032/75‑1 (MQ=255) atatCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGAAAGCACCTa > 1:1157530/1‑76 (MQ=255) atatCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGAAAGCACCTa < 1:631957/76‑1 (MQ=255) tatCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGAAAGCACCTAc < 1:559755/76‑1 (MQ=255) atCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGGCACGCATACCGACGAAAGCACCTAcc > 2:957712/1‑76 (MQ=255) atCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCCCGCATACCTACGAAAGCACCTAcc > 2:572362/1‑76 (MQ=255) cGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGAAAGCACCTACCgg > 1:534052/1‑76 (MQ=255) cGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGAAAGCACCTACCGGATCGa < 2:1389347/76‑1 (MQ=255) ccACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGAAAGCACCTACCGGATCGATACGtt < 1:47374/76‑1 (MQ=255) cACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTAACGGATCGAAACGttt > 2:13697/1‑76 (MQ=255) cAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGAAAGCACCTACCGGATCGATACGTTTATCg < 2:1188200/76‑1 (MQ=255) ttCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGAAAGCACCTACCGGATCGATACGTTTATCGTTg < 1:28015/76‑1 (MQ=255) cAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGAAAGCACCTACCGGATCGATACGTTTATCGTTGGt > 2:162469/1‑76 (MQ=255) tAGACACCGCCTGAACACGCGCGCCACGCATACCTACGAAAGCACCTACCGGATCGATACGTTTATCGTTGGtttt < 1:898654/76‑1 (MQ=255) gCCTGAACACGCGCGCCACGCATACCTACGAAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACcgcgc > 1:831689/1‑75 (MQ=255) tGAACACGCGCGCCACGCATACCTACGAAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGCttt > 2:1148346/1‑76 (MQ=255) tGAACACGCGCGCCACGCATACCTACGAAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGAttt > 1:827870/1‑76 (MQ=255) gAACACGCGCGCCACGCATACCTACGAAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGAtttt < 1:633121/76‑1 (MQ=255) gcgcCACGCATACCTACGAAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACgag < 2:324032/76‑1 (MQ=255) cgcCACGCATACCTACGAAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACgaga < 2:1225599/76‑1 (MQ=255) ccACGCATACCTACGAAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAAc > 2:784260/1‑76 (MQ=255) ccACGCATACCTACGAAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAAc < 1:1046073/76‑1 (MQ=255) caggCATACCTACGAAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAAcc > 1:784577/4‑76 (MQ=255) cACGCATACCTACGAAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAAcc > 1:59867/1‑76 (MQ=255) cACGCATACCTACGAAAGCACCTACCGGATCGATACGTTTATCGCTGGTTTTCACCGCGATTTTCGCACGCGAAcc > 2:369094/1‑76 (MQ=255) cGCATACCTACGAAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCg > 2:404834/1‑76 (MQ=255) cATACCTACGAAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGa > 1:895917/1‑76 (MQ=255) cATACCTACGAAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGa > 1:484302/1‑76 (MQ=255) aCCTACGAAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATcg < 2:300657/76‑1 (MQ=255) aCCTACGAAAGCACCTACCGGATAGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATcg > 2:603914/1‑76 (MQ=255) tACGAAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGa > 2:475811/1‑76 (MQ=255) | AGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGA > REL606/3251937‑3252080 |
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: ‑ |