Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,216,130 | 0 | G | T | 48.5% | 31.0 / 11.6 | 33 | C156F (TGT→TTT) | dld | D‑lactate dehydrogenase |
Reads supporting (aligned to +/- strand): ref base G (7/10); new base T (16/0); total (23/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.97e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTGGAAAAAGCCCTCAAACCGCTGGGACGCGAACCGCACTCAGTGATTGGATCATCGTGTATAGGCGCATCGGTCATCGGCGGTATTTGTAACAACTCCGGCGGCTCGCTGGTGCAACGTGGCCCGGCGTAT‑ACCGAAATGTCGTTATTCGCGCGTATAAATGAAGACGGCAAA > NZ_CP009273/2216042‑2216215 | cTGGAAAAAGCCCTCAAACCGCTGGGACGCGAACCGCACTCAGTGATTGGATCATCGTGTATAGGCGCATCGGTCATCGGCGGTATTTGt > 1:1396288/1‑90 (MQ=255) gCTGGGACGCGAACCGCACTCAGTGATTGGATCATCGTGTATAGGCGCATCGGTCATCGGCGGTATTTGTAACAACTCCGGCGGCTCGCt < 1:1040550/90‑1 (MQ=255) gcgAACCGCACTCAGTGATTGGATCATCGTGTATAGGCGCATCGGTCATCGGCGGTATTTTTAAAAACTCCCGCGGCTCGCTGGGGGAAc > 2:655257/1‑90 (MQ=255) gcgAACCGCACTCAGTGATTGGATCATCGTGTATAGGCGCATCGGTCATCGGCGGTATTTTTAAAAACACCCGCGGCCCGCCGGGGCAAc > 2:303553/1‑90 (MQ=255) gAACCGCACTCAGTGATTGGATCATCGTGTATAGGCGCATCGGTCATCGGCGGTATTTGTAACAACTCCGGCGGCTCGCTGGTGCAACGt < 1:958609/90‑1 (MQ=255) gCACTCAGTGATTGGATCATCGTGTATAGGCGCATCGGTCATCGGCGGTATTTTTAAAAACTCCGGCGGCTCGCTGGGGCAACGGGGgcc > 1:403900/1‑90 (MQ=255) gCACTCAGTGATTGGATCATCGTGTATAGGCGCATCGGTCATCGGCGGTATTTGTAACAACTCCCGCGGCTCGCTGGTGGAACGTGGgcc > 1:1086580/1‑90 (MQ=255) gCACTCAGTGATTGGATCATCGTGTATAGGCGCATCGGTCATCGGCGGGATTTTTAAAAACTCCGGCGGCTCGCTGGGGGAACGGGGgcc > 1:1007604/1‑90 (MQ=255) cAGTGATTGGATCATCGTGTATAGGCGCATCGGTCATCGGCGGTATTTTTAACAACTCCGGCGGCTCGCTGGTGCAACGTGGCCCcgcgg > 1:674326/1‑89 (MQ=255) cAGTGATTGGATCATCGTGTATAGGCGCATCGGTCATCGGCGGTATTTTTAACAACTCCGGCGGCTCGCTGGGGGAACGGGGCCCcgcgg > 1:145694/1‑89 (MQ=255) cAGTGATTGGATCATCGTGTATAGGCGCATCGGTCATCGGCGGTATTTTTAACAAATCCGGCGGCTCGCTGGTGGAAAGTGGCCCcgcgg > 1:479505/1‑89 (MQ=255) aGTGATTGGATCATCGTGTATAGGCGCATCGGTCATCGGCGGTATTTGTAAAAAATCCCGCGGCTCCCTGGGGGAACGTGGCCCGGCgga > 1:1293469/1‑88 (MQ=255) ttGGATCATCGTGTATAGGCGCATCGGTCATCGGCGGGATTTTTAAAAACTCCGGCGGCTCGCTGGGGGAACGGGGCCCCGCGGAT‑ACCg > 2:1156578/1‑90 (MQ=255) ggATCATCGTGTATAGGCGCATCGGTCATCGGCGGTATTTTTAACAACTCCCGCGGCTCGCTGGGGCAACGTGGCCCCGCGGAT‑ACCGaa > 1:1283539/1‑90 (MQ=255) ggATCATCGTGTATAGGCGCATCGGTCATCGGCGGTATTTGTAACAACTCCGGCGGCTCGCTGGGGCAACGTGGCCCCGCGGAT‑ACCgga > 1:333667/1‑88 (MQ=255) catcGTGTATAGGCGCATCGGTCATCGGCGGTATTTTTAAAAACTCCCGCGGCTCCGGGGGGGAACGGGGCCCCGGGGATAACCGaaagt > 2:785308/1‑88 (MQ=255) gtATAGGCGCATCGGTCATCGGCGGGATTTGTAACAACCCCGGCGGCTCGCTGGGGGAACGGGGCCCGGCGTAT‑ACCGAAATGTCGTTAt > 1:817163/1‑90 (MQ=255) ggCGCATCGGTCATCGGCGGTATTTTTAACAACCCCCGCGGCTCGCTGGTGCAACGTGGCCCCGCGTAT‑ACCGAAATGTCGTTTTTtgcg > 1:1011243/1‑90 (MQ=255) gcgcATCGGTCATCGGCGGTATTTGTAACAACTCCGGCGGCTCGCTGGTGCAACGTGGCCCGGCGTAT‑ACCGAAATGTCGTTATTcgcgc < 2:246495/90‑1 (MQ=255) gcATCGGTCATCGGCGGTATTTTTAAAAACTCCCGCGGCTCGCTGGGGGAAAGTGGCCCCGCGGAT‑AACGAAAAGTCGTTATTcgcgcgg > 2:614215/1‑89 (MQ=255) gcATCGGTCATCGGCGGTATTTGTAACAACTCCGGCGGCTCGCTGGTGCAACGTGGcc < 2:1158089/58‑1 (MQ=255) gcATCGGTCATCGGCGGTATTTGTAACAACTCCGGCGGCTCGCTGGTGCAACGTGGcc > 1:1158089/1‑58 (MQ=255) gcATCGGTCATCGGCGGGATTTTTAACAACTCCGGCGGCTCGCTGGGGCAACGTGGCCCGGCGGAT‑ACCGAAATGTCGTTCTTtgcgcgg > 2:1427025/1‑89 (MQ=255) gcATCGGTCATCGGCGGGATTTTTAACAACCCCCGCGGCCCGCTGGGGCGACGGGGCCCCGCGGAT‑ACCGAAATGTCGTTATTtgcgcgg > 2:1399583/1‑89 (MQ=255) aTCGGTCATCGGCGGTATTTTTAACAACTCCGGCGGCTCGCTGGGGCAACGTGGCCCCGCGGAT‑ACCGAAATGTCGGTATTCGCGCGgat > 1:490615/1‑90 (MQ=255) aTCGGTCATCGGCGGTATTTGTAACAACTCCGGCGGCTcac < 2:1053431/41‑3 (MQ=255) aTCGGTCATCGGCGGTATTTGTAACAACTCCGGCGGCTcac > 1:1053431/1‑39 (MQ=255) tCGGCGGTATTTGTAACAACTCCGGCGGCTCGCTGGTGCAACGTGGCCCGGCGTAT‑ACCGAAATGTCGTTATTCGCGCGTATAAATGAAg < 1:255433/90‑1 (MQ=255) tCGGCGGTATTTGTAACAACTCCGGCGGCTCGCTGGTGCAACGTGGCCCGGCGTAT‑ACCGAAATGTCGTTATTCGCGCGTATAAATGAAg < 1:1399583/90‑1 (MQ=255) cggTATTTGTAACAACTCCGGCGGCTCGCTGGTGCAACGTGGCCCGGCGTAT‑ACCGAAATGTCGTTATTCGCGCGTATAAATGAAGACgg < 2:1283539/90‑1 (MQ=255) cggTATTTGTAACAACTCCGGCGGCTCGCTGGTGCAACGTGGCCCGGCGTAT‑ACCGAAATGTCGTTATTCGCGCGTATAAATGAAGACgg < 2:333667/90‑1 (MQ=255) gTATTTTTAAAAACTCCCGCGGCTCGCTGGGGCAACGTGGCCCGGCGGAT‑ACCGAAAAGTCGTTTTTTGCGCGGATAAAAGAAAACGGCa > 2:219676/1‑90 (MQ=255) aTTTGTAACAACTCCGGCGGCTCGCTGGTGCAACGTGGCCCGGCGTAT‑ACCGAAATGTCGTTATTCGCGCGTATAAATGAAGACGGCaaa < 2:1396288/90‑1 (MQ=255) | CTGGAAAAAGCCCTCAAACCGCTGGGACGCGAACCGCACTCAGTGATTGGATCATCGTGTATAGGCGCATCGGTCATCGGCGGTATTTGTAACAACTCCGGCGGCTCGCTGGTGCAACGTGGCCCGGCGTAT‑ACCGAAATGTCGTTATTCGCGCGTATAAATGAAGACGGCAAA > NZ_CP009273/2216042‑2216215 |
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: ‑ |