| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 4,501,865 | 0 | A | C | 53.8% | 10.3 / 15.5 | 13 | L123V (TTG→GTG) | fecD | Fe(3+) dicitrate ABC transporter permease subunit FecD |
| Reads supporting (aligned to +/- strand): ref base A (5/1); new base C (0/7); total (5/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.66e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAG > NZ_CP009273/4501776‑4501945 | cAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGGCTTTGCCAGCATCTTCAGCAATATCaa > 1:71637/1‑90 (MQ=255) gCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCGAAATCAACCCCGCCATgcc > 1:570127/1‑90 (MQ=255) aaaGCCCCCCGCGGGTGAGCGCCAGCTTCATCGGCGGGTGGGTCTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCaa < 1:33899/90‑1 (MQ=255) aaGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCaaa > 2:603686/1‑90 (MQ=255) ccAGCTTCATCGGCTGGTGGGTCTTTGCCACCATCTTCAGCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCa < 1:17813/90‑1 (MQ=255) cTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGa > 2:568850/1‑90 (MQ=255) gggtggGTTTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGAc < 1:414355/89‑1 (MQ=255) gggggggCCTTTGCAGGCATCCTCAAAAAAACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGAc < 1:494566/86‑1 (MQ=255) cTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATaag > 2:580290/1‑90 (MQ=255) ccAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTa < 1:580290/90‑1 (MQ=255) aGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTaga < 1:373309/90‑1 (MQ=255) gCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTagag > 1:361841/1‑90 (MQ=255) gaaaTATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACag < 1:340085/88‑1 (MQ=255) gaaaTACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACag < 1:219964/88‑1 (MQ=255) | CAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAG > NZ_CP009273/4501776‑4501945 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |