Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,447,744 | 0 | A | C | 52.9% | 13.2 / 19.7 | 17 | G52G (GGT→GGG) | paaZ | phenylacetic acid degradation bifunctional protein PaaZ |
Reads supporting (aligned to +/- strand): ref base A (8/0); new base C (0/9); total (8/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.11e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACGCTCTTTTTCACTCAGCAGATGTTTAGCGACCGCTTTAAGCATCGCCGCACGTTCGATAAAGGTCATAGCGCGAAGGGCGGGGGCACCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCT > NZ_CP009273/1447657‑1447821 | aCGCTCTTTTTCACTCAGCAGATGTTTAGCGACCGCTTTAAGCATCGCCGCACGTTCGATAAAGGTCATAGCGCGAAGGGCGGGGGCAcc > 1:832210/1‑90 (MQ=255) gCTCTTTTTCACTCAGCAGATGTTTAGCGACCGCTTTAAGCATCGCCGCACGTTCGATAAAGGTCATAGCGCGAAGGGCGGGGGCACCtt > 1:580332/1‑90 (MQ=255) gCGACCGCTTTAAGCATCGCCGCACGTTCGATAAAGGTCATAGCGCGAAGGGCGGGGGCACCTTTTTCAATGGCAAACTGGCGGGCAgcc > 2:245311/1‑90 (MQ=255) cccGTTTAAAGCACCGCCGCCCGTTCGAAAAAGGTCATAGCGCGAGGGGGGGGGGCCCTTTTTTAAAGGGAAACCGGGGGGGCAgccgcc < 1:169753/89‑1 (MQ=255) cccGCTTTAAGCTCCGCCGCACGTTCGAAAAAGGTCATAGCGCGAAGGGGGGGGGCCCTTTTTTAAATGCAAAACTGGGGGGCAgccgcc < 1:17645/89‑1 (MQ=255) aCCGCTTTAAGCATCGCCGCACGTTCGATAAAGGTCATAGCGCGAAGGGCGGGGGCACCTTTTTCAATGGCAAACTGGCGGGCAgccgcc > 1:910338/1‑90 (MQ=255) ccGCTTTAAGCATCGCCGCACGTTCGATAAAGGTCATAGCGCGAAGGGCGGGGGCACCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCa > 1:999700/1‑90 (MQ=255) ttAAGCATCGCCGCACGTTCGAAAAAGGTCATAGCGCGAAGGGCGGGGGCCCCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCa < 2:908058/90‑1 (MQ=255) tgaaGCCTCGCCGCCCGTCAAAAAAAGGTCATAGGGCGAAGGGCGGGGGCCTTTTTTACAAGGGCAAACGGGCGGGCCGCCGCCATATCa < 1:859527/88‑1 (MQ=255) cATCGCCGCACGTTCGATAAAGGTCATAGCGCGAAGGGCGGGGGCCCCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGAcc < 2:372578/90‑1 (MQ=255) tCGCCGCACGTTCGATAAAGGTCATAGCGCGAAGGGCGGGGGCACCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGACCtt > 1:582573/1‑90 (MQ=255) cgccgcACGTTCGATAAAGGTCATAGCGCGAAGGGCGGGGGCACCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGACCTTc > 2:450096/1‑90 (MQ=255) ccgcCCGTTCGAAAAAGGTAATAGCGGGAGGGGGGGGGGCCCTTTTTTCAATGGAAAACGGGGGGGCAGCCGCCATATCAAGACCTTCAc < 2:999700/90‑1 (MQ=255) cATAGCGCGAAGGGCGGGGGCACCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGc > 2:406850/1‑90 (MQ=255) tAGCGCGAGGGGCGGGGGCCCCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCt < 2:148537/90‑1 (MQ=255) gcgcGAGGGGGGGGGCCCCCTTTTTAAATGGAAAACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTcg < 2:954510/90‑1 (MQ=255) aaGGGCGGGGGCCCCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCt < 2:353987/90‑1 (MQ=255) | ACGCTCTTTTTCACTCAGCAGATGTTTAGCGACCGCTTTAAGCATCGCCGCACGTTCGATAAAGGTCATAGCGCGAAGGGCGGGGGCACCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCT > NZ_CP009273/1447657‑1447821 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |