Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,380,341 | 0 | G | A | 59.1% | ‑2.0 / 18.1 | 22 | intergenic (+31/‑72) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
Reads supporting (aligned to +/- strand): ref base G (2/7); new base A (13/0); total (15/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.11e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.61e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAAC > NZ_CP009273/2380260‑2380399 | aTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGc < 2:512053/90‑1 (MQ=255) ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGACGCGGGTCCAGAAAAAGCGACATGcac > 2:467439/1‑90 (MQ=255) tCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACatt < 1:175409/90‑1 (MQ=255) cAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACaatt > 2:74905/1‑90 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 1:459670/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 2:161267/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 2:331464/90‑1 (MQ=255) ggtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 1:63877/89‑1 (MQ=255) ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAAAACATTTTTTTTTTTCTCATCtt > 1:356170/1‑89 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCCTCTG‑tttttaa > 2:81657/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACATTTTTTTTTTCCTCCTCtt‑tttttaa > 2:562516/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACAATTTTTTTTTTCTCCTCtt‑tttttaa > 2:157004/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGCACAATTTATTTTTCCTCCTCTG‑tttttaa > 1:594200/1‑87 (MQ=255) ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAAGATGCGa < 2:828572/90‑1 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACAAGGAAATTTTATTTTTTCTCATCCG‑CTTTTAGATGAGaaaaaa > 1:384091/1‑87 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACAATTTTTTTTTCCTCATCTT‑TTTTTAAATGCGgaaaaa > 2:242256/1‑87 (MQ=255) aGTAAATAATAAAAATTTCTCGAGGCGGGGGCAGAAAAAAAGACAAGCACAATTTTTTTTTTCTCCTCTT‑TTTTTAGATGGGAAACaaaa > 1:637946/1‑87 (MQ=255) aGTAAATAATAAAAATTTCTCGAGGCGGGGGCAGAAAAAAAAACAAGCAAATTTTTTTTTTTCCCCTCTG‑TTTTTAGATGCGGAAAATGa > 1:332376/1‑90 (MQ=255) aGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACAATTTTTTTTTCCTCATTTGTTTTTAAAATGGGaaaaaaa > 2:652876/1‑86 (MQ=255) taataaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACAAGCACATTTTATTTTTCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCa > 1:57638/1‑90 (MQ=255) ataaAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTG‑TTTTTAGATGCGAAAAATGAAAATCaaa > 2:372841/1‑89 (MQ=255) ataaAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACATTTTTTTTTTCCTCCTCTT‑TTTTTAGATGCGGAAAATGAAAATTAAc > 2:692851/1‑90 (MQ=255) | ATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAAC > NZ_CP009273/2380260‑2380399 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |