Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,380,362 | 0 | C | T | 51.8% | 7.7 / 37.5 | 27 | intergenic (+52/‑51) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
Reads supporting (aligned to +/- strand): ref base C (8/5); new base T (14/0); total (22/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.59e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.03e-05 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑A‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAA > NZ_CP009273/2380273‑2380451 | tGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑A‑TTTTATTTTc < 1:1017782/90‑1 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCAA‑A‑ATTTTTTTTTCCTCAt > 2:989389/1‑90 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑A‑TTTTATTTTCCCTCATc < 1:1422190/90‑1 (MQ=255) ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCAC‑A‑TTTTTTTTTTCCTCATCtttt > 2:1246247/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACATACACAA‑TTTTTTTTTCCCCCTTTTTTTTTaa > 2:356477/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCAC‑A‑TTTTTTTTTTCCTCATCtttttttaa > 2:1057100/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAAAAAAAACGACAAGCAC‑A‑TTTTTTTTTTCCTCCTCTGtttttaa > 2:1602514/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCCGAAAAAAAGAAAAGCAC‑A‑ATTTTTTTTTTCTCATCTGttttttt > 1:1517736/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAAGAAAAACAC‑A‑ATTTTTTTTTTCTCCTCtttttttaa > 1:1510369/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGAAAAGCAC‑A‑TTTTTTTTTTTCTCATCtttttttaa > 2:285053/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGCAC‑A‑TTTTTTTTTTCCTCCTCtttttttaa > 2:334653/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGAAC‑A‑TTTTTTTTTTCCTCATCTGtttttaa > 2:1323797/1‑87 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCAAAAAAAAAGACAAGCAC‑A‑TTTTTTTTTTCCCCCTCTTTTTTTAAATGCGaaaaaa > 2:1223496/1‑87 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGAAAAGAAC‑A‑TTTTTTTTTTTCTCCTCTGTTTTTAGAAGGGAAACAAGaaaa > 1:957771/1‑89 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAAACACAGAC‑ATTTTTTTTTTTCCTC‑TTTGTTTTTAAATAGGAAAGATGaaaa > 1:1564917/1‑89 (MQ=255) aaaaTTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCAC‑A‑TTTTTTTTTTCCTCCTCTTTTTTTAAATGGGAAAAAAGAAAATCAACgga > 2:418572/1‑90 (MQ=255) aaaaTTTCTCGAGGCGGGGCCAGAAAAAGCCACATTTCC‑A‑TTTTATTTTTCCTCGTCTGTTTTAATAAGGGGAAAATCAAATTCAACcga > 2:1584276/1‑90 (MQ=255) aGGCGGGTCCAGAAAAAGCGACATGCAC‑A‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAaaataaa < 1:704870/90‑1 (MQ=255) ggTCCAGAAAAAGCGACATGCAC‑A‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGa < 2:685036/90‑1 (MQ=255) tttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTg > 1:109950/1‑90 (MQ=255) tttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTg > 1:18164/1‑90 (MQ=255) tattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGAc < 2:1565601/90‑1 (MQ=255) ttttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTc > 2:859965/1‑90 (MQ=255) tttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCa > 1:325042/1‑90 (MQ=255) tCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa < 2:1510369/90‑1 (MQ=255) tCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa > 2:1323750/1‑90 (MQ=255) tCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa > 1:59096/1‑90 (MQ=255) tCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa > 2:832931/1‑90 (MQ=255) cccTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGaa < 2:20435/90‑1 (MQ=255) | TGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑A‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAA > NZ_CP009273/2380273‑2380451 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |