Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,380,357 | 0 | A | T | 54.0% | 28.6 / 54.7 | 37 | intergenic (+47/‑56) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
Reads supporting (aligned to +/- strand): ref base A (2/15); new base T (20/0); total (22/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.45e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-08 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGC > NZ_CP009273/2380276‑2380435 | aTGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCTGAAAAAAAAACAAGCACaaattttttttctt > 2:1651353/1‑86 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCAGAAAAAAAGACAAACACAATTTTTTTTTTCTCAt > 2:1584864/1‑90 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACATGCACATTTTTTTTTTCCTcct > 1:670656/1‑88 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 1:1060488/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 2:37078/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 2:1507131/90‑1 (MQ=255) ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCTGGAAAAAAGAAATGCACATTTTTTTTTTTCTCATCtt > 2:750746/1‑89 (MQ=255) ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACATGCACCCTTTTTTTTTTCTCctctt > 1:778196/1‑89 (MQ=255) cTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTaa < 1:1584864/90‑1 (MQ=255) cTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAAGAAAAACACAATTTTTTTTTTCTCATCtttttttt > 2:352976/1‑88 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAg < 1:1650453/90‑1 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGGCAAGAACCTTTTTTTTTTCCTCTTCTCttttttg > 1:1313346/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGAAAAAAACAATTTTTTTTTCCTCCTCTGtttttaa > 2:246077/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAAAAAAAAAGAAAAAAACATTTTTTTTTTTCTCCTCTTTTTTTAg > 2:689045/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAAGAAAAACAAATTTTTTTTTTTCTCTTCtttttttaa > 2:593597/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGCACAATTTTTTTTTCCTCTTCtttttttta > 2:1060488/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCCTCTGTTTTTAg > 1:1191226/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAAAAGAACAATTTTTTTTTTCTCATCtttttttaa > 2:1542684/1‑87 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCAAAAAAAAAACAAGCACAATTTTTTTTTCCTCCTCtttttttaaa > 1:662329/1‑86 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAAAAAAAAAGAAAAGAACAATTTTTTTTTTTTTATCTGttttttga > 2:1168071/1‑86 (MQ=255) ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCTGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGa < 1:750746/90‑1 (MQ=255) ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCTCATTTTATTTTCCCTCATCTGTTTTAAGATGCGa < 1:569437/90‑1 (MQ=255) acaacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCTACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGaaa < 1:273856/90‑1 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACATGCACATTTTTTTTTTCCTCATCTGTTTTTAAATGCGaaaaaa > 1:1448437/1‑87 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACATTTTTTTTTTTCTCCTTTTTTTTTAAATGCGaaaaaa > 1:799542/1‑87 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAACACATTTTTTTTTTTCTCCTCCTTTTTTAAATGGGaaaaaa > 2:622352/1‑87 (MQ=255) aGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGa < 1:1020455/90‑1 (MQ=255) aGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGa < 1:540878/90‑1 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAAt < 1:743566/90‑1 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGGCCCAAAAAAAAGAAATGAACAGTTTATTTTTCCTCCTCTTTTTTTAAATGCGAAACATGaaaa > 1:366562/1‑89 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACATGTACAATTTTTTTTTTCTCCTCTTTTTTTAAATGCGGAAAATaaaaa > 1:1718502/1‑89 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAAAAGCACAATTTTTTTTTCCTCCTCTTTTTTTAAAAGGGGAAAATGaaaa > 1:1102002/1‑89 (MQ=255) ataaAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACAATTTTTTTTTCCTCCTTTGTTTTTAGATGGGAAAAAAGAAAATTAAc > 1:40163/1‑90 (MQ=255) taaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACg < 2:1448437/90‑1 (MQ=255) tCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAt < 1:1370657/90‑1 (MQ=255) aGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTa > 1:30264/1‑90 (MQ=255) cGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTAtt < 2:352026/90‑1 (MQ=255) aTGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGc < 2:40163/90‑1 (MQ=255) | ATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGC > NZ_CP009273/2380276‑2380435 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |