| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 341,039 | 0 | C | A | 26.6% | 182.4 / 26.8 | 65 | intergenic (+337/‑76) | BW25113_RS01680/BW25113_RS01685 | hypothetical protein/hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base C (43/4); new base A (0/17); total (43/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.34e-12 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCA > NZ_CP009273/340952‑341125 | aTAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCaa > 1:295671/1‑90 (MQ=255) tttGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGc > 1:345607/1‑90 (MQ=255) tgtgtgGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCtc > 1:334896/1‑90 (MQ=255) gtgtgGGGTTTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAAAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCtca < 1:112656/90‑1 (MQ=255) gtgtgGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCtca > 1:82980/1‑90 (MQ=255) ggggTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAg > 1:387381/1‑90 (MQ=255) gggTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAgg > 2:449442/1‑90 (MQ=255) gggTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAgg > 2:391592/1‑90 (MQ=255) gtttttttCTGTATCTTCCGCATCGCACCCAAGCCGGAAAAAAAATATGCTGAAAGGCTCATAAAAAAAACCGCCATTAGCTCATCAGGa < 2:243573/86‑1 (MQ=255) ggTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGa > 2:182286/1‑90 (MQ=255) cTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAaga > 2:264965/1‑90 (MQ=255) ttttCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAagag > 2:223295/1‑90 (MQ=255) ttttCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAagag > 2:177715/1‑90 (MQ=255) ttttCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAagag > 1:67129/1‑90 (MQ=255) tGTTTCTTCCGCTCCGCAATCAAACCTGAAAAAAAATACGCAAGACGCCCAAAATTAAAAGCCGCCATTAGCTCATCAGGAAGAGCAgac < 1:474925/90‑1 (MQ=255) tATCTTCCGCACCGCCCTAAAGCCGGAAAAAAAATTTGCTGTAAGCCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacga < 1:264965/90‑1 (MQ=255) tATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacga > 2:337837/1‑90 (MQ=255) tATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacga > 2:291649/1‑90 (MQ=255) tATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacga > 2:466820/1‑90 (MQ=255) tATCTTACGCACCCCACACAACCCTAAAAAAAAATATGCTGAAAGCCTAAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacga < 2:345607/90‑1 (MQ=255) aTCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacgac > 2:431801/1‑90 (MQ=255) tCTTACGCACCGCACTCAGCCCTAAAAAAAAATATGGTGTAGGGCTCATACAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACa < 1:356066/90‑1 (MQ=255) cTTACGCACCGCACTCAAGCCTGAAAAAAAATATGCTGAAAGGCTAATAAAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACaa < 2:12002/90‑1 (MQ=255) ttACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc > 2:356066/1‑90 (MQ=255) ttACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc > 1:164834/1‑90 (MQ=255) ttACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc > 1:425724/1‑90 (MQ=255) ttACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc > 1:204447/1‑90 (MQ=255) tccGCACCCCACCCAAGCCTAAAAAAAAATATGCTGAAAGCCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAcc < 2:425724/88‑1 (MQ=255) tACGCACCGCCCCCAACCCTGAAAAAAAATATGCTGAAGGCCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAcc < 2:60649/90‑1 (MQ=255) cATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGa > 2:61563/1‑90 (MQ=255) cATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGa > 2:178135/1‑90 (MQ=255) accGCACCCAAGCCTAAAAAAAAATATGTGGAAAGCCTAAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAt < 2:359162/88‑1 (MQ=255) tCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCg > 2:48355/1‑51 (MQ=255) tCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt > 1:338731/1‑90 (MQ=255) tCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt > 1:182460/1‑90 (MQ=255) tCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt > 1:301381/1‑90 (MQ=255) tCGCACTCAAGCCTGACAAAAAATATGCTGTAAGGCTCATATAAAAAACCg < 1:48355/51‑1 (MQ=255) ccACTCAAGCCTAAAAAAAAATATGCTGAAAGCCTCATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgt < 2:73661/89‑1 (MQ=255) aagcTCAAGCCTAAAAAAAAATTTGCTGAAAGCCTCAAAAAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgt < 2:333316/87‑1 (MQ=255) tCAAGCCTGAAAAAAAATTTGTTGAAAGGCTAAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTa < 1:424491/90‑1 (MQ=255) cAAGCCAGAAAAAAAATATGTTGAAAGCCTAAAAAAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAt < 2:491216/90‑1 (MQ=255) cAACCCTGAAAAAAAATATGCTGAAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAt < 1:167924/90‑1 (MQ=255) gCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGt > 2:206363/1‑90 (MQ=255) gCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGt > 1:455975/1‑90 (MQ=255) tGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACg > 2:78392/1‑90 (MQ=255) tGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACg > 1:246215/1‑90 (MQ=255) aCAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTAcctg > 1:333316/1‑87 (MQ=255) atatGTTGAAAGCCTAATAAAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAg < 2:387381/90‑1 (MQ=255) tatGTTGTAGGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAgg < 1:223295/90‑1 (MQ=255) tatGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAgg < 1:178135/90‑1 (MQ=255) atGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgttg > 2:297747/1‑70 (MQ=255) atGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGc > 1:24070/1‑90 (MQ=255) atGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGc > 2:276529/1‑90 (MQ=255) atGCTGTAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgttg < 1:297747/70‑1 (MQ=255) gCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCt > 1:87475/1‑90 (MQ=255) gTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGa > 1:324843/1‑90 (MQ=255) aaGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAt < 1:481191/66‑1 (MQ=255) aaGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAt > 2:481191/1‑66 (MQ=255) tCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGt > 1:388261/1‑90 (MQ=255) tCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGt > 2:141343/1‑90 (MQ=255) tCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGt > 1:264195/1‑90 (MQ=255) cATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacga < 1:357817/42‑1 (MQ=255) cATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacga > 2:357817/1‑42 (MQ=255) cATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTc > 1:190468/1‑90 (MQ=255) tatCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCa > 1:298095/1‑90 (MQ=255) | ATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCA > NZ_CP009273/340952‑341125 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |