| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 4,501,865 | 0 | A | C | 56.2% | 6.8 / 16.6 | 16 | L123V (TTG→GTG) | fecD | Fe(3+) dicitrate ABC transporter permease subunit FecD |
| Reads supporting (aligned to +/- strand): ref base A (5/2); new base C (0/9); total (5/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.81e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-02 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGG > NZ_CP009273/4501776‑4501953 | cAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCaa > 1:177948/1‑90 (MQ=255) aTGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATgccgc < 1:124944/90‑1 (MQ=255) gCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGccgccc > 1:614729/1‑90 (MQ=255) aaGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCaaa > 1:251438/1‑90 (MQ=255) gcgggTGAGCGCCAGCTTCATCGGCGGGGGGGTTTTTGCCAGCATCTTCAGAAATATCACCCCCGCCATGCCGCCCGCAAAGGCcagcag < 2:697945/87‑1 (MQ=255) ccAGCTTCATCGGCTGGGGGGTTTTTGCCAGCATCTTCAGAAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCa < 2:614729/90‑1 (MQ=255) cATCGGCGGGGGGGTTTTTGCCAGCATCTTCAGAAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACggg < 2:294608/90‑1 (MQ=255) gCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTg > 1:382169/1‑90 (MQ=255) cTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGa > 1:696415/1‑90 (MQ=255) tggtggGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGAc < 2:696415/90‑1 (MQ=255) gggTCTTTGCCAGCATCTTCAGAAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCa < 2:521467/90‑1 (MQ=255) gtttttGCCAGCATCTTCAGAAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATa < 2:177948/87‑1 (MQ=255) gTCTTTGCCAGCATCTTCAGAAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATa < 2:307329/90‑1 (MQ=255) tGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATaagaag > 1:452038/1‑90 (MQ=255) accATCTCCAGAAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTaga < 1:380391/88‑1 (MQ=255) aGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCcaca < 2:452038/90‑1 (MQ=255) gaaaTACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACag < 2:253828/88‑1 (MQ=255) aaCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAgg > 2:177171/1‑90 (MQ=255) | CAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGG > NZ_CP009273/4501776‑4501953 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |