Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,126,661 | 0 | T | G | 60.0% | 13.7 / 31.0 | 25 | H219P (CAC→CCC) | yghS | dTMP kinase |
Reads supporting (aligned to +/- strand): ref base T (1/9); new base G (15/0); total (16/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.89e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGG > NZ_CP009273/3126573‑3126747 | cGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGgtg < 1:748688/90‑1 (MQ=255) aaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGGTGTCGcc < 1:1473030/90‑1 (MQ=255) tcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGGTGTCGCCCATcc < 1:1200783/90‑1 (MQ=255) tcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGGTGTCGCCCATcc < 1:650676/90‑1 (MQ=255) tcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGGGGTCGCCCCTcc > 1:841673/1‑90 (MQ=255) tcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGGGGGTGCCCCTcc > 2:1535602/1‑90 (MQ=255) cTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGGGGGGGCCCCTCCCGtt > 2:152782/1‑90 (MQ=255) ttttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGGGGGCGCCCCCCCCGTTCAAGGAt > 2:86694/1‑90 (MQ=255) tGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGGTGTCGCCCATcc < 2:956214/76‑1 (MQ=255) tGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGGTGGCGCCCATcc > 1:956214/1‑76 (MQ=255) gTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGGGGGCGCCCATCCAGGTCAAGGATCtttt > 1:636063/1‑89 (MQ=255) gATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGGGGGCGCCCCCCCCCGTCAAAGATCTTTtcg > 2:900523/1‑90 (MQ=255) gATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGGGGGCGCCCCCCCAGTTCTAGGATCTTTtcg > 2:387640/1‑90 (MQ=255) aTCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGGTGTCGCCCATCCAGTTCAAGGATCTTTgcgc < 2:1040309/90‑1 (MQ=255) aGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCAt < 2:1226452/90‑1 (MQ=255) gcgcGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGGGGTCGCCCCTCCAGTTCAAGGGTCTTTTCGCCCTTGAATg > 1:403095/1‑90 (MQ=255) gcgcGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGGGGGCGCCCCCCCAGTTCCAAGAGCTTTGTGCCCTTGaagg > 2:664783/1‑88 (MQ=255) gcggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGgtg > 1:985076/1‑50 (MQ=255) gcggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGgtg < 2:985076/50‑1 (MQ=255) gcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGGGGGTGCCCCTCCAGTTCAAAGGTCTTTTTGCCCTTTGATGTCAGt > 2:303313/1‑90 (MQ=255) gcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGGGGGCGCCCCTCCCGTTCAAAGATCTTTTTGCCCTTTTATGTCAGt > 2:341709/1‑90 (MQ=255) gCTTGCAGAATTTCATCGGCGGGGGGGTCCCCCTCCCCGTTC‑AGGGTCTTTTGGCCACTGAGATTCAGGTGCGGGGTGACAGGGAatttt > 2:1231437/1‑90 (MQ=255) tGCAGAATTTCATCGGC‑GGGGGGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTcc > 1:418108/1‑90 (MQ=255) tGCAGAATTTCATCGGC‑GGGGGGGCGCCCCTCCAGTTTAAGGGTCTTTTGGCCCTTTAATGTCAGTTGGGGGGGAACGGCGGtttttttc > 2:40214/1‑88 (MQ=255) cAGAATTTCATCGGC‑GGGGGGTCGCCCATCCCGTTCAAGGATCTTTGCGGCATTGAATGTCAGTTGGGGGGTAACGGCGGGTTTTTTCTg > 1:1127134/1‑90 (MQ=255) ttgaTCGGC‑GGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGc < 1:311066/87‑1 (MQ=255) ggggTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTtg < 1:40214/90‑1 (MQ=255) gggTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTtgg < 2:841673/90‑1 (MQ=255) | CGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGC‑GGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGG > NZ_CP009273/3126573‑3126747 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |