| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | REL606 | 3,251,818 | C→T | A315T (GCG→ACG) | nusA ← | transcription elongation factor NusA |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 3,251,818 | 0 | C | T | 94.3% | 119.8 / ‑3.4 | 35 | A315T (GCG→ACG) | nusA | transcription elongation factor NusA |
| Reads supporting (aligned to +/- strand): ref base C (0/2); new base T (16/17); total (16/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.89e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01 | |||||||||||
ATCACGTTGAGTTCCCAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTG > REL606/3251745‑3251889 | atcaCGTTGAGTTCCCAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCa < 1:217571/76‑1 (MQ=255) tcaCGTTGAGTTCCCAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAg > 2:1614583/1‑76 (MQ=255) caCGTTGAGTTCCCAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAgg < 1:671804/76‑1 (MQ=255) cGTTGAGTTCCCAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGtt < 1:783242/76‑1 (MQ=255) ttGAGTTCCCAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTAc < 1:1579518/76‑1 (MQ=255) tGAGTTCCCAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTAcc > 1:595433/1‑76 (MQ=255) cAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCaa > 2:769625/1‑76 (MQ=255) cAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCaa < 1:881901/76‑1 (MQ=255) aCCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACg > 1:784929/1‑76 (MQ=255) ccGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACgg < 2:371440/76‑1 (MQ=255) cGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGc > 2:494632/1‑76 (MQ=255) cAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGata < 1:922756/76‑1 (MQ=255) cTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATcc > 2:901408/1‑76 (MQ=255) cGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATg > 1:942610/1‑76 (MQ=255) aaGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGgt > 2:993617/1‑76 (MQ=255) aGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGgtg > 1:824233/1‑76 (MQ=255) cACGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCt > 2:698451/1‑76 (MQ=255) aCGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCtt < 1:1142195/76‑1 (MQ=255) aCGTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCtt > 1:1559153/1‑76 (MQ=255) gTTCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCa < 2:158747/76‑1 (MQ=255) ttCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGGGTGTTTATCTTGAt > 2:419486/1‑76 (MQ=255) tCTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATc > 1:1299451/1‑76 (MQ=255) cTGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATcc > 1:1291076/1‑76 (MQ=255) tGACCGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATcca < 2:924675/76‑1 (MQ=255) ccGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATccacca < 1:1157100/76‑1 (MQ=255) cGTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATNTTCATCcaccac < 1:20781/76‑1 (MQ=255) gTTACGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACg < 1:1113789/76‑1 (MQ=255) cGGCCAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAg < 1:333231/76‑1 (MQ=255) ccAATCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAg > 1:379655/1‑76 (MQ=255) aaTCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCa > 1:831382/1‑76 (MQ=255) aTCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCaa < 1:695357/76‑1 (MQ=255) aTCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCaa < 1:442565/76‑1 (MQ=255) tCGCCTGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAAc < 1:350869/76‑1 (MQ=255) cTGCGCCAGGTTACCGGTTTCAACGGCGATATCCGTGGTGTGTTTATCTCCATCCACCACGATAGAAGCAACGTCt < 1:177993/76‑1 (MQ=255) tGCGTCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATTCACCACGATAGAAGCAACGTCTg < 2:1491437/76‑1 (MQ=255) | ATCACGTTGAGTTCCCAACCGCTCAGCTGCGAAGCCAGACGCACGTTCTGACCGTTACGGCCAATCGCCTGCGCCAGGTTACCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTG > REL606/3251745‑3251889 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |