Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 349,999 | 0 | T | G | 52.8% | 27.2 / 50.5 | 36 | V528G (GTG→GGG) | prpE | propionate‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base T (2/15); new base G (19/0); total (21/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.44e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-07 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCAC > NZ_CP009273/349916‑350084 | tGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGa > 2:866761/1‑90 (MQ=255) tGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGa > 1:900007/1‑90 (MQ=255) tGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGa > 1:516213/1‑90 (MQ=255) gCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAt < 1:484113/90‑1 (MQ=255) gCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGAAAAAt > 2:503155/1‑90 (MQ=255) cGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGaaa < 2:237679/90‑1 (MQ=255) gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAAGGGCTGGAAGGGGAgggggg > 1:181015/1‑89 (MQ=255) gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGGGCGGAAAGGGCAgggggg > 1:830616/1‑89 (MQ=255) cgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtggcggt < 2:919275/90‑1 (MQ=255) gagagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTc < 1:278100/90‑1 (MQ=255) gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGGTTGTCAt > 2:278100/1‑90 (MQ=255) gTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAAGGGAGGGGGCGGGGGCGGTTGTCAttt > 2:825571/1‑89 (MQ=255) ctcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGaa < 2:287247/90‑1 (MQ=255) ccAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTTGGGGGAAAGAAGGGCTGGAAGGGCAGGGGGCGGTGGGGGTTGGCATTTCGCAag > 2:303660/1‑90 (MQ=255) ccAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAag < 1:833873/90‑1 (MQ=255) ccAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGGGCTTAAAGGGCAGGGGGCGGTGGGGTTTGTTATTGCGAAag > 2:219112/1‑90 (MQ=255) aGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGAAGGGCTGAAAGGGCAGGGGGGGGGGGGGGTTGTCATTCCGGAaggg > 2:981431/1‑88 (MQ=255) aGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGGGCTGGAAGGGGAGGGGGGGGGGGGGGTTGTCATTCCGAAagag > 1:911912/1‑90 (MQ=255) aGCCATCCGGGCGTTGCCGAAGT‑GCGGGGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGGGGGGGCGGTTGTCATTTCGAAagaga > 1:649642/1‑90 (MQ=255) ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc < 2:396403/90‑1 (MQ=255) ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc < 1:164690/90‑1 (MQ=255) ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc < 2:756458/90‑1 (MQ=255) ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAAATTCGCTGAAAGGGGAGGGGGCGGGGGCGTTTGTCATTCCCAAAGAGAGCGAc > 2:640839/1‑90 (MQ=255) tGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGTGGGGGGGGGGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGa > 2:164690/1‑90 (MQ=255) tGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGa > 2:341107/1‑90 (MQ=255) tGCCGAAGTGGCGGGGGGTGGGGGGGAAAATGCGCTGGAAGGGGAGGGGGGGGGGGGGGTTGTTATTCCGAAAGAGAGCGACAGTCTggg > 1:874685/1‑89 (MQ=255) gAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGAt < 1:827520/90‑1 (MQ=255) gAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGGAGGGGGGGGGGGGGTTTGTCATTCCGAAAGAGAGGGACAGTCTGGAAGAt > 1:228399/1‑90 (MQ=255) gAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGGGGTTGTCATTCCGAAAAAGAGGGACAGTCTGGGAGAt > 1:173474/1‑90 (MQ=255) aGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCg < 2:911912/90‑1 (MQ=255) tGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCa < 2:33969/90‑1 (MQ=255) ggggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCAc < 1:113457/90‑1 (MQ=255) ggggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCAc < 2:874685/90‑1 (MQ=255) ggggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCAc > 2:591365/1‑90 (MQ=255) gggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTc > 2:584337/1‑44 (MQ=255) gggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTc < 1:584337/44‑1 (MQ=255) | TGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCAC > NZ_CP009273/349916‑350084 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |