Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 341,015 | 0 | G | A | 26.7% | 137.6 / 16.1 | 60 | intergenic (+313/‑100) | BW25113_RS01680/BW25113_RS01685 | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (39/5); new base A (0/16); total (39/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.36e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGG > NZ_CP009273/340928‑341101 | gttgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGaa > 2:280642/1‑90 (MQ=255) gACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAat > 2:395649/1‑90 (MQ=255) aCATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAata > 2:241365/1‑90 (MQ=255) cATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAatat < 2:402685/90‑1 (MQ=255) cATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAatat < 1:17092/90‑1 (MQ=255) ttCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGc < 1:241365/90‑1 (MQ=255) cATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTg > 1:12002/1‑90 (MQ=255) cATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTg > 2:114143/1‑90 (MQ=255) cATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTg > 1:462247/1‑90 (MQ=255) tatGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTa < 1:66083/90‑1 (MQ=255) atGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTaa < 2:93149/90‑1 (MQ=255) aaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCat > 1:36191/1‑90 (MQ=255) aaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCat > 1:420362/1‑90 (MQ=255) 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:466820/1‑90 (MQ=255) tATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacga > 2:337837/1‑90 (MQ=255) tATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacga > 2:291649/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 > 1:164834/1‑90 (MQ=255) ttACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc > 1:204447/1‑90 (MQ=255) ttACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc > 1:425724/1‑90 (MQ=255) ttACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc > 2:356066/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 > 1:246215/1‑90 (MQ=255) tGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACg > 2:78392/1‑90 (MQ=255) aCAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTAcctg > 1:333316/1‑87 (MQ=255) | GTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGG > NZ_CP009273/340928‑341101 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |