Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,826,720 | 0 | C | A | 63.0% | 24.8 / 30.1 | 27 | R12R (CGC→CGA) | xthA | exodeoxyribonuclease III |
Reads supporting (aligned to +/- strand): ref base C (6/4); new base A (9/8); total (15/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCATTAACAACCATCGATCAAATCACTTAACAACAGGCGGTAAGCAACGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGCGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGATATGTTTCCGCTCGAAGAGGTGGCGAAGCTCGGCTACAACGTGTTTTATC > NZ_CP009273/1826584‑1826853 | tcatTAACAACCATCGATCAAATCACTTAACAACAGGCGGTAAGCAACGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTgcgcgc < 2:475779/139‑1 (MQ=255) cATCGATCAAATCACTTAACAACAGGCGGTAAGCAACGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCAcc < 2:259425/139‑1 (MQ=255) cATCGATCAAATCACTTAACAACAGGCGGTAAGCAACGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCAc > 1:221112/1‑138 (MQ=255) ttAACAACAGGCGGTAAGCAACGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGCGCCAGACCTCACCAGCTTGAAGCCAtcg > 1:154248/1‑139 (MQ=255) caGGCGGTAAGCAACGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGaa > 1:44845/1‑137 (MQ=255) aaCGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGa < 2:166685/139‑1 (MQ=255) aCGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGCGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGAt > 1:33167/1‑139 (MQ=255) cgcgAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGAtg > 1:232172/1‑139 (MQ=255) cGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGCGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACaa > 2:3430/1‑138 (MQ=255) tATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGCGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCa < 2:33167/136‑1 (MQ=255) aaTAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGCGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGatat > 2:333577/1‑139 (MQ=255) aaTAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGCGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGatat > 2:475412/1‑139 (MQ=255) aaTAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCa > 1:515408/1‑115 (MQ=255) aaTAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCa < 2:515408/115‑1 (MQ=255) aaTAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGatat > 2:503692/1‑139 (MQ=255) tAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGAt < 1:455190/104‑1 (MQ=255) tAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGAt > 2:455190/1‑104 (MQ=255) gCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGCGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGATATGTTTCCg < 1:333577/138‑1 (MQ=255) ctTTTAATATCAACGGCCTGCGCGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGATATGTTTc > 2:299553/1‑114 (MQ=255) ctTTTAATATCAACGGCCTGCGCGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGATATGTTTc < 1:299553/114‑1 (MQ=255) atatCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGATATGTTTCCGCTCGAAGAGGTGGCGAAGCTCGGCTACa < 2:232172/138‑1 (MQ=255) atatCAACGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGATATGTTTCCGCTCGAAGAGGTGGCGAAGCTCGGCTACa < 2:438977/138‑1 (MQ=255) aaCGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCa > 1:563601/1‑50 (MQ=255) aaCGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCa < 2:563601/50‑1 (MQ=255) aCGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGATATGTTTCCGCTCGAAGAGGTGGCGAAGCTCGGCTACAACg < 2:416804/135‑1 (MQ=255) aCGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGATATGTTTCCGCTCGAAGAGGTGGCGAAGCTCGGCTACAACg > 1:416804/1‑135 (MQ=255) cGGCCTGCGAGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGATATGTTTCCGCTCGAAGAGGTGGCGAAGCTCGGCTACAACGTGttt > 1:87112/1‑139 (MQ=255) cTGCGCGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGATATGTTTCCGCTCGAAGAGGTGGCGAAGCTCGGCTACAACGTGTTTTATc < 1:498932/139‑1 (MQ=255) | TCATTAACAACCATCGATCAAATCACTTAACAACAGGCGGTAAGCAACGCGAAATTCTGCTACCATCCACGCACTCTTTATCTGAATAAATGGCAGCGACTATGAAATTTGTCTCTTTTAATATCAACGGCCTGCGCGCCAGACCTCACCAGCTTGAAGCCATCGTCGAAAAGCACCAACCGGATGTGATTGGCCTGCAGGAGACAAAAGTTCATGACGATATGTTTCCGCTCGAAGAGGTGGCGAAGCTCGGCTACAACGTGTTTTATC > NZ_CP009273/1826584‑1826853 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |