Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_017635 | 3,668,239 | 0 | G | A | 22.9% | 100.2 / 31.2 | 48 | G318G (GGC→GGT) | tuf | elongation factor Tu |
Reads supporting (aligned to +/- strand): ref base G (23/14); new base A (9/2); total (32/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.93e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.68e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TACGGAAGTAGAACTGCGGACGGTAGCCTTTGAAGAACGGAGTATGACGGCCGCCTTCATCTTTGGACAGAATGTACACTTCAGATTCGAACTTGG > NC_017635/3668190‑3668285 | taCGGAAGTAGAACTGCGGACGGTAGCCTTTGAAGAACGGAGTATGACgg > 1:1843159/1‑50 (MQ=32) cGGAAGTAGAACTGCGGACGGTAGCCTTTGAAGAACGGAGTATGACGgcc < 1:184581/50‑1 (MQ=38) ggAAGTAGAACTGCGGACGGTAGCCTTTGAAGAACGGAGTATGACGaccg < 1:1455789/50‑1 (MQ=12) ggAAGTAGAACTGCGGACGGTAGCCTTTGAAGAACGGAGTATGACGaccg > 2:472474/1‑50 (MQ=12) gAAGTAGAACTGCGGACGGCAGCCTTTGAAGAACGGAGTATGACGaccgc > 2:411231/1‑50 (MQ=11) aaGTAGAACTGCGGACGGTAGCCTTTGAAGAACGGAGTATGACGgccgcc < 1:355982/50‑1 (MQ=255) aaGTAGAAATGCGGACGGTAGCCTTTGAAGAACGGAGTATGACGaccgcc < 1:912143/50‑1 (MQ=11) aGTAGAACTGCGGACGGTAGCCTTTGAAGAACGGAGTATGACGGCCGCCt > 2:842947/1‑50 (MQ=255) tAGAACTGCGGACGGTAGCCTTTGAAGAACGGAGTATGACGGCCGCCtt > 2:611275/1‑49 (MQ=255) tAGAACTGCGGACGGTAGCCTTTGAAGAACGGAGTATGACGGCCGCCTTc < 1:2061334/50‑1 (MQ=255) tAGAACTGCGGACGGTAGCCTTTGAAGAACGGAGTATGACGACCGCCtt > 1:1667135/1‑49 (MQ=255) tAGAACTGCGGACGGTAGCCTTTGAAGAACGGAGTATGACGACCGCCtt > 2:1791923/1‑49 (MQ=255) gAACTGCGGACGGTAGCCTTTGAAGAACGGAGTATGACGGCCGCCTTCAt < 1:195980/50‑1 (MQ=255) aaCTGCGGACGGTAGCCTTTGAAGAACGGAGTATGACGGCCGCCTTCATc < 2:233244/50‑1 (MQ=255) aCTGCGGACGGTAGCCTTTGAAGAACGGAGTATGACGGCCGCCTTCATCt < 1:250908/50‑1 (MQ=255) aCTGCGGACGGTAGCCTTTGAAGAACGGAGTATGACGGCCGCCTTCATCt > 2:1288782/1‑50 (MQ=255) cTGCGGACGGTAGCCTTTGAAGAACGGAGTATGACGACCGCCTTCATCtt > 2:998020/1‑50 (MQ=255) cTGCGGACGGTAGCCTTTGAAGAACGGAGTATGACGACCGCCTTCATCtt > 2:964483/1‑50 (MQ=255) cGGACGGTAGCCTTTGAAGAACGGAGTATGACGGCCGCCTTCATCTTTg < 2:906868/49‑1 (MQ=255) cGGACGGTAGCCTTTGAAGAACGGAGTATGACGGCCGCCTTCATCTTTg < 1:394682/49‑1 (MQ=255) ggACGGTAGCCTTTGAAGAACGGAGTATGACGGCCGCCTTCATCTTTGGa > 1:1961483/1‑50 (MQ=255) ggACGGTAGCCTTTGAAGAACGGAGTATGACGGCCGCCTTCATCTTTGGa < 1:860873/50‑1 (MQ=255) gACGGTAGCCTTTGAAGAACGGAGTATGACGACCGCCTTCATCTTTGGa > 1:553308/1‑49 (MQ=255) gACGGTAGCCTTTGAAGAACGGAGTATGACGACCGCCTTCATCTTTGGa > 1:1325814/1‑49 (MQ=255) aCGGTAGCCTTTGAAGAACGGAGTATGACGGCCGCCTTCATCTTTGGa > 1:396694/1‑48 (MQ=255) aCGGTAGCCTTTGAAGAACGGAGTATGACGGCCGCCTTCATCTTTGGa > 1:799167/1‑48 (MQ=255) cGGTAGCCTTTGAAGAACGGAGTATGACGGCCGCCTTCATCTTTGGa > 2:696288/1‑47 (MQ=255) cGGTAGCCTTTGAAGAACGGAGTATGACGGCCGCCTTCATCTTTGGACAg < 2:547256/50‑1 (MQ=255) ggTAGCCTTTGAAGAACGGAGTATGACGGCCGCCTTCATCTTTGGACAg > 1:547135/1‑49 (MQ=255) ggTAGCCTTTGAAGAACGGAGTATGACGGCCGCCTTCATCTTTGGACAg > 2:506932/1‑49 (MQ=255) ggTAGCCTTTGAAGAACGGAGTATGACGGCCGCCTTCATCTTTGGACAg > 2:1996825/1‑49 (MQ=255) ggTAGCCTTTGAAGAACGGAGTATGACGACCGCCTTCATCTTTGGACAGa < 2:980307/50‑1 (MQ=255) gTAGCCTTTGAAGAACGGAGTATGACGGCCGCCTTCATCTTTGGACAGaa > 2:807001/1‑50 (MQ=255) gCCTTTGAAGAACGGAGTATGACGGCCGCCTTCATCTTTGGACAGAATg < 1:2060931/49‑1 (MQ=255) cTTTGAAGAACGGAGTATGACGGCCGCCTTCATCTTTGGACAGAATGTa > 1:161025/1‑49 (MQ=255) tttGAAGAACGGAGTATGACGGCCGCCTTCATCTTTGGACAGAATGTaca > 2:106442/1‑50 (MQ=255) ttGAAGAACGGAGTATGACGACCGCCTTCATCTTTGGACAGAATGTacac > 2:1725945/1‑50 (MQ=255) tGAAGAACGGAGTATGACGGCCGCCTTCATCTTTGGACAGAATGTaca < 1:1708351/48‑1 (MQ=255) tGAAGAACGGAGTATGACGGCCGCCTTCATCTTTGGACAGAATGTACACt > 1:1223431/1‑50 (MQ=255) tGAAGAACGGAGTATGACGGCCGCCTTCATCTTTGGACAGAATGTACACt > 1:264940/1‑50 (MQ=255) gaagaaCGGAGTATGACGGCCGCCTTCATCTTTGGACAGAATGTACACt < 1:484918/49‑1 (MQ=255) aagaaCGGAGTATGACGGCCGCCTTCATCTTTGGACAGAATGTACACTTc < 1:149517/50‑1 (MQ=255) aagaaCGGAGTATGACGACCGCCTTCATCTTTGGACAGAATGTACACtt > 2:697029/1‑49 (MQ=255) agaaCGGAGTATGACGGCCGCCTTCATCTTTGGACAGAATGTACACtt > 1:1650148/1‑48 (MQ=255) aaCGGAGTATGACGACCGCCTTCATCTTTGGACAGAATTTACACtt > 2:1999854/1‑46 (MQ=11) aCGGAGTATGACGACCGCCTTCATCTTTGGACAGAATGTACACTTCAGAt < 1:1108687/50‑1 (MQ=255) cGGAGTATGACGGCCGCCTTCATCTTTGGACAGAATGTACACTTCAGAtt < 1:325903/50‑1 (MQ=255) ggAGTATGACGGCCGCCTTCATCTTTGGACAGAATGTACACTTCAGATTc > 1:1672399/1‑50 (MQ=255) ggAGTATGACGGCCGCCTTCATCTTTGGACAGAATGTACACTTCAGATTc > 2:287934/1‑50 (MQ=255) gAGTATGACGGCCGCCTTCATCTTTGGACAGAATGTACACTTCAGATTc > 2:1204466/1‑49 (MQ=255) gAGTATGACGGCCGCCTTCATCTTTGGACAGAATGTACACTTCAGATTCg < 1:43938/50‑1 (MQ=255) aGTATGACGGCCGCCTTCATCTTTGGACAGAATGTACACTTCAGATTCg > 1:97395/1‑49 (MQ=255) aGTATGACGGCCGCCTTCATCTTTGGACAGAATGTACACTTCAGATTCg > 1:580549/1‑49 (MQ=255) aTGACGGCCGCCTTCATCTTTGGACAGAATGTACACTTCAGATTCGAACt > 1:546039/1‑50 (MQ=255) gacgacCGCCTTCATCTTTGGACAGAATGTACACTTCAGATTCGAACTTg > 1:1967481/1‑50 (MQ=255) aCGGCCGCCTTCATCTTTGGACAGAATGTACACTTCAGATTCGAACTTgg > 1:231737/1‑50 (MQ=38) | TACGGAAGTAGAACTGCGGACGGTAGCCTTTGAAGAACGGAGTATGACGGCCGCCTTCATCTTTGGACAGAATGTACACTTCAGATTCGAACTTGG > NC_017635/3668190‑3668285 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |