Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 328,702 | 0 | G | T | 35.3% | 68.0 / 12.2 | 34 | L207F (TTG→TTT) | pdeL | DNA‑binding transcriptional activator/c‑di‑GMP phosphodiesterase PdeL |
Reads supporting (aligned to +/- strand): ref base G (3/19); new base T (12/0); total (15/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.30e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGT > NZ_CP009273/328617‑328789 | aGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTGGCAg > 1:166428/1‑90 (MQ=255) gACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGGAGc > 2:326486/1‑90 (MQ=255) aCTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCg > 1:6472/1‑90 (MQ=255) cTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGGAGCgg > 1:250701/1‑90 (MQ=255) ggATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGAtt < 1:125533/90‑1 (MQ=255) tGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATAATGGCATCAACGTCTCGGCGGGGTGTTTTTTTGGAGCGGGGTTTTtaaaaa > 1:445747/1‑89 (MQ=255) gATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGGTTTTaaaaaaa > 2:454814/1‑88 (MQ=255) aTGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAagag < 2:234683/90‑1 (MQ=255) aTGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAagag < 2:124536/90‑1 (MQ=255) cGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGATTTGGAAAaaaggggc > 2:461688/1‑86 (MQ=255) ggTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGGGGGATTTGGAAAAGAGGGTCt > 2:470827/1‑90 (MQ=255) aaCATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAAcc < 2:59971/90‑1 (MQ=255) aaCATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAAcc < 1:203549/90‑1 (MQ=255) aCATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCt < 2:460871/90‑1 (MQ=255) cATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTg < 2:168773/90‑1 (MQ=255) tgctgcCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTta < 2:382600/90‑1 (MQ=255) ccGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGGTTTTAAAAAAAGGGTCTGAAACTGGGTAATaaa > 2:38651/1‑90 (MQ=255) tttCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTaa < 2:119293/90‑1 (MQ=255) atTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGGTTTTAAAAAAAGGGTCTGAAACTGGGTAATAAATTAGGGAAAGATaaa > 2:46813/1‑90 (MQ=255) ggCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGGTTTGAAAAAGAGGGGCTGAACCTGGTTAATAAATTTGGGAACGATaaaaac > 1:303452/1‑88 (MQ=255) gCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCa < 2:290463/90‑1 (MQ=255) cAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCt < 2:363407/90‑1 (MQ=255) tctcGGCGGGTTGTTTTTTTGCAGCGGGGTTTGAAAAAGAGTGTCTGAACCTGGGTAATAAATTAGGGAACGATAAAATCAAGCTGGTtc > 1:341493/1‑90 (MQ=255) tctcGGCGGGGTGTTTTTTTGCAGCGGGGTTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGGAAAGATAAAATCAAGCTgggtt > 2:207379/1‑87 (MQ=255) tcGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTtctc < 1:207379/90‑1 (MQ=255) tcGGCGGGGTGTTTTTTTGCAGCGGGGTTTGAAAAAAAGGGTCTGAAACTGGTTAATAAATTAGGGAAAGATAAAATCAAAGTGGTtttc > 2:340364/1‑90 (MQ=255) ggcggGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGa < 2:6472/90‑1 (MQ=255) cggGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGc < 2:195477/90‑1 (MQ=255) gggTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTTGGGAACGATAAAATCAAGCTGGTTCTCGAGCt > 1:351018/1‑90 (MQ=255) ttGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGt < 2:407843/90‑1 (MQ=255) ttGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGt < 2:18004/90‑1 (MQ=255) ttGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGt < 2:175051/90‑1 (MQ=255) ttGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGt < 1:46813/90‑1 (MQ=255) ttGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGt < 1:110023/90‑1 (MQ=255) | AGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGT > NZ_CP009273/328617‑328789 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |