Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,548,257 | 0 | T | G | 60.0% | 25.9 / 34.5 | 35 | Y557S (TAC→TCC) | maeA | malate dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (4/10); new base G (21/0); total (25/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.45e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.13e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCC‑TGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCT > NZ_CP009273/1548174‑1548335 | gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGAGGCGGGATTCgg > 1:106520/1‑90 (MQ=255) gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGCATCGCGGGAttatg > 1:676619/1‑87 (MQ=255) gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCCGGATTCgg > 1:472780/1‑90 (MQ=255) gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGCGGGATTTgg > 1:197221/1‑90 (MQ=255) aTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGc < 2:142025/90‑1 (MQ=255) gAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGGCCGGGATTCGGCGTGc > 2:680887/1‑90 (MQ=255) agGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATc < 1:1142955/90‑1 (MQ=255) gCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGGATTCGGGTTGCCAGAAAATATTGTCAAAGGc > 1:589442/1‑90 (MQ=255) gCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGCGGGATTTGGGTTGGCCGAAAATATTGTCAAAGGc > 1:620078/1‑90 (MQ=255) cccGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGcc < 1:1105235/90‑1 (MQ=255) cGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTTCCAGAAATTATCGTCAATGGcccg > 2:784880/1‑88 (MQ=255) gCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCC‑TGTTGCa < 1:50662/90‑1 (MQ=255) gCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCC‑TGTTGCa < 2:669902/90‑1 (MQ=255) cTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTGGTCGCGGTATTCGGGTGGCCCGAAAATATTGGCAATGGGC‑TGGTggggg > 1:1077508/1‑86 (MQ=255) aCCGGGCGCAGGCTTAGATAGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCC‑CGTTGCAGGGCTTCgg > 2:440683/1‑90 (MQ=255) cGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTTGGCTTTCCCGAAAATATTGGCAATGGGC‑CGTTGCCGGGCTTCGGcc > 1:1103394/1‑89 (MQ=255) gggCGCAGGCTTAGATGGAGGTACGGGGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCC‑TGTTGCAGGGCTTCGGCag > 1:102756/1‑90 (MQ=255) gggCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCC‑TGTTGCAGGGCTTCGGCag < 1:891461/90‑1 (MQ=255) gggCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCC‑TGTTGCAGGGCTTCGGCag < 1:940369/90‑1 (MQ=255) gggCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCC‑TGTTGCAGGGCTTCGGCag > 2:77955/1‑90 (MQ=255) gggCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATAGTCAATGGCC‑TGTTGTAGGGGTTCGGCag > 1:379890/1‑90 (MQ=255) gggCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAAATATCGTCAATGGCC‑TGTTGCAGGGCTTCGGCag > 2:729265/1‑90 (MQ=255) gcgcAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTTGGCTTTCCAGAAATTATCGTCAATGGCC‑TGTTGCAGGGCTTCGGCagag > 2:714115/1‑90 (MQ=255) gcgcAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTTCCAGAAATTATCGTCAATGGCC‑TGTTGCAGGGCTTCGGCagag > 2:905891/1‑90 (MQ=255) cAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGGTTTCCCGAAAATATCGTCAATGGGC‑CGTTGCCGGGGTTCGGGAGAgggt > 1:912131/1‑88 (MQ=255) ggCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCCTTCCAGAAATTATCGTCAATGGCC‑TGTTGCAGGGGTTCGGCCGAGGgttt > 1:48513/1‑90 (MQ=255) ggAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCC‑TGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCAc < 2:396329/90‑1 (MQ=255) gAGGTACGGCGGTAGTCGCGGGATTCGGCTTGGCAGAAATTATCGTCAATGGGC‑TGTTGAAGGGCTTCGGGGGAGGGTTTCACCGCCCCg > 2:864299/1‑90 (MQ=255) gAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCC‑TGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCCCg > 2:1082095/1‑90 (MQ=255) aCGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAAATATCGTCAATGGCA‑TGTTGCAGGGGTTCGGGAGAGGGTTTTACCGCCACCCCCtg > 2:572655/1‑90 (MQ=255) cggTGGTCGCGGTATACGAGTTACCATAAGTTATC‑TTATTGGCCGTTTGGCGGGGTTTGGGAAGAGTTTTCCACCGCCACGCCTtgctgc < 1:992437/90‑1 (MQ=255) ggTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCaa > 2:147070/1‑39 (MQ=255) ggTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCaa < 1:147070/39‑1 (MQ=255) ggTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCC‑TGTTGCAGGGCTTCGGCAGAGGGTTTTACCGCCACGCCCtgctgct > 1:277417/1‑90 (MQ=255) ggTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCC‑TGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCCCGCCCtgctgct > 2:736859/1‑90 (MQ=255) | GATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCC‑TGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCT > NZ_CP009273/1548174‑1548335 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |