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% | 21.4 / 37.6 | 32 | L123V (TTG→GTG) | fecD | Fe(3+) dicitrate ABC transporter permease subunit FecD |
Reads supporting (aligned to +/- strand): ref base A (11/3); new base C (0/18); total (11/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.82e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCC > NZ_CP009273/4501776‑4501950 | cAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGGGGGTCTTTGCCAGCATCTTCAGCAATATCaa > 1:163213/1‑90 (MQ=255) gCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGGGGGGCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCCTGccgccc > 1:138937/1‑90 (MQ=255) aGCGCCAGCTCCATCGGCTGGGGGGTCTTGGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGgcagc < 1:116859/90‑1 (MQ=255) cgcCAGCTTCATCGGCTGGTGGGTTTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCAc < 2:163213/90‑1 (MQ=255) cgcCAGCTTCACCGGCTGGGGGGTCTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCAc < 1:375544/90‑1 (MQ=255) cgcCAGCTTCACCGGCGGGGGGGTTTTGGCCAGCATCTTCAGCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGAAGCAc < 1:960609/90‑1 (MQ=255) aGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATc > 1:421922/1‑90 (MQ=255) cTTCATGGGCTGGTGGGTCTTTGCCACCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCAc < 2:63368/90‑1 (MQ=255) cTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCAc < 2:39578/90‑1 (MQ=255) cTTCATCGGCTGGGGGGTTTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCAc < 2:573575/90‑1 (MQ=255) cATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACggg > 2:789586/1‑90 (MQ=255) gCTGGTGGTTTTTTGCCAGCATCTTCAGAAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTg < 1:424613/90‑1 (MQ=255) gCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGc > 1:960439/1‑61 (MQ=255) gCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTg < 2:421922/90‑1 (MQ=255) gCTGGGGGGTTTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGc < 2:960439/61‑1 (MQ=255) tggtggGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGAc > 2:16165/1‑90 (MQ=255) gtggGTTTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCAc < 2:916723/72‑1 (MQ=255) gtggGTTTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACgg < 2:1072345/90‑1 (MQ=255) gtggGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCAc > 1:916723/1‑72 (MQ=255) gtttttGCCAGCATCTTCAGAAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATa < 1:335206/87‑1 (MQ=255) tttGCCAGCATTTCCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATaaga < 2:838950/90‑1 (MQ=255) ccAGCATCTCCAGAAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACTATCACGGGCAGTGACGGCATAAGAAGTa < 1:581850/90‑1 (MQ=255) aGCATCTTCAGAAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTaga < 1:955254/90‑1 (MQ=255) gCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTagag > 2:989367/1‑90 (MQ=255) cTTCAGAAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGcccc < 1:192484/90‑1 (MQ=255) gaaaTATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATa < 1:924505/69‑1 (MQ=255) gaaaTACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACag < 1:747529/88‑1 (MQ=255) gaaaTACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACag < 1:434116/88‑1 (MQ=255) gCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACag < 2:781826/90‑1 (MQ=255) gCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATa > 2:924505/1‑71 (MQ=255) atatCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAgg > 2:521286/1‑90 (MQ=255) atCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGcc > 2:548191/1‑90 (MQ=255) atCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGcc > 1:768972/1‑90 (MQ=255) | CAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCC > NZ_CP009273/4501776‑4501950 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |