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.7% | 11.0 / 18.4 | 28 | Y557S (TAC→TCC) | maeA | malate dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (4/7); new base G (17/0); total (21/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.79e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.88e-02 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCT > NZ_CP009273/1548174‑1548335 | gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCgg > 1:29065/1‑90 (MQ=255) gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGGATTCgg > 2:449941/1‑90 (MQ=255) tGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGGATTCGGCTTg > 1:1230991/1‑90 (MQ=255) gagGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGCGGGATTTGGCTTGCCCGAAAatat > 2:171532/1‑90 (MQ=255) gagGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGCGGGATTTGGCTTGCCAGAAAatat > 2:821330/1‑90 (MQ=255) gggAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTTGGCTTGCCCGGAAatattg > 2:1066551/1‑88 (MQ=255) gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGCGGGATTCGGCTTTCCCGAAAATATTGTc > 2:441693/1‑90 (MQ=255) aaaTAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCa > 1:741144/1‑90 (MQ=255) ccGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCt < 1:172954/90‑1 (MQ=255) aGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGc < 1:262198/90‑1 (MQ=255) cTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTAt < 2:679180/51‑1 (MQ=255) cTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGAt > 1:679180/1‑51 (MQ=255) cTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTACCGTCAATGGCCTGTTGCAgg > 1:462940/1‑90 (MQ=255) cTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTc > 1:51586/1‑90 (MQ=255) gggCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGGGTTCGGGTTGCCCGGAAATATCGTCAATGGGCTGTTTCAGGGCGTCCGcgg > 2:967877/1‑88 (MQ=255) gcgcAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAAAGGCCTGTTGCAGGGCTTCGGCagag > 1:580022/1‑90 (MQ=255) tGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAAGGGCCTGTTGCAGGGTTTCGGCAGAGGTTTTCACCGCCa < 2:1171707/90‑1 (MQ=255) ggcggTGCGGCGGTAGTCGCGGTATTCGCTTTGCCAAAAATTTTCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGTTTTTCACCGCCAc < 2:227061/87‑1 (MQ=255) ggAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCAc < 1:449941/90‑1 (MQ=255) gAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCCGAAAATATCGTCAATGGGCTGTTGGAGGGGTTCGGCAGAGGGTTTTACCGccccc > 2:215476/1‑87 (MQ=255) gAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCACg > 1:366785/1‑90 (MQ=255) gAGGTACGGCGGTAGTCGCGGGATTCGGCTTTCCCGAAATTATCGTCAAAGGCCTGTTGCAGGGGTTTGGCAGAGGGTTTTACCGccccc > 2:508846/1‑87 (MQ=255) gAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCCCg > 1:1263410/1‑90 (MQ=255) gAGGCACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACg > 2:1259141/1‑90 (MQ=255) aGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTAAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGc < 1:205004/90‑1 (MQ=255) ggTACGGCGGGAGTCGCGGGATTTGGCTTGCCAGAAAATATCGTCAATGGGCCGTTTGCGGGGTTCGGCAGAGGGTTTCAACCCCACGcc > 2:687447/1‑90 (MQ=255) ggcggTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCACGCCTtgct > 1:78163/1‑90 (MQ=255) ggTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCACGCCTtgctgct > 1:335478/1‑90 (MQ=255) | GATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCT > 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: ‑ |