Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,527,926 | 0 | A | C | 52.2% | 43.4 / 76.7 | 46 | G219G (GGT→GGG) | nanM | N‑acetylneuraminate epimerase |
Reads supporting (aligned to +/- strand): ref base A (21/1); new base C (0/24); total (21/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.17e-12 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-09 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAG > NZ_CP009273/4527799‑4528061 | cGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGG‑TTTTATCCCCTTTATTcac < 2:261714/139‑1 (MQ=255) cTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATTCACAAccgc > 2:942233/1‑139 (MQ=255) cTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATTCACAAccgc > 2:119772/1‑139 (MQ=255) gTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGTTTTGGCTTCGCCATTAATAAGCCAGG‑TTTTATCCCCTTTATTCACAAccg < 1:260746/135‑1 (MQ=255) gagaCGGGAGCAAGTTTATTCCATTTTAAATTATTCCCGGTAAAATCAAGTTCAAATACGGCATCCGTCCGCAATCCTGTTTTGGCTCCGCCATTAATAAGCCAGG‑TTTTATCCCTTTTATTCACAACCGCCGCACCAg < 1:876169/139‑1 (MQ=255) agccGGGACCAAGTTTATCCCATTTTAAATTATTCCCGGTGAAATCAATTTCAAATCCGCCACCCGTCCCCAACCCTGTTTTGCCTTCCCAATTAATAACCCAGT‑TTTTACCCCCTTTATTCACAACCGCCGCACCAg < 1:276967/135‑1 (MQ=255) agaCGGGAGCAAGCTTTTTCCATTTTAAATTATTACCGGTAAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGT‑TTTTATCACCTTTATTCACAACCGCCGCACCAg < 1:342178/138‑1 (MQ=255) agaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGc > 2:851487/1‑139 (MQ=255) agaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGc > 1:267916/1‑139 (MQ=255) agaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGc > 1:311048/1‑139 (MQ=255) agaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGc > 2:164409/1‑139 (MQ=255) aggggAGCAAGTTTATTCCATTTAAAATTATACCCGGGAAAATAAATTAAAAATACGCCATCCTTTCGCAACCCTGTTTTGGCTTCGCCATAAATAAGCCAGT‑TTTTATCCCTTTTATCCACAACCGCCGCACCAGCCg < 1:614544/137‑1 (MQ=255) gggAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGtt > 2:276967/1‑139 (MQ=255) ggaaGCAAGCTTTTTCCTTTTAAAATTATCCCCGTGAAAACCAATTAAAAATACGGCCTCCTTCCCCATCCCTGTTTTGCTTTCCCTATAAAAAAGCCAGT‑TTTTACCCCCTTTTTTCACACCCGCCGCACCAGCCGtt < 2:815254/136‑1 (MQ=255) ggAGCAATTTTATTCCTTTTTAAATTTTTCCCGGGAAAATCAATTACAAATACGGCACCCGTCCGCAATCCTGGTTGGGCTTCGCCTTTAATAACCCTTT‑TTTTACCCCCTTTATCCACACCCGCCGCACCAGCCGTTc < 2:770419/139‑1 (MQ=255) aGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCg > 2:545090/1‑139 (MQ=255) aGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCg > 1:747931/1‑139 (MQ=255) cctttttAAATTATTCCCGGTGAAATAAAGTTAAAATACGGCACCCGTTCGCAATCCTGGTTTGCCTTCCCCATTAAAAAGCCGGT‑TTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGa < 1:164409/135‑1 (MQ=255) ccATTTTAAATTATTCCCGGTAAAACCAACTTCAATTACGCCATCCGTTCCAACCCCTGTTTTGCCTTCCCCATTAAAAACCCATT‑TTTTACCCCTTTTTTCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGa < 1:851487/138‑1 (MQ=255) tttAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc > 1:710699/1‑139 (MQ=255) ttAAATTATTCCCGGTGAAATCAAGTACAAATACGGCATCCGTTCGAAATCCTGTTTTGGCTCCGCCATAAATAAGCCAGT‑TTTTATCCCTTTTATACACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc < 1:913229/138‑1 (MQ=255) taaaaTTATTCCGGGTAAAATAAATTTAAAATAGGGAACCCGTTCGAAACCCTGGTTGGGCTCCGCAATTAAAAAGCAAGG‑TTTTATCCCTTTTATTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc < 1:452511/136‑1 (MQ=255) taaaaTTATTACGGGGAAAATAATGTAAAAATACGCCTTTTGTCCGCACTCCGGGTTTGCCTCCCCCATAAATAAGCCAGG‑TTTAATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc < 2:313160/136‑1 (MQ=255) ttattccGGGGAAAACAATTTTAAAATCCGCCTCCCGTCCCCATCCCTGTTTTGCTTCCCCCATAAAAAGCCCGGTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt < 2:267916/134‑1 (MQ=255) cGGTAAAACAAAGTAAAAATACGGCATCCGTTCGCAATCTTGTTTGGCCTTCGCCATTAATAAGCCAGT‑TTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCa < 2:311048/139‑1 (MQ=255) tGAAATCAATTAAAAACCCGGCTCCCGTCCCAACTCCTGTTTTGTCTTC‑CCATTAAAAAACCAGGTTTTTATCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTg < 2:710699/139‑1 (MQ=255) tGAAATCAAGTTAAAATCCGGCATCCTTTCGCATCCTTGTTTTGGCTTCGCCATTAATAAGCCAGT‑TTTTACCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTg < 1:545090/139‑1 (MQ=255) tGAAATAAAGTTAAAATCCGCCCCCCGTCCGCATTCCGGGTTTGTTTTCGCCATTAATAAGCCAGG‑TTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTg < 1:422490/139‑1 (MQ=255) aaaTCAATTTCAAATACGGCACCCGTTCGCAACCCTGTTTGGTCTCCCCCATTAATAAGCCAGT‑TTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCt < 1:884880/139‑1 (MQ=255) aaaTCAAGTTCAAATCCGGCATCCGTTCCCAATCCTGGTTTGGCTTCGCAATTAAAAAGCAAGT‑TTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCt < 1:556827/139‑1 (MQ=255) aaaTAAAGTTCAAATCCGGCTCCCGTTCGCAATCCTGTTTTGGCTCCCCCATTAATAAGCCAGT‑TTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCt < 1:473707/139‑1 (MQ=255) aaTACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt > 1:125880/1‑139 (MQ=255) aCGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa > 1:56148/1‑138 (MQ=255) gCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa > 2:709658/1‑136 (MQ=255) ccGTTCGCAATCCTGGTTTGGCTCCCCCATTAATAAGCCAGT‑TTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaa < 1:251597/139‑1 (MQ=255) gtccGCACTCCGGTTTGGGCTTCCCAATTAAAAAGCCAGT‑TTTAACCCCCTTTATCCACATCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaaca < 2:807002/136‑1 (MQ=255) ttCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaa > 2:547651/1‑136 (MQ=255) tCCTGGTTTGGCTCCCCCATTAATAAGCCAGT‑TTTTATCCCTTTTATTCGCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTa < 1:79545/139‑1 (MQ=255) gCTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGaaa > 2:561401/1‑138 (MQ=255) gCTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGaaa > 1:919635/1‑138 (MQ=255) gCTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGaaa > 2:776680/1‑138 (MQ=255) aGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGA‑GGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt > 1:945939/1‑139 (MQ=255) ccAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGCGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg > 2:46395/1‑139 (MQ=255) ccAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg > 1:46368/1‑139 (MQ=255) attttttATCCCTTTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc < 2:837649/136‑1 (MQ=255) aTCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAg > 1:716488/1‑139 (MQ=255) | CGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAG > NZ_CP009273/4527799‑4528061 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |