Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 529,787 | 0 | T | G | 52.4% | 17.5 / 32.3 | 21 | V139G (GTG→GGG) | gcl | glyoxylate carboligase |
Reads supporting (aligned to +/- strand): ref base T (3/7); new base G (11/0); total (14/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.03e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTCCTATTCTGTGCATTACCGGCCAGGCACCGCGCGCCCGTCTGCATAAAGAAGATTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGG‑TTACAGTTCGTGAAGCGGCGCTGGTGCCTCGCGTGCTGCAACAGGCATTTCACCTGATGCGTTCTGGTCGT‑CCGGGTCCGGTACTGGTGGATTTACCGTTCGACGTTCAGGTTGCGGAAATCGAGTTTGATCCTGACATGTACGAACCGCTGCCG > NZ_CP009273/529655‑529917 | ttccTATTCTGTGCATTACCGGCCAGGCACCGCGCGCCCGTCTGCATAAAGAAGATTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGG‑TTACAGTTCGTGAAGCGGCGCTGGTGCCTcg < 2:356073/139‑1 (MQ=255) tgCATTACCGGCCAGGCACCGCGCGCCCGTCTGCATAAAGAAGATTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGG‑TTTTATTTCTTTAAGAGGTGCCTGTGCCCCCCCTGCTGCaa > 1:742313/1‑138 (MQ=255) tgCATTACCGGCCAGGCACCGCGCGCCCGTCTGCATAAAGAAGATTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGG‑GTACAGTTTCTGAAGCGGGGCTGGGGGCTCCCCtgcgccca > 1:437992/1‑133 (MQ=255) cgcgcCCGTCTGCATAAAGAAGATTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGG‑TTACAATTCGGGGAAGGGCGCGGGGGGCCCGCCGGCTGCCACAAGCATTTTCCCTTGTGCGtt > 2:856623/1‑139 (MQ=255) gaagaTTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGG‑TTACAGTTCGTGAAGCGGCGCTGGTGCCTCGCGTGCTGCAACAGGCATTTCACCTGATGCGTTCTGGTCgt‑ccgggtccgg < 2:18183/139‑1 (MQ=255) aTTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGG‑TTTCCGTTTGTGAAACGGCGCTTGGGCGCCTCGGGGTGCCACAAGGATTTTTCCTGATGCGGTTTGGGGGT‑CCGGGGCCGGGACt > 2:121636/1‑139 (MQ=255) aTTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGG‑TTTAAGTTTGTGGAGCGGCGCTGGGGGATCGCGGGCCGCAACAGGCAATTCACCTGCTGGGTTTTTGTCGT‑CCGGGGTCGGtgag > 1:440169/1‑136 (MQ=255) cGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGGTTTAAAGTT‑GTGAAGAGGGGCTGGGGGCCCGCGTGCTGCAACAGGCCTTTCACCTTATGCGTTTTGGTCGT‑TCGGGGCCGGGACTGGGGGGTTTa > 2:731248/1‑139 (MQ=255) tatTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGG‑TTACAGTTCGTGAAGCGGCGCTGGTGCCTCGCGTGCTGCAACAGGCATTTCACCTGATGCGTTCTGGTCGT‑CCGGGTCCGGTACTGGTGGATTTACCGTTc < 2:749015/139‑1 (MQ=255) gatTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGG‑TTACAGTTCGTGAAGCGGCGCTGGTGCCTCGCGTGCTGCAACAGGCATTTCACCTGATGCGTTCTGGTCGT‑CCGGGTCCGGTACTGGTGGATTTACCGTTc < 2:15074/138‑1 (MQ=255) tGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGG‑TTACAATTCGGGGAGGGGGGCCGGGGGCCCCCGGGCTGCAAAAAGCATTTTTCCTGATGCGGTTTGGGCGG‑CCGGGGCCGGGGATGGGGGGTTTTCCGtttggc > 1:436960/1‑135 (MQ=255) aTTGCTAAACCGGTCAGCAAAATGGCGG‑TTACAATTCGTGAAGCGGCGCTGGGGCCTCCCGTGCTGCCACAGGCATTTTACCTGATGCGTTCTGGGCGT‑CCGGGTCCGGGACTGGGGGGTTTACCGTTCGACGTTCAgg > 1:42978/1‑139 (MQ=255) gTCAGCAAAATGGCGG‑TTTCAGTTTGCGAACGGGCTATGGGGGCTCGCGTGCCGGAAAAGGCATTTTACCCCATGCGTGCTTGTTGTGCCGGGCCCCG‑ACTGGTGGAATTTTCGTTCTACGGTCTGGGTGGGGgaaat > 1:40617/1‑136 (MQ=255) gTCAGCAAAATGGCGG‑TTTACGTTCGTGAAGCGGGGCTGGGGGCCCGCGTGGTGCAAAAGGCATTTCTCCTGATGCGTTCTGGTTGT‑CCGGGGCCGGGACTGGGGGGTTTACCGTTTGGCGTTCAGGGTGCGGGAAACg > 1:714825/1‑139 (MQ=255) gTCAGCAAAATGGCGG‑TTACAGTTTGTGAAGCGGCGCTGGGGCCTCGCGGGCTGCAAAAAGCCTTTTCCCCGATGCGGTCTGGTCGT‑CCGGGGGCGGGACTGGGGGGTTTTCCGTTCGGCGGTTAGGGTGCGGaaaaag > 1:766025/1‑136 (MQ=255) gTCAGCAAAATGGCGG‑TTACAGTTCGTGAAGCGGCGCTGGGGGCTCGCGTGCTGCAACAGGCATTTTTCCCGATGCGGTCTGGGCGT‑CCGGGGCCGGGACTGGGGGGTTTACCGTTCGACGTTCAGGGTGCGGGAATCg > 1:197714/1‑139 (MQ=255) cAGTTCGTGAAGCGGCGCTGGTGCCTCGCGTGCTGCAACAGGCATTTCACCTGATGCGTTCTGGTCGT‑CCGGGTCCGGTACTGGTGGATTTACCGTTCGACGTTCAGGTTGCGGAAATCGAGTTTGATCCTGACATGTa < 1:856623/139‑1 (MQ=255) gcggcgCTGGTGCCTCGCGTGCTGCAACAGGCATTTCACCTGATGCGTTCTGGTCGT‑CCGGGTCCGGTACTGGTGGATTTACCGTTCGACGTTCAGGTTGCGGAAATCGAGTTTGATCCTGACATGTACGAACCGCTGc < 2:197714/139‑1 (MQ=255) ggggcTGGTGCCTCGCGTGCTGCAACAGGCATTTCACCTGATGCGTTCTGGTCGT‑CCGGGTCCGGTACTGGTGGATTTACCGTTCGACGTTCAGGTTGCGGAAATCGAGTTTGATCCTGACATGTACGAACCGCTGCCg < 2:98493/136‑1 (MQ=255) ggcgCTGGTGCCTCGCGTGCTGCAACAGGCATTTCACCTGATGCGTTCTGGTCGT‑CCGGGTCCGGTACTGGTGGATTTATCGTTCGACGTTCAGGTTGCGGAAATCGAGTTTGATCCTGACATGTAGGAACCGCAGCCg > 2:287308/1‑139 (MQ=255) ggcgCTGGTGCCTCGCGTGCTGCAACAGGCATTTCACCTGATGCGTTCTGGTCGT‑CCGGGTCCGGTACTGGTGGATTTACCGTTCGACGTTCAGGTTGCGGAAATCGAGTTTGATCCTGACATGTACGAACCGCTGCCg > 2:100945/1‑139 (MQ=255) | TTCCTATTCTGTGCATTACCGGCCAGGCACCGCGCGCCCGTCTGCATAAAGAAGATTTTCAGGCCGTAGATATTGAAGCAATTGCTAAACCGGTCAGCAAAATGGCGG‑TTACAGTTCGTGAAGCGGCGCTGGTGCCTCGCGTGCTGCAACAGGCATTTCACCTGATGCGTTCTGGTCGT‑CCGGGTCCGGTACTGGTGGATTTACCGTTCGACGTTCAGGTTGCGGAAATCGAGTTTGATCCTGACATGTACGAACCGCTGCCG > NZ_CP009273/529655‑529917 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |