| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,380,347 | 0 | T | A | 56.5% | 10.9 / 26.8 | 23 | intergenic (+37/‑66) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
| Reads supporting (aligned to +/- strand): ref base T (3/7); new base A (13/0); total (16/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.89e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.06e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTG > NZ_CP009273/2380260‑2380434 | 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) cATGCACATTTTATTTTCCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTg > 1:558145/1‑90 (MQ=255) | ATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTG > NZ_CP009273/2380260‑2380434 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |