Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,548,257 | 0 | T | G | 51.1% | 45.9 / 55.0 | 45 | Y557S (TAC→TCC) | maeA | malate dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (11/11); new base G (23/0); total (34/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.95e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCAT > NZ_CP009273/1548169‑1548342 | gATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTa > 1:65413/1‑90 (MQ=255) gagaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGGATTc > 1:65007/1‑90 (MQ=255) gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCgg > 1:717985/1‑90 (MQ=255) gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTTgg > 2:1769953/1‑90 (MQ=255) gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCgg > 1:383291/1‑90 (MQ=255) gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTTgg > 2:1467567/1‑90 (MQ=255) gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTTgg > 1:780782/1‑90 (MQ=255) gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCCGGATTTgg > 2:1031199/1‑90 (MQ=255) gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGGATTTgg > 2:217034/1‑90 (MQ=255) gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGGATTCgg > 2:1192319/1‑90 (MQ=255) gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGGACGGCGGTGGGCGCGGGATTccg > 1:1012903/1‑88 (MQ=255) gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGAAGCCTTCACTACCGGGCGCAGGCTTGGATGGAGGTACGGCGGAAGACGGGGGAtatgg > 2:1293652/1‑86 (MQ=255) tGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTg > 1:322999/1‑90 (MQ=255) aGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGGATTTGGCTTGCCAGaa > 1:1033818/1‑90 (MQ=255) ggAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAAt < 2:717985/90‑1 (MQ=255) ggAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAAt < 2:1308667/90‑1 (MQ=255) agGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATc < 1:85396/90‑1 (MQ=255) gggAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCg < 2:1320345/90‑1 (MQ=255) cccGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGcc < 1:929588/90‑1 (MQ=255) ccGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCCGAAAATATCGTCAATGGCCt > 2:908988/1‑90 (MQ=255) acagtCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCa < 1:710874/86‑1 (MQ=255) ccGGGCGCAGGCTTAGATGGAGGTACGGCGGCGGTCGCGGGAATCGGCTTGCCCCCAAATATCGGCCATGGGCCGTTGCCGGGGTTCGGc > 2:1500651/1‑90 (MQ=255) gggCGCAGGCTTAGATGGAGGTACGGCGGTTGACGCGGTATTCGGCTTTGCCGAAAATATCGTTAACGGGCTGTTGCAGGGCTTCTGcgg > 2:1199846/1‑88 (MQ=255) gggCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCag < 2:65413/90‑1 (MQ=255) gggCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGCGGGATTCGGCTTGCCCGAAATTATTGTCAAAGGGCTGTTTCAGGGGGTCgggcg > 2:832156/1‑87 (MQ=255) ggCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCaga < 1:1881678/90‑1 (MQ=255) gcgcAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGGCTGTTGGAGGGGTTCGGcctgt > 2:710874/1‑86 (MQ=255) gcgcAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCagag > 1:1451497/1‑90 (MQ=255) cAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCCTGCCCGAAAATATCGGCAAAGGGCCGTTTGAGGGGTTCGGGAGAgggt > 2:894196/1‑88 (MQ=255) tGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCa < 1:537133/90‑1 (MQ=255) gAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCACg > 1:244746/1‑90 (MQ=255) gAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCACg > 1:691963/1‑90 (MQ=255) gAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTTACCGCCACg > 1:1583490/1‑90 (MQ=255) gAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAAATATTGTCAATGGGCTGTTGGAGGGGTTCGGCAGAGGGTTTTACCGCCCCg > 2:1379701/1‑90 (MQ=255) gAGGTACGGCGGTAGGCGCGGTATTTGGCTTTCCAGAAATTATTGTCAATGGCCCGTTGCAGGGGTTCGGGCGAGGGTTTTACCCCCACg > 2:900656/1‑90 (MQ=255) gAGGTACGGCGGTAGGCGCGGGATTTGGCTTGCCAGAAAATATTGTCAAAGGGCCGTTGGAGGGGGTCGGGAGAGGGTTTTACCCCCCCg > 2:919348/1‑90 (MQ=255) gAGGTACGGCGGGAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAAAGGCCGGTTGCAGGGGTTCGGCAGAGGGTTTCACCGCCACg > 2:1367408/1‑90 (MQ=255) ggTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGcc > 1:316267/1‑90 (MQ=255) tACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCACGCCtt > 1:755238/1‑90 (MQ=255) aCGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCag < 2:696696/68‑1 (MQ=255) aCGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCag > 1:696696/1‑68 (MQ=255) cggcggTAGTCGCGGGATTTGGCTTGCCAGAAAATATCGTCAATGGCCTGTTGCAGGGGTTCGGCAGAGGGTTTTACCCCCCCGcccttc > 2:1846702/1‑86 (MQ=255) ggTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAAt > 1:1406699/1‑40 (MQ=255) ggTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAAt < 2:1406699/40‑1 (MQ=255) gcgGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCAt > 1:695620/1‑90 (MQ=255) | GATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCAT > NZ_CP009273/1548169‑1548342 |
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: ‑ |