| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 4,501,865 | 0 | A | C | 54.2% | 13.0 / 24.1 | 24 | L123V (TTG→GTG) | fecD | Fe(3+) dicitrate ABC transporter permease subunit FecD |
| Reads supporting (aligned to +/- strand): ref base A (8/3); new base C (0/13); total (8/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.24e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGG > NZ_CP009273/4501782‑4501953 | ggCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGGGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGc > 2:37058/1‑90 (MQ=255) gCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGGGGGGCTTTTCCAGCATCTTCAGCAATATCAACCCCGCCATgcc > 2:238529/1‑90 (MQ=255) cATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATgccg < 1:136318/90‑1 (MQ=255) gCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTTCCAGCATCTTCAGCAATATCAACCCCGCCATGccgccc > 2:459967/1‑90 (MQ=255) aaaGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCaa < 1:108181/90‑1 (MQ=255) cgcCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGc > 2:106612/1‑90 (MQ=255) cGCCGGGGGGCGCCGGCTTCACCGGCGGGTGGGTCTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCcagca < 2:180331/90‑1 (MQ=255) ggTGAGCGCCAGCTTCATCGGCTGGGGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCgg > 2:187199/1‑90 (MQ=255) gggAGCGCCAGCTTCATCGGTGGGTGGGTTTTTGCCAGCATCTTCAGAAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGgc < 2:396978/88‑1 (MQ=255) gTGAGCGCCAGCTTCATCGGCTGGTGGGTTTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGgc < 1:106612/90‑1 (MQ=255) gCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGGCAGCAGCGGCAGCACCATCa > 2:212774/1‑90 (MQ=255) cGGCTGGGGGGTTTTTGCCAGCATCTTCAGCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAg < 2:186701/90‑1 (MQ=255) gtggGTCTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACgg < 1:37058/90‑1 (MQ=255) gggTTTTTGCCAGCATCTTCAGAAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCa < 1:238529/90‑1 (MQ=255) gtttttGCCAGCATCTCCAGAAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATa < 2:54458/87‑1 (MQ=255) ccAGCATTTTCAGAAATCTCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTa < 2:138047/90‑1 (MQ=255) aGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTaga < 2:105647/90‑1 (MQ=255) aGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTaga < 1:143624/90‑1 (MQ=255) aGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTaga < 1:384287/90‑1 (MQ=255) aGCATCTTCAGAAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTaga < 2:56139/90‑1 (MQ=255) cAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCcac > 2:387618/1‑90 (MQ=255) gCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACag < 2:287258/90‑1 (MQ=255) gCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACag < 2:138835/90‑1 (MQ=255) atatCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAgg > 2:125236/1‑90 (MQ=255) aaCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAgg > 1:192155/1‑90 (MQ=255) aaCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACACAGGCCAgg > 2:456120/1‑90 (MQ=255) | GGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGG > NZ_CP009273/4501782‑4501953 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |