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.1% | 20.5 / 28.2 | 32 | Y557S (TAC→TCC) | maeA | malate dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (5/10); new base G (17/0); total (22/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.65e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.04e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGC > NZ_CP009273/1548172‑1548339 | gagaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTc > 1:821469/1‑90 (MQ=255) agaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTAtttg > 2:891835/1‑88 (MQ=255) gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCgg > 1:1227232/1‑90 (MQ=255) gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCgg > 1:670164/1‑90 (MQ=255) aTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGc < 2:837912/90‑1 (MQ=255) gAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCTGCTTGc > 1:1447394/1‑90 (MQ=255) aGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTTGGCTTTCCAGaa > 1:897198/1‑90 (MQ=255) gagGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAAatat > 1:241378/1‑90 (MQ=255) gagGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGGATTTGGCTTGCCAGAAAatat > 1:888202/1‑90 (MQ=255) aTAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAAt < 2:1227232/90‑1 (MQ=255) aTAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGCGGGATTCGGCTTGCCCGAAATTATCCTCaaa > 2:818839/1‑89 (MQ=255) gCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCCGAAATTTTCGTCAATTGc > 1:326569/1‑90 (MQ=255) gCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGGATCGGGGGGTTTGGGTTGGCCGAAATTATCTTCAATggg > 2:1204672/1‑89 (MQ=255) ccGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCt < 1:1126416/90‑1 (MQ=255) ccGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGACGCGGTATTCGGGTTGGCCGGAATTATCGGCAATTGCCt > 1:653731/1‑90 (MQ=255) ccGGTAGCCTTCACTACCGGGCGCAGGCGTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCt < 2:1295952/90‑1 (MQ=255) cGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTTGGGTTGTCCGAAAATATCTTCAATGGCCTg > 2:731372/1‑90 (MQ=255) cACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCt < 2:1013689/90‑1 (MQ=255) gggCGCAGGCTTAGATGGAGGTACGGCGGGAATCGGGGTATTCGGCTTGCCCGGAAATAAAGTCAAGGGGCTGGTGGAAGGGTTCGGCaa > 1:139814/1‑89 (MQ=255) gcgcAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCag > 1:1547677/1‑88 (MQ=255) gcgcAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCag < 2:1547677/88‑1 (MQ=255) ggAGGTACGGGGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCAc < 1:478250/90‑1 (MQ=255) gAGGTACGGCGGTAGTCGCGGGATTTGGCTTTCCAGAAAATATCGTCAATGGGCTGTTGGAGGGGTTCGGCAGAGGGTTTTACCGCCCCg > 2:1253471/1‑90 (MQ=255) gAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATTGTCAATGGCCCGTTGCAGGGGTTCGGGAGGGGGTTTTACCCCCCCg > 2:478250/1‑90 (MQ=255) gAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTTACCGCCACg > 1:808700/1‑90 (MQ=255) gAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCAcc > 1:578410/1‑89 (MQ=255) aGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGc < 2:241378/90‑1 (MQ=255) cggTAGTCGCGTTATTCGGCTTGCCAGAAATTATCGTCAATGGCGTGTTGCAGGGCTTCGGCAGGGGTTTTCACCGCCACGCCTtgctgc < 1:495898/90‑1 (MQ=255) ggTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCCCCACGCCCtgctgct > 1:576710/1‑90 (MQ=255) ggTAGTCGCGGGATTCGGCTTGCCAGAAAATATCGTCAATGGCCTGTTGGCGGGCTTCGGCAGAGGGTTTCACCGCCCCGCCTtgctgct > 2:85371/1‑90 (MQ=255) aGTCGCGGGATTCGGCTTGCCAGAAATTATCGGCAATGGCCCGTTGGCGGGCTTCGGCAGAGGTTTTCACCCCCACGCCCTGCTGCTgcg > 1:105313/1‑90 (MQ=255) gTCGCGGTATTGGGCTTGCCAAAAATTATCGTCAATGGCCTGTTGCGGGGCTTCGGAAGAGTTTTTCACCGCCACGCCTTGCTGCTgcgc < 2:1552897/90‑1 (MQ=255) | GAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGC > NZ_CP009273/1548172‑1548339 |
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: ‑ |