Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,527,926 | 0 | A | C | 55.0% | 17.9 / 30.2 | 20 | G219G (GGT→GGG) | nanM | N‑acetylneuraminate epimerase |
Reads supporting (aligned to +/- strand): ref base A (9/0); new base C (0/11); total (9/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.95e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAG > NZ_CP009273/4527814‑4528061 | ggaCGGGATCATTTTTATTCCATTTTAAATATTTACGGGGAAAATAAAGTTAAATCACGGCTTCCGTTCGCAACCCTGTTTTGGCTCCGCAATTAATACGCCGGGTTTTACCCCTTTTTTCCACAA‑CCGCCGCACCAGc < 1:557812/138‑1 (MQ=255) agaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAA‑CCGCCGCACCAGc > 1:518178/1‑139 (MQ=255) gggAGCAAGCTTATTCCTTTTTAAATTATTACCGGTAAAATCAAGTTCAAATACGGCACCCGTTCACAATCCTGGTTTGGCTCCGCCATTAAAAAGCCAGTTTTTATCCCCTTTATTCACAA‑CCGCCGCACCAGCCGtt < 1:352915/139‑1 (MQ=255) gCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAA‑CCGCCGCACCAGCCGTTCCGTAc > 1:476397/1‑137 (MQ=255) ttCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGc > 2:524655/1‑138 (MQ=255) tttAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc > 1:21916/1‑139 (MQ=255) taaatttattcCCGGTAAAACAAAGTAAAAATCCGCCATCCGTCCGAAATCCTGTTTTGGCTCCCCCATAAATAAGCCGGTTTTTATCCCTTTTTTCCACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc < 2:21916/133‑1 (MQ=255) taaaaTTATTCCGGTTAAAAGCAAGTTAAAATCCGGCATCTGTCCGCAATCCTGTTTTGCCTCCCCCATTAATAGGCCGGTTTTTACCCCCTTTATCCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc < 2:265807/136‑1 (MQ=255) ttaCCGGTGAAATAAAGTTCAAATACGGCATCCTTTCGCAACCCTGGTTGGGCTTCGCCATTAAAAAGCCAGGTTTTATCCCCTTTATCCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTaa < 1:186875/138‑1 (MQ=255) cGGTGAAATCAATTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCCCCATAAATAAGCCGGGTTTTATCCCCTTTATTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCa < 1:411646/139‑1 (MQ=255) cGGGGAAATCAATTTCAAATACGCCATCCTTCCGCAACCCTGGTTTGGCTCCCCCATAAATAAGCCAGGTTTTACCCCCTTTATCCCCAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCa < 1:525112/139‑1 (MQ=255) aaaTCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATAAATAAGCCAGGTTTTACCCCCTTTATCCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCt < 1:139863/139‑1 (MQ=255) aaTACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt > 1:415942/1‑139 (MQ=255) ccATCCGTCCGCAATCTTGTTTTGGCTCCGCCATTAATAAGCCGGGTTTTATCCCCTTTATCCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAg < 1:436504/138‑1 (MQ=255) tCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac > 1:439986/1‑138 (MQ=255) ggttttGTCTCCGCCATTAATAAGCCGGTTTTTATCCCTTTTATCCACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg < 2:476397/136‑1 (MQ=255) ccAGGTTTTATCACCTTTATTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg > 1:512940/1‑139 (MQ=255) caagTTTTATCCCCTTTTATCCCAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg < 2:518178/136‑1 (MQ=255) gTTTTATCACCTTTATTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa > 2:510935/1‑138 (MQ=255) aTCACCTTTATTCACAA‑CCGCCGCACCAGCCGTTCCGGACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAg > 2:225997/1‑139 (MQ=255) | AGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAG > NZ_CP009273/4527814‑4528061 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |