| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 728,670 | 0 | A | C | 73.0% | 46.3 / 28.7 | 37 | D1211A (GAC→GCC) | rhsC | rhs element protein RhsC |
| Reads supporting (aligned to +/- strand): ref base A (7/3); new base C (14/13); total (21/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.61e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCA > NZ_CP009273/728540‑728803 | gCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGCCCCGCTGc < 2:289195/139‑1 (MQ=17)gCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTg > 2:595897/1‑138 (MQ=32) cAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGCCCCGCTGCa < 1:550844/139‑1 (MQ=17) aTGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGCCCCGCTGCAGGGGCGatatat > 2:285739/1‑138 (MQ=21) aTGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGCCCCGCTGCAGGGGCGatatat > 2:417504/1‑138 (MQ=21) aTGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGCCCCGCTGCAGGGGCGatatat > 2:529019/1‑138 (MQ=21) aTGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGAcc < 1:512984/119‑1 (MQ=32) aTGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGAcc > 2:512984/1‑119 (MQ=32) ggttAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGCCCCGCTGCAGGGGCGATATATCACt < 2:270794/136‑1 (MQ=18) ctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGCCCCGCTGCAGGGGCGATATATCACTCAGGATCCGa > 2:109846/1‑139 (MQ=17) tgaatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGCCCCGCTGCAGGGGCGATATATCACTCAGGATCCGAt > 1:487422/1‑139 (MQ=17) aaatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGCCCCGCTGCAGGGGCGATATATCACTCAGGATCCGAtt < 2:64503/138‑1 (MQ=17) aaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCGGCCGGGGCAGCAGTATGAGGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGCCCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCt < 2:564754/139‑1 (MQ=12) agaACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGa > 1:494601/1‑138 (MQ=32) gCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGCCCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGa < 1:417504/139‑1 (MQ=21) cTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGCCCCGCTGCAGGGGCGATATATCACTCAGGATCCGa > 2:312637/1‑115 (MQ=17) cTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGCCCCGCTGCAGGGGCGATATATCACTCAGGATCCGa < 1:312637/115‑1 (MQ=17) ccGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGACAGGGGGATGGAATTTTTATCAGTATCCg > 1:533414/1‑139 (MQ=21) cTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGCCCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTg > 2:83232/1‑138 (MQ=21) tgatgaGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGCCCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAata > 1:120558/1‑137 (MQ=255) aggagTCCGGCCTGTATTACAACCGCCACCGCTATTATGCCCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATc > 2:428967/1‑138 (MQ=255) gagTCCGGCCTGTATTACAACCGCCACCGCTATTATGCCCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCAt < 2:360819/139‑1 (MQ=255) agTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGAt < 1:610243/121‑1 (MQ=38) agTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGAt > 2:610243/1‑121 (MQ=38) tCCGGCCTGTATTACAACCGCCACCGCTATTATGCCCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAg > 2:530291/1‑139 (MQ=255) tCCGGCCTGTATTACAACCGCCACCGCTATTATGCCCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAg > 1:381927/1‑139 (MQ=255) gCCTGTATTACAACCGCCACCGCTATTATGCCCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGAtt < 1:585966/139‑1 (MQ=255) gCCACCGCTATTATGCCCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATg > 1:102742/1‑139 (MQ=255) gCCACCGCTATTATGCCCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATg > 2:489505/1‑139 (MQ=255) gCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCg > 2:477253/1‑98 (MQ=38) gCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCg < 1:477253/98‑1 (MQ=38) gCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATg > 2:50127/1‑139 (MQ=255) cGCTATTATGCCCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTa < 1:530291/139‑1 (MQ=255) ttatGCCCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGAt < 2:589120/103‑1 (MQ=255) ttatGCCCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGAt > 1:589120/1‑103 (MQ=255) ttatGCCCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCa < 1:489505/139‑1 (MQ=255) ttatGCCCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCa < 2:487422/139‑1 (MQ=255) | GCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGGCTGAAAGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGATCTCAAATATAGATCCATTAGGATTAGAAACACTAAAATGCATTAAGCCA > NZ_CP009273/728540‑728803 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |