Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,126,661 | 0 | T | G | 55.6% | 12.2 / 24.6 | 18 | H219P (CAC→CCC) | yghS | dTMP kinase |
Reads supporting (aligned to +/- strand): ref base T (0/8); new base G (10/0); total (10/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.29e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGG > NZ_CP009273/3126573‑3126747 | cGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgtg < 2:114425/90‑1 (MQ=255) tcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATcc < 1:217487/90‑1 (MQ=255) aTCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCCCCCGTTCa > 2:682669/1‑90 (MQ=255) aTTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAGTTCAAggg > 1:136743/1‑89 (MQ=255) ttttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGcc < 2:39615/74‑1 (MQ=255) ttttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGcc > 1:39615/1‑74 (MQ=255) gATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTtcg > 1:98200/1‑90 (MQ=255) gcgcGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATg < 2:22222/90‑1 (MQ=255) gcgcGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGGGCCCCCCCCGTTTAAGGATCTTTTCGCCCTTTAATg > 2:213410/1‑90 (MQ=255) gcgcGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGGGCCCCCCCCCTTCAAGGAGTTTTTCGCCCTTGTATg > 2:364461/1‑90 (MQ=255) gcggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTc < 1:389657/90‑1 (MQ=255) acgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGGGCCCCTCCCGTTCAAAGATCTTTTCGCCATTTAATGTCAGt > 2:725404/2‑90 (MQ=255) tGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTgc < 2:83737/90‑1 (MQ=255) gacagTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTa < 1:338774/86‑1 (MQ=255) tAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCGAGGATCTTTTCGCCATTGAATGTCAGTTGCGGCGGAACGGc > 1:178108/1‑90 (MQ=255) tGCAGAATTTCATCGGCGGGGGGGGGCCCCTCCCGGTTAAGGGTCTTTTCGCCCTTTAATGTTAGGTTGGGCGGGAAGGCGAtttttttt > 2:338774/1‑88 (MQ=255) gCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTTTGCCCTTGAATGTTAGTTGGGGGGGAAAGGCGGTTTTTTTCt > 2:208793/1‑90 (MQ=255) gggTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTtgg < 1:592351/90‑1 (MQ=255) gggTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTtgg < 1:472637/90‑1 (MQ=255) | CGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGG > 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 < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |