Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,380,341 | 0 | G | A | 56.4% | ‑1.8 / 49.9 | 39 | intergenic (+31/‑72) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
Reads supporting (aligned to +/- strand): ref base G (1/16); new base A (22/0); total (23/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.51e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-08 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTA > NZ_CP009273/2380260‑2380429 | aTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGc < 2:466565/90‑1 (MQ=255) aTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGc < 1:1138791/90‑1 (MQ=255) aTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGc < 2:170188/90‑1 (MQ=255) 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:540878/90‑1 (MQ=255) aGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGa < 1:1020455/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) | ATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTA > NZ_CP009273/2380260‑2380429 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |