Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,588,734 | 0 | A | C | 53.3% | 11.3 / 18.0 | 15 | V36G (GTG→GGG) | livM | branched chain amino acid ABC transporter permease LivM |
Reads supporting (aligned to +/- strand): ref base A (3/4); new base C (0/8); total (3/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.69e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAAGCCGGTCGCAAAAGCTGGAAGAAAAAGACCACCGCCGTGCCGATAAACACCCACTGCCAACGGACATCCGAA‑GCCGTGTCGACCACCAGTTTGGTGCCATCCAGCTCCAGTTGCACGCCCATAAAGACGCCCGCCAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGC > NZ_CP009273/3588647‑3588822 | aaaGCGGGTCGAAAAAGCTGGAAAAAAAAGCCCCCCGCCGTGCCGATAAACACCCACTGCCAACGGACATCCGAA‑GCCGTGTCGaccacc < 2:167898/90‑1 (MQ=255) cGGTCGCAAAAGCTGGAAGAAAAAGCCCACCGCCGTGCCGATAAACACCCACTGCCAACGGACATCCGAA‑GCCGTGTCGACCACCAGttt < 1:662294/90‑1 (MQ=255) gCTGGAAAAAAAAGACCCCCGCCGTGCCGATAAACACCCACTGCCAACGGACATCCGAA‑GCCGTGTCGACCACCAGTTTGGTGCCATCCa < 2:511481/90‑1 (MQ=255) gACCACCGCCGTGCCGATAAACACCCACTGCCAACGGACATCCGAA‑GCCGTGTCGACCACCAGTTTGGTGCCATCCAGCTCCAGTTGCAc < 2:110823/90‑1 (MQ=255) ccgccgTGCCGATAAACACCCACTGCCAACGGACATCCGAA‑GCCGTGTCGACCACCAGTTTGGTGCCATCCAGCTCCAGTTGCACGCCCa > 2:25433/1‑90 (MQ=255) ccccccCTGCCACCGGCCACCCGAA‑GCCGTGTCGACCCCCAGTTTGGTGCCATCCAGCTCCAGTTGCACGCCCATAAAGACGCCCGCCAg < 2:167526/88‑1 (MQ=255) cccccACTGCCAACGGACATCCGAA‑GCCGTGTCGACCACCAGTTTGGTGCCATCCAGCTCCAGTTGCACGCCCATAAAGACGCCCGCCAg < 1:25433/88‑1 (MQ=255) cACTGCCAACGGACATCCGAA‑GCCGTGTCGACCACCAGTTTGGTGCCATCCAGCTCCAGTTGCACGCCCATAAAGACGCCCGCCAGcaca > 1:430829/1‑90 (MQ=255) cACTGCCAACGGACATCCGAA‑GCCGTGTCGACCACCAGTTTGGTGCCATCCAGCTCCAGTTGCACGCCCATAAAGACGCCCGCCAGcaca > 1:321382/1‑90 (MQ=255) cAACGGACACCCGAA‑GCCGTGTCGCCCCCCAGTTTGGTCCCATCCAGCTCCAGTTGCACGCCCAAAAAGACGCCCGCCAGCACAAAGAAc < 2:321382/90‑1 (MQ=255) aaCGGACACCCGAACCCCGTGTCCCCCCCCAGTTTGGTCCCACCCAGCTCCAGTTGCCCGCCCAAAAAGACGCCCGCCAGCACAAAGAAc < 2:561324/90‑1 (MQ=255) gACATCCGAA‑GCCGTGTCGCCCCCCAGTTTGGTGCCATCCAGCTCCAGTTGCACGCCCATAAAGACGCCCGCCAGCACAAAGAACATcgc < 1:398014/90‑1 (MQ=255) gACATCCGAA‑GCCGTGTCGACCCCCAGTTTGGGCCCACCCAGCCCCAGTTGCACGCCCAAAAAGACGCCCGCCAGCACAAAGAACATcgc < 2:491212/90‑1 (MQ=255) gACATCCGAA‑CCCGTGTCGACCCCCAGTTTGGTGCCACCCAGCTCCAGTTGCCCGCCCATAAAGACGCCCGCCAGCACAAAGAACATcgc < 2:14094/90‑1 (MQ=255) gACACCCGAA‑GCCGTGTCGCCCCCCAGTTGGGGGCCACCCAGCCCCAGTTGCACGCCCAAAAAGCCGCCCGCCAGCACAAAGAACATcgc < 2:430994/90‑1 (MQ=255) gACCCCCAGTTTGGTGCCATCCAGCTCCAGTTGCACGCCCATAAAGACGCCCGCCAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCAt < 1:116349/90‑1 (MQ=255) caccaGTTTGGTCCCACCCAGCTCCAGTTGCACCCCCAAAAAGCCGCCCGCCAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGc < 2:19459/90‑1 (MQ=255) | AAAGCCGGTCGCAAAAGCTGGAAGAAAAAGACCACCGCCGTGCCGATAAACACCCACTGCCAACGGACATCCGAA‑GCCGTGTCGACCACCAGTTTGGTGCCATCCAGCTCCAGTTGCACGCCCATAAAGACGCCCGCCAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGC > NZ_CP009273/3588647‑3588822 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |