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.6% | 7.3 / 26.9 | 19 | G219G (GGT→GGG) | nanM | N‑acetylneuraminate epimerase |
Reads supporting (aligned to +/- strand): ref base A (9/0); new base C (0/10); total (9/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.08e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTT‑CAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTG > NZ_CP009273/4527792‑4528052 | ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTT‑CAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCtt > 1:607073/1‑139 (MQ=255) gggAGCAATCTTTTTCCTTTTAAATTTATTCCCGGGAAAATCAATTT‑AAAACCCGGCATCCTTCCGCAATCCTGGTTTGGCTCCCCCATTAAAAACCCAGTTTTTACCCCCTTT‑ATCCACAACCGCCGCACCAGCCGtt < 2:607073/139‑1 (MQ=255) gggAGCAAGTTTATTCCATTTTAAATTATTACCGGTAAAATAAAGTT‑CAAATCCGGCACCCGTTCGAAACCCTGGTTTGGCTCCGCCATAAATAAGCCAGGTTTTATCCCCTTT‑ATTCACAACCGCCGCACCAGCCGtt < 1:469241/139‑1 (MQ=255) gggACCAAGCTTATCCCATTTTAATTTATGCCCGGAGAAAAAAACTA‑AAAACACGACCCCCTTCCGCAATCCTGGTTGGGCTCCGCCATTAAAAGGCCAGTTTTTATCCCTTTA‑ATCCCCACCCGCCGCACCAGCCGtt < 2:150553/139‑1 (MQ=255) gCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTT‑CAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTAc > 1:49887/1‑137 (MQ=255) cctttttAAATTATTCCCGGTGAAATAAATTT‑CAAATCCGCCACCCTTCCGCAATCCTGGTTTGGCTTCGCCTTTAATAACCCAGTTTTTATCCCTTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGa < 1:577052/135‑1 (MQ=255) ttttAAATTATTACCGGTGAAATCAAGTT‑CAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGctgtctc > 2:2897/1‑133 (MQ=255) taaaaTTATTACCGGTAAAATAAAGTA‑AAAATACGGCACCCTTGCGCAATCCTGTTTTGGCTCCGCCATAAATAAGCCATGTTTTACCCCCTTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc < 1:810632/136‑1 (MQ=255) taaaaTTATT‑CCGGGAAAAATAAGTTCCAAATCCGTACTCCGTCCCCAATCTTGGTTTGCCTCCGCCATAAATAGGCCGGGTTTT‑TCCCCTTTAATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc < 2:316811/136‑1 (MQ=255) ttcccGGTAAAACCAAGTT‑AAAATCCGCCATCCGTCCGCATTCCTGGTTTGCCTCCGCCATAAATAACCCAGTTTTTATCCCCTTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTaa < 2:423659/135‑1 (MQ=255) gTGAAATCAAGTT‑CAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt > 1:423659/1‑139 (MQ=255) aaaTCAATTA‑CAATTACGGCATCCGTTCGCAATCCTGGTTTGGCTCCGCCTTTAATAACCCAGTTTTTATCCCCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCt < 1:694014/139‑1 (MQ=255) ttt‑CAAATCCGCCATCCGTTCACAACCCTGGTTTGGCTTCCCTTTTAAAAAGCCAGTTTTTATCCCTTTT‑ATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGa < 1:394374/138‑1 (MQ=255) tt‑CAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCCTTGCTGCGTTGAg > 2:182986/1‑139 (MQ=255) ggCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa > 2:25025/1‑137 (MQ=255) ggTTTGGCTCCCCCTTAAATAAGCCTGTTTTTATCCCTTTT‑ATCAACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg < 2:694236/138‑1 (MQ=255) gCTTCGCCATTAATAAGCCAGGTTTTATCACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGaaa > 2:596541/1‑138 (MQ=255) gCTTCGCCATTAATAAGCCAGGTTTTATCACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGaaa > 2:102533/1‑138 (MQ=255) ccAAGTTTTATCACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg > 1:126389/1‑139 (MQ=255) | CCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTT‑CAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTG > NZ_CP009273/4527792‑4528052 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |