| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 108,789 | 0 | T | C | 63.2% | 2.2 / 13.6 | 19 | S100G (AGC→GGC) | zapD | cell division protein ZapD |
| Reads supporting (aligned to +/- strand): ref base T (5/2); new base C (0/12); total (5/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.81e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-03 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATT‑AATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATCCAGGTCTGGAGTTTACGTTGC > NZ_CP009273/108706‑108877 | gATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATT‑AATACGCTCCccgc < 2:332291/90‑1 (MQ=255) gCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATT‑AATACGCTCCCcgccgc > 1:801117/1‑90 (MQ=255) gACCCACCAGAGAAATCAAACGATCTTCACGAAAAAATTGCCCGATACGCGGCGCGAAAATT‑AATCCGCCCCCCGCCGCTTTTAActgct < 2:844259/90‑1 (MQ=255) aacaaaaCGATCTTCACGCAAAAATTGCCCGATCCGCGGGGCGGAAATT‑AATCCCCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTc < 2:437812/86‑1 (MQ=255) aaCGATCTTCCCGAAAAATTTCCCGGATACGCGGCGCGAAAATT‑AACCCCCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATAc < 1:193367/90‑1 (MQ=255) aaCGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATT‑AATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATAc < 2:459139/90‑1 (MQ=255) tCTTCACGCAAAAATTGCCCGATACGCGGCGCGGAAATT‑AATACGCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGctc < 2:98032/90‑1 (MQ=255) cTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATT‑AATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGctct > 1:221752/1‑90 (MQ=255) cGCAAAATTTGCCCGATCCGCGGCGGGAAAATT‑AATCCGCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTcc < 1:189513/90‑1 (MQ=255) aaaaTTGCCCGATACGCGGCGCGAAAATT‑AATACGCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGc < 2:509388/89‑1 (MQ=255) tGCCCAATCCGCGGCGCGAAAATT‑AATCCGCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGc < 2:792964/90‑1 (MQ=255) gATCCGCGGCGCGGAAATT‑AATACGCCCCCCGCCGCTTTTAACTGCTGAATTAAt < 2:134350/55‑1 (MQ=255) gATACGCGGCGCGGAAATT‑AATACGCTCCCCGCCGCTTTTAACTGCTGAATTAAt > 1:134350/1‑55 (MQ=255) ggggcgggaAAATT‑AATCCCCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATcc < 1:553125/88‑1 (MQ=255) gcggcg‑ggaAATTAAATACGCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATcc < 1:722460/90‑1 (MQ=255) cgggcgcggAAAAT‑ATGTCGCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATcc < 1:387508/88‑1 (MQ=255) gcgaAAATT‑AATCCGCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATCCAGGTc < 2:142584/90‑1 (MQ=255) cGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATCCAGGTCTGGAGTTTACGtt > 1:371829/1‑90 (MQ=255) cTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATCCAGGTCTGGAGTTTACGTtgc > 1:583241/1‑90 (MQ=255) | GATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATT‑AATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATCCAGGTCTGGAGTTTACGTTGC > NZ_CP009273/108706‑108877 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |