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.6 / 36.9 | 28 | intergenic (+31/‑72) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
Reads supporting (aligned to +/- strand): ref base G (4/9); new base A (15/0); total (19/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.04e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-07 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGAC‑ATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTA > NZ_CP009273/2380252‑2380429 | aTCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAg < 2:88385/90‑1 (MQ=255) aTCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAg < 2:1162973/90‑1 (MQ=255) cTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGAC‑ATGCACat < 2:991215/90‑1 (MQ=255) cAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAC‑AAGCACattt > 2:2268970/1‑90 (MQ=255) ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAC‑AAGCACaatttt > 2:1372404/1‑89 (MQ=255) gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAC‑AAGCACattttttttt > 2:1914670/1‑90 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGGAACATTTTTTTTTcctcgt > 2:1166954/1‑86 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGAC‑ATGCACATTTTATTTTCCCTCATc < 2:1657211/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGAC‑ATGCACATTTTATTTTCCCTCATc < 1:915470/90‑1 (MQ=255) ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAC‑ATACACATTTTTTTTTTTCTCATCTGtt > 1:1136474/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCGGAAAAAAAGAG‑AGGCACATTTTTTTTTTTCTCCTCTGTTTTTAg > 2:915470/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGAC‑ATGCACATTTTTTTTTTCCTCATCTGtttttaa > 2:1458789/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGAA‑AAGCAAAATTTTTTTTTTCTCATCTGtttttaa > 1:961257/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAC‑ATGCACATTTTTTTTTTCCTCCTCTGTTTTTAg > 2:1023497/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGGAAAAACGGC‑ACGCACATTTTTTTTTTCCTCACCTGTTTTaat > 2:1352489/1‑89 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAAAAAAAAAAAC‑AAGCACAATTTTTTTTTTCTCCTCTGtttttat > 1:1641871/1‑87 (MQ=255) acaacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGAC‑ATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGaaa < 2:2329034/90‑1 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAC‑AAGCACATTTTTTTTTTTCTCCTCTGTTTTTAGATGGGaaaaaa > 1:2254905/1‑87 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAA‑AAGCACATTTTTTTTTTTCTCCTCTGTTTTTAGATTGGaaaaaa > 1:1211608/1‑87 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAA‑AAGCACATTTTTTTTTTTCTCCTCTTTTTTTATATGGGAAAAAAGaaaa > 1:2777980/1‑89 (MQ=255) ataaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGAC‑ATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAAc < 1:1166954/90‑1 (MQ=255) tttCTCGAGGCGGGGCCAGAAAAAACGAC‑AAGCACATTTTTTTTTTCCTCCTCTGTTTTTAAAAGCGAAAAAAGAAAATTAAAGGAGGaa > 1:36057/1‑90 (MQ=255) ctcGAGGCGGGTCCAGAAAAAGCGAC‑ATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGaaaaa < 2:2254905/90‑1 (MQ=255) gCGGGTCCAGAAAAAGCGAC‑ATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATg < 1:2079079/90‑1 (MQ=255) aaaGCGAC‑ATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGc < 1:1668041/90‑1 (MQ=255) aaaGCGAC‑ATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGc < 2:30906/90‑1 (MQ=255) aaGCGAC‑ATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCt > 1:96748/1‑90 (MQ=255) aaGCGAC‑ATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCt > 1:589863/1‑90 (MQ=255) aaGCGAC‑ATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCt > 1:1699848/1‑90 (MQ=255) aGCGAC‑ATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTa > 1:2007634/1‑90 (MQ=255) | ATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGAC‑ATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTA > NZ_CP009273/2380252‑2380429 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |