Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,548,257 | 0 | T | G | 53.3% | 28.3 / 35.8 | 30 | Y557S (TAC→TCC) | maeA | malate dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (4/10); new base G (16/0); total (20/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.33e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTC‑AATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCT > NZ_CP009273/1548169‑1548335 | gATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTa > 2:1310824/1‑90 (MQ=255) gagaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTc < 2:1212209/90‑1 (MQ=255) gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCgg > 2:1043491/1‑90 (MQ=255) gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCgg > 2:824933/1‑90 (MQ=255) gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCgg > 2:1461554/1‑90 (MQ=255) gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGGATTTgg > 1:568915/1‑90 (MQ=255) gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGCGGGATTCgg > 1:1468219/1‑90 (MQ=255) tCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCt < 1:1043491/90‑1 (MQ=255) aaaGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAg < 1:645851/90‑1 (MQ=255) aGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGCGGTTTTTCGCGTGccggga > 1:1022278/1‑86 (MQ=255) agGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATc < 1:791603/90‑1 (MQ=255) agGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGCGGGAATCGGCTTGCCAGAAAatatt > 1:1521110/1‑89 (MQ=255) aTAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTC‑AAt < 1:1437549/90‑1 (MQ=255) ccGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTTGGCTTGCCAGAAAATATCGTC‑AATGGCCt > 2:519949/1‑90 (MQ=255) cGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGGATTCGGCTTGCCAGAAATTATTGTC‑AATgggccg > 2:702875/1‑86 (MQ=255) cTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGCGGGATTCGGCTTGCCCGAAAATATTGTCTAATGGCCGTTTGcgg > 1:1255938/1‑88 (MQ=255) aCTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTC‑AATGGCCTGTTGCAGGGCtt < 1:271605/90‑1 (MQ=255) aCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTTGGCTTGCCCGAAATTATCGTC‑AATGGCCTGTTGCAGGGGTTCgg > 1:1501963/1‑90 (MQ=255) gggCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTC‑AATGGCCTGTTGCAGGGCTTCGGCag < 1:1310824/90‑1 (MQ=255) gggCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTACTCGGCTTGCCAGAAATTATCGTC‑AATGGCCTGTTGCAGGGCTTCGGCag < 2:1521110/90‑1 (MQ=255) gcgcAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTC‑AATGGCCTGTTGCAGGGCTTCGGCagag > 2:31127/1‑90 (MQ=255) gcgcAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAAATATCGTC‑AATGGCCTGTTGCAGGGCTTCGGCagag > 2:511017/1‑90 (MQ=255) gggggTACGGCGGTCGTCGCGGTTTTCGGCTTGCCAAAAATTATCGTC‑AATGGCCTGTTGCGGGGCTTGGGCAGAGTTTTTCACCGCCAc < 1:443956/87‑1 (MQ=255) ggAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAAAAATTATCGTC‑AATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCAc < 2:999080/90‑1 (MQ=255) gAGGTACGGCGGTAGTCGCGGGATTTGGCTTGCCAGAAAATATTGTC‑AATGGGCTGTTGCAGGGGGTCGGGGGGGGGTTTTACCGccccc > 1:999080/1‑87 (MQ=255) gAGGTACGGCGGTAATCGCGGGATTCGGCTTGCCCGAAAATATCGTC‑AAAGGGCCGGTGGCGGGGGTCGGGGGAGGGTTTTACCCCCAcc > 1:1000905/1‑89 (MQ=255) tACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTC‑AATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTTACCGCCACGccct > 2:233143/1‑88 (MQ=255) aCGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTC‑AATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCACGCCTtg > 2:723860/1‑90 (MQ=255) aCGGCGGCAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTC‑AATGGCCCGTTGCACGGCTTCGGCAGAGGGTTTTACCCCCACGCCTtg > 2:443956/1‑90 (MQ=255) ggTAGTCGCGGGATTCGGCTTTCCAGAAATTATTGGC‑AATGGCCCGTTGCAGGGCTTCGGCAGAGGGTTTTACCGCCTCGCCTttctgct > 2:1610074/1‑90 (MQ=255) | GATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTC‑AATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCT > NZ_CP009273/1548169‑1548335 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |