Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,492,959 | 0 | T | G | 55.6% | 18.5 / 22.0 | 18 | G51G (GGT→GGG) | yhfL | YhfL family protein |
Reads supporting (aligned to +/- strand): ref base T (2/6); new base G (10/0); total (12/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.51e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACA‑CCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAGGGCGTATTTGTT > NZ_CP009273/3492870‑3493047 | ggTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGt < 2:569287/90‑1 (MQ=255) ggTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGt < 1:462163/90‑1 (MQ=255) tCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTcc < 1:725837/90‑1 (MQ=255) aaaaCCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTTCGGGCCTACTGcac > 1:81311/1‑90 (MQ=255) aaaaCCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTTCGGGCCCTCTGccc > 2:147637/1‑88 (MQ=255) gTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTGCTACTGCCCAGGaaa > 2:755262/1‑89 (MQ=255) gTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGGCCCACTTCACAGTAAg > 2:224535/1‑90 (MQ=255) aaCTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTTCGGGCCCACTGCCCAAGAAGGCTCTGAc > 1:530278/1‑90 (MQ=255) gCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTTCGGGCCTACTCCCCAGTAAATCTCTTACAGaa > 2:641964/1‑90 (MQ=255) gCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGGCCTACTGCACAGTAAGTCTCTTACAGaa > 2:175066/1‑90 (MQ=255) ctTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATc < 2:19439/90‑1 (MQ=255) tACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGcac < 2:412596/67‑1 (MQ=255) tACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGcac > 1:412596/1‑67 (MQ=255) cTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGGTTCGGGCCCTCCGCCCAATAAATTTCTTACAGAAAACCGCTaaa > 2:264725/1‑89 (MQ=255) tACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGGCCCCCTTCCCAGTAAATCTCTGACAGAAAACGGCTAaca > 2:587157/1‑90 (MQ=255) tgctgcACCCGGCAATTTCTATTTCTAAAATCATTGGCGGGTGCGGGCCTACTTCACAATAAATCTCTTTAAGAAAAAGGCCAACA‑CCgg > 1:424033/1‑90 (MQ=255) aCCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTTCGGGCCCACTCCACAGTAAATCTCTTACAAAAAACGGCTAACACCCGGTTTat > 2:628646/1‑89 (MQ=255) ggTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACA‑CCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAGGGCGTATTTGt < 1:637551/90‑1 (MQ=255) ggTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACA‑CCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAGGGCGTATTTGt < 1:395812/90‑1 (MQ=255) gTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACA‑CCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAGGGCGTATTTGtt < 1:147637/90‑1 (MQ=255) | GGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACA‑CCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAGGGCGTATTTGTT > NZ_CP009273/3492870‑3493047 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |