| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 3,594,920 | 0 | A | C | 53.8% | 12.8 / 16.1 | 13 | V176G (GTG→GGG) | ftsX | permease‑like cell division protein FtsX |
| Reads supporting (aligned to +/- strand): ref base A (6/0); new base C (0/7); total (6/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.83e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCA > NZ_CP009273/3594835‑3594997 | cACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATcaccac > 2:447089/1‑90 (MQ=255) ctttATTCTGCGTGATACGATCCCCCTGCTTTTCCAGTGTTCCCGCCCCTTGAAACTCGAGTTGGGGAACCCCCACCGCCACTGCCGGaa < 1:371462/88‑1 (MQ=255) ttAATCTGCGTGATACGACCCCGCGGCTTATTAAGTGTTTCCGCCCCCTGAAAATCGATTTTGGGGATCCCCACCGCCACTGCCGGAAGc < 2:273110/90‑1 (MQ=255) gCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGtttt > 1:41266/1‑90 (MQ=255) gCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGtttt > 2:227275/1‑90 (MQ=255) aTCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCat > 1:98197/1‑90 (MQ=255) ccAGCGTTTCCAGTGTTCCCGCCCCCTGAAAATCGAGTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCa < 2:41266/89‑1 (MQ=255) ttttCAGTGCTCCCGTCCCGTGAAATTCGGTTTTCGGCACCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGcac < 1:227275/88‑1 (MQ=255) tATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGcac < 2:199353/90‑1 (MQ=255) tATTCAGTGATTCCGCCCCCTGAAAATCGAGTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGcac < 2:217033/90‑1 (MQ=255) atttCAGTGATTCCGCCCCCTGAAAATCGATTTTGGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGcac < 2:313071/88‑1 (MQ=255) tCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAacca > 1:326778/1‑90 (MQ=255) ccGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAaccct > 1:313071/1‑88 (MQ=255) aGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCa > 2:247005/1‑90 (MQ=255) | CACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCA > NZ_CP009273/3594835‑3594997 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |