| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 341,039 | 0 | C | A | 32.4% | 78.2 / 12.3 | 37 | intergenic (+337/‑76) | BW25113_RS01680/BW25113_RS01685 | hypothetical protein/hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base C (23/2); new base A (0/12); total (23/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.91e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCA > NZ_CP009273/340960‑341125 | atcatGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGcc > 1:62171/1‑90 (MQ=255)atcatGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGcc > 1:260603/1‑90 (MQ=255) tttGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGc > 2:1767/1‑90 (MQ=255) ttGTGTGGGTTTTTTTTTGTATCTTCCCCACCCCACCAAACCCTAAAAAAAAATATGCTGTAAGGCTCATATAAAAAACCGCCATTAGCt < 1:152564/90‑1 (MQ=255) ttGTGTGGGTTTTTTTCTGTTTCTTACCCACCCCCCCCAAGCCTGAAAAAAAATATGCTGTAAGCCTCATATAAAAAACCGCCATTAGCt < 1:243883/90‑1 (MQ=255) ttGTGTGGGGTTTTTTCTTTATCTTCCCCCCCCCCCCCAACCCTAAAAAAAAATATGCTGTAAGGCTCATATAAAAACCCGCCATTAGCt < 1:209352/90‑1 (MQ=255) ttGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCt > 2:260954/1‑90 (MQ=255) tgtgtgGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCtc > 2:73179/1‑90 (MQ=255) tgtgtgGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCtc > 1:196699/1‑90 (MQ=255) ggggTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAg > 1:238736/1‑90 (MQ=255) gggTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAgg > 1:127855/1‑90 (MQ=255) cTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAAGAAGAGCAga > 2:230971/1‑90 (MQ=255) tATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacga > 1:237504/1‑90 (MQ=255) tATCTTACGCACCGCACTAAAGCCTGAAAGAAAATATGTTGTAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacga < 2:62171/90‑1 (MQ=255) tATCTTACGCACCGCACTAAAGCCGAAAAAAAAATATGCTGAAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacga < 2:237504/90‑1 (MQ=255) ttACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc > 2:289020/1‑90 (MQ=255) ttACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc > 1:83496/1‑90 (MQ=255) ttACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc > 2:164839/1‑90 (MQ=255) ttACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc > 2:225315/1‑90 (MQ=255) tACGCACCCCACTCAACCCTAAAAAAAAATATGCTGAAAGGCTAATAAAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAAcc < 1:120801/90‑1 (MQ=255) tACCCACCGCACACAACCCCGAAAAAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAcc < 2:248921/90‑1 (MQ=255) gCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAgg < 1:293348/90‑1 (MQ=255) accGCACTAAAGCCTGAAAAAAAATATGCTGTAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAt < 2:238736/88‑1 (MQ=255) tCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt > 1:187107/1‑90 (MQ=255) tCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt > 2:120801/1‑90 (MQ=255) gCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgt > 2:34952/1‑90 (MQ=255) gCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgt > 2:205273/1‑90 (MQ=255) gCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgt > 1:110667/1‑90 (MQ=255) gCACTCAACCCTGAAAAAAAATATGCTGAAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgt < 2:110667/90‑1 (MQ=255) aCTAAAGCCTGAAAAAAAATATGCTGAAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgttg < 2:237785/90‑1 (MQ=255) aaGCCTAAAAAAAAATATGTTGAAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATg < 2:139922/90‑1 (MQ=255) aaaaTATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAt > 1:57303/1‑70 (MQ=255) aaaaTATGCTGTAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAt < 2:57303/70‑1 (MQ=255) aaTATGCTGAAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGa < 2:113079/90‑1 (MQ=255) cATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTc > 1:236342/1‑90 (MQ=255) cATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTc > 1:118419/1‑90 (MQ=255) tatCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCa > 2:301346/1‑90 (MQ=255) | ATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCA > NZ_CP009273/340960‑341125 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |