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.3% | 9.7 / 15.3 | 15 | L123V (TTG→GTG) | fecD | Fe(3+) dicitrate ABC transporter permease subunit FecD |
Reads supporting (aligned to +/- strand): ref base A (6/1); new base C (0/8); total (6/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.40e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGG > NZ_CP009273/4501791‑4501953 | tGCAGAAAGCGCCCCGCCGGTGAGCGCCAGCTTCACCGGCGGGGGGGTCTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATgccgcc < 2:329075/90‑1 (MQ=255) gCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGccgccc > 2:479911/1‑90 (MQ=255) aaaGCGCCACGCGGGTGAGCGCCAGCTTCATCGGCGGGGGGGTTTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCaa < 2:306755/90‑1 (MQ=255) aaGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGGGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCaaa > 2:426484/1‑90 (MQ=255) gAGCGCCAGCTTCATCGGCTGGTGGGGCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGgcag > 1:8245/1‑90 (MQ=255) aGCGCCGGCTTCATGGGCTGGGGGGTTTTGGCCAGCATCTTCAGCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGgcagc < 1:508482/90‑1 (MQ=255) cgcCGGCTCCAGCGGCTGGGGGTTTTTGGCCAGCATCTCAAGCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCAc < 1:322277/90‑1 (MQ=255) aGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCCGCCGCGGCAGCACCATc > 1:177300/1‑90 (MQ=255) aTCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGAAGCACCATCACGGGc < 2:16760/90‑1 (MQ=255) gggTCTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCa < 2:16079/90‑1 (MQ=255) gtttttGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATa < 1:277869/87‑1 (MQ=255) tGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATaagaag > 1:51005/1‑90 (MQ=255) gCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGt > 1:319943/1‑90 (MQ=255) ccAGCATCTTCACCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTa < 2:471219/90‑1 (MQ=255) atCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGcc < 1:426484/90‑1 (MQ=255) aaCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAgg > 1:331033/1‑90 (MQ=255) | TGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGG > NZ_CP009273/4501791‑4501953 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |