Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,992,510 | 0 | G | A | 43.3% | 25.6 / 19.6 | 30 | D390N (GAT→AAT) | uvrD | DNA helicase II |
Reads supporting (aligned to +/- strand): ref base G (4/13); new base A (13/0); total (17/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.17e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTGCTCGAAGAGGCGTTATTGCAGGCCAGTATGCCGTACCGTATTTACGGCGGGATGCGCTTCTTC‑GAACGCCAGGAAATCAAAGATGCGCTCTCGTATCTGCGCCTGATTGCCAACCGCAACGACGACGCGGCCTTTGAGCGTGTGGTGAATACGCCAACGCGG > NZ_CP009273/3992426‑3992590 | gTGCTCGAAGAGGCGTTATTGCAGGCCAGTATGCCGTACCGTATTTACGGCGGGATGCGCTTCTTC‑GAACGCCAGGAAATCAAAGATgcg < 2:507163/90‑1 (MQ=255) gCTCGAAGAGGCGTTATTGCAGGCCAGTATGCCGTACCGTATTTACGGCGGGATGCGCTTCTTC‑GAACGCCAGGAAATCAAAGATGCGct < 2:84689/90‑1 (MQ=255) cTCGAAGAGGCGTTATTGCAGGCCAGTATGCCGTACCGTATTTACGGCGGGGGGGGCTTCTTC‑GTACGCCCGGAAAAAAAAAAAGCGctc > 2:107424/1‑90 (MQ=255) aGGCCAGTATGCCGTACCGTATTTACGGCGGGTGGGGCTTTTTCGGAACCCCCGGAAA‑CAAAGAGGTGCGCTAGCATCTTGGCCTGTTTg > 2:352674/1‑90 (MQ=255) aGGCCAGTATGCCGTACCGTATTTACGGCGGGGGGGGCTTTTTT‑TAACGCCCGGAAAACAAAAAAGGGCTCTCGTATCTGCGCCCGAttt > 2:576352/1‑89 (MQ=255) gTATGCCGTACCGTATTTACGGCGGGGTGCGCTTCTTC‑GAACGCCAGGGAAAAAAAAAAGCGCCCTCCTATCTGCGCCTGATTGCCCAcc > 1:333889/1‑90 (MQ=255) gTATGCCGTACCGTATTTACGGCGGGATGGGCCTCTTT‑TAACGCCCGGAAAACAAAAAAGGGCCCTCGGATTTGCGCCCGCTTTCCAAcc > 1:419133/1‑90 (MQ=255) ccgtaccgtaTTTACGGCGGGAGGCGCTTCTTC‑GCACGCCAGGAAAACAAAAATGGGGTCTCGTATCTGCGCCTGATTTCCAACCCCAac > 2:325062/1‑90 (MQ=255) gtaccgtaTTTACGGCGGGGTGCGCTTCTTT‑TCACCCCCCGAAAACAAAAAAGCGCTCTCCCATCTTCGCCCGGTTTCCCAACGCAacga > 2:329195/1‑90 (MQ=255) ccgtaTTTACGGCGGGGTGGGCTTCTTC‑GTACGCCAGGAAAACAAAAATGCGCTCTCGTATCTTCGCCTGGTTGCCAACCGCCaccacga > 2:258534/1‑90 (MQ=255) gtaTTTACGGCGGGGTGGGCTTCTTT‑GAACGCCAGGAAAACAAAAAAGGGCTCTCGTATCTGCGCCTGATTGCCCACCGCAacgacgacg > 2:137154/1‑90 (MQ=255) gtaTTTACGGCGGGGTGCGCTTTTTT‑GAACGCCACGAAAAAAAAAAAGGGCTCTCCTGTCTTCCCCCGGGTTGCAACCGCAaccacgacg > 2:461622/1‑90 (MQ=255) gtaTTTACGGCGGGGTGCGCTTCTTT‑GAACGCCCGGGAAACAAAAAAGCGCTCTCGTATCTTCGCCTCATTTCCCACCCCAacgacgacg > 2:18594/1‑90 (MQ=255) gtaTTTACGGCGGGATGCGCTTCTTC‑GAACGCCCGGAAAACAAAAATGCGCTCTCGTATCTGCGCCCGATTGCCAACCGCAacgactgtc > 1:78612/1‑86 (MQ=255) gtaTTTACGGCGGGAAGCGCTTCTTC‑TTACGCCCGGAAAACAAAAAAGCGCTCTCCTATCTGCGCCTGTTTTCCAACCGCAacgacgacg > 2:126450/1‑90 (MQ=255) taTTTACGGCGGGATGCGCTTCTTC‑GAACGCCAGGAAATCAAAGATGCGCTCTCGTATCTGCGCCTGATTGCCAACCGCAACGACGAcgc < 1:576352/90‑1 (MQ=255) taTTTACGGCGGGATGCGCTTCTTC‑GAACGCCAGGAAATCAAAGATGCGCTCTCGTATCTGCGCCTGATTGCCAACCGCAACGACGAcgc < 1:264883/90‑1 (MQ=255) tttACGGCGGGATGCGCTTCTTC‑GAACGCCAGGAAATCAAAGATGCGCTCTCGTATCTGCGCCTGATTGCCAACCGCAACGACGACGCgg < 1:461622/90‑1 (MQ=255) cggGATGCGCTTCTTC‑GAACGCCAGGAAATCAAAAATGCGCTCTCGTATCTGCGCCTGATTGCCAACCGCAACGACGACGCGGCCtttta > 1:245745/1‑88 (MQ=255) ggATGCGCTTCTTC‑GAACGCCAGGAAATCAAAGATGCGCTCTCGTATCTGCGCCTGATTGCCAACCGCAACGACGACGCGGCCTTTGAGc < 1:352674/90‑1 (MQ=255) gcTTCTTC‑GAACGCCAGGAAATCAAAGATGCGCTCTCGTATCTGCGCCTGATTGCCAACCGCAACGACGACGCGGCCTTTGAGCGTgtgg < 1:42914/90‑1 (MQ=255) ttc‑GAACGCCAGGAAATCAAAGATGCGCTCTCGTATCTGCGCCTGATTGCCAACCGCAACGACGACGCGGCCTTTGAGCGTGTGGTGAAt < 2:611749/90‑1 (MQ=255) tc‑GAACGCCAGGAAATCAAAGATGCGCTCTCGTATCTGCGCCTGATTGCCAACCGCAACGACGACGCGGCCTTTGAGCGTGTGGTGAATa > 1:34840/1‑90 (MQ=255) ccAGGAAATCAAAGATGCGCTCTCGTATCTGCGCCTGATTGCCAACCGCAACGACGACGCGGCCTTTGAGCGTGTGGTGAATAc < 1:562837/84‑1 (MQ=255) ccAGGAAATCAAAGATGCGCTCTCGTATCTGCGCCTGATTGCCAACCGCAACGACGACGCGGCCTTTGAGCGTGTGGTGAATAc > 2:562837/1‑84 (MQ=255) gAAATCAAAGATGCGCTCTCGTATCTGCGCCTGATTGCCAACCGCAACGACGACGCGGCCTTTGAGCGTGTGGTGAATACGCCAACGCgg < 1:137154/90‑1 (MQ=255) gAAATCAAAGATGCGCTCTCGTATCTGCGCCTGATTGCCAACCGCAACGACGACGCGGCCTTTGAGCGTGTGGTGAATACGCCAACGCgg < 2:419133/90‑1 (MQ=255) gAAATCAAAGATGCGCTCTCGTATCTGCGCCTGATTGCCAACCGCAACGACGACGCGGCCTTTGAGCGTGTGGTGAATACGCCAACGCgg < 1:329195/90‑1 (MQ=255) tCAAAGATGCGCTCTCGTATCTGCGCCTGATTGCCAACCGCAACGACGACGCGGCCTTTGAGCGTgtggtg > 2:334865/1‑71 (MQ=255) tCAAAGATGCGCTCTCGTATCTGCGCCTGATTGCCAACCGCAACGACGACGCGGCCTTTGAGCGTgtggtg < 1:334865/71‑1 (MQ=255) | GTGCTCGAAGAGGCGTTATTGCAGGCCAGTATGCCGTACCGTATTTACGGCGGGATGCGCTTCTTC‑GAACGCCAGGAAATCAAAGATGCGCTCTCGTATCTGCGCCTGATTGCCAACCGCAACGACGACGCGGCCTTTGAGCGTGTGGTGAATACGCCAACGCGG > NZ_CP009273/3992426‑3992590 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |