| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,380,341 | 0 | G | A | 53.6% | 2.5 / 31.8 | 28 | intergenic (+31/‑72) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
| Reads supporting (aligned to +/- strand): ref base G (3/10); new base A (15/0); total (18/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.18e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-05 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCT > NZ_CP009273/2380260‑2380428 | aTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGc < 1:770737/90‑1 (MQ=255) cTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACat < 2:621449/90‑1 (MQ=255) cAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACattt < 2:822281/90‑1 (MQ=255) ggATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACaatttttttttc > 2:331846/1‑88 (MQ=255) ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAAGACAAGGACATTTTTTTTTTCCTCa > 2:614446/1‑90 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 2:552132/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 1:101278/90‑1 (MQ=255) ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACATTTTATTTTTCCTCATCTGtt > 2:619860/1‑90 (MQ=255) cTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTaa < 1:1185206/90‑1 (MQ=255) cTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTaa < 1:566789/90‑1 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACAAGCACATTTTTTTTTTCCTCCTCTGtttttaa > 2:796970/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGAAATGCACATTTTTTTTTTTCTCATCTGtttttaa > 1:233853/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACa < 1:400790/58‑1 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACa > 2:400790/1‑58 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTACCGAAAAAGCCACGAGCACATTTTTTTTTTCCTTCGCTGtttttga > 2:783485/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGGACATTTTTTTTTTCCTCCTCtttttttaa > 1:538920/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACATTTTTTTTTTTCTCCTCtttttttaa > 1:32116/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACATTTTTTTTTTCCTCCTCTGtttttaa > 2:887280/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAAAAGGACAATTTATTTTTCCTCCTCCGtttttaa > 1:495791/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAAAAGCACATTTTTTTTTTTCTCCTCTGTTTTTAg > 1:112809/1‑90 (MQ=255) ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAAAAGCACATTTTTTTTTTCCTCCTCTTTTTTTAAATGCGa > 2:62642/1‑90 (MQ=255) caacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAAc > 2:125126/1‑90 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGCGCAGAAAAAATGACATGGACATTTTATTTTCCCTTATCTGTTTTTAGAAGCGaaaaaa > 2:1144405/1‑87 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACAATTTTTTTTTTCTCCTCTGTTTTTAGATGGGAAAAAAGaaaa > 1:1015164/1‑89 (MQ=255) aTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACggagga < 1:783485/90‑1 (MQ=255) tttCTCGAGGCGGGTCCAGAAAAAACGACATGCACATTTTATTTTTCCTCCTCTGTTTTAAGAAGCGGAAAAAGAAAATCAACAGAGGaa > 2:780912/1‑90 (MQ=255) aGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAaaataaa < 2:112809/90‑1 (MQ=255) aaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCt > 2:52605/1‑90 (MQ=255) | ATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCT > NZ_CP009273/2380260‑2380428 |
| 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: ‑ |