| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,363,696 | 0 | T | C | 71.4% | 21.0 / 13.4 | 21 | F16S (TTC→TCC) | pspD | phage shock protein PspD |
| Reads supporting (aligned to +/- strand): ref base T (2/4); new base C (5/10); total (7/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTTGCT > NZ_CP009273/1363616‑1363781 | tACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGc < 1:28605/90‑1 (MQ=255) tAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGc < 1:187126/90‑1 (MQ=255) tAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGc < 2:159671/90‑1 (MQ=255) ccGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAg < 1:248710/90‑1 (MQ=255) ttCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTgg > 2:649242/1‑90 (MQ=255) tCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTActt < 1:689371/90‑1 (MQ=255) gtgAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTAcc > 2:111955/1‑90 (MQ=255) ggAAAGTTATGAATACTCGGTGGGAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCAc < 2:223760/90‑1 (MQ=255) ggAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCAc < 1:649242/90‑1 (MQ=255) ggAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCAc < 2:815901/90‑1 (MQ=255) aaaGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTg < 2:304995/90‑1 (MQ=255) aaaGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTg < 2:124495/90‑1 (MQ=255) aaGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTgc < 2:581363/90‑1 (MQ=255) aTACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGcc > 2:211266/1‑90 (MQ=255) cTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCcgg > 2:204680/1‑90 (MQ=255) cTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCggg < 2:541571/90‑1 (MQ=255) ggCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGgtg < 1:816368/90‑1 (MQ=255) aGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGGGGCgg > 2:812378/1‑90 (MQ=255) gggCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTActt > 2:235415/1‑48 (MQ=255) gggCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTActt < 1:235415/48‑1 (MQ=255) ggTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGGCGGGGGGTGAAATAAGTTggc > 1:428698/1‑88 (MQ=255) | TACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTTGCT > NZ_CP009273/1363616‑1363781 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |