| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,380,331 | 0 | T | G | 58.8% | 32.6 / 40.1 | 34 | intergenic (+21/‑82) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
| Reads supporting (aligned to +/- strand): ref base T (3/11); new base G (20/0); total (23/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.27e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-06 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAA > NZ_CP009273/2380251‑2380419 | tATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGaaaaa > 1:900614/1‑90 (MQ=255)tATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGaaaaa > 2:592347/1‑90 (MQ=255) gCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAATGcaca > 2:256898/1‑90 (MQ=255) tCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACatt < 1:549404/90‑1 (MQ=255) aTGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGGCAAGGCCATTTTTTTTTTCCt > 2:71136/1‑90 (MQ=255) ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAAGACAAGCACAATTTTTTTTTCCTcc > 2:512611/1‑89 (MQ=255) ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCACAAAAAAAGAAAAACACAATTTTTTTTTCCTcc > 2:1149770/1‑89 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 2:900614/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 2:785688/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 2:286857/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 1:877374/90‑1 (MQ=255) ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGAAAAAAACATTTTATTTTTTCTCTTCTg > 1:135378/1‑90 (MQ=255) aCTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTa < 2:378808/90‑1 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACATTTTTTTTTTTCTCCTCtt‑tttttta > 2:877374/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTG‑tttttaa > 2:375838/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGCACATTTTTTTTTTTCTCCTCTG‑gttttta > 1:467908/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGCACATTTTTTTTTTTCTCATCtt‑tttttaa > 1:93876/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAAAACATTTTTTTTTTTCTCATCTG‑tttttaa > 1:726668/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACATGCACATTTTTTTTTTTCTCCTCTG‑tttttaa > 2:505331/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGAAAATTTTATTTTTTCTCATCTT‑TTTTTAg > 1:892796/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAAAAATTCACATTTTTTTTTTTCTCCTCtt‑tttttat > 2:419146/1‑87 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACATGCACATTTTTTTTTTCCTCCTCtt‑tttttaaa > 2:209156/1‑86 (MQ=255) ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAAGATGCGa < 1:512611/90‑1 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACAATTTTTTTTTTCTCCTCTG‑TTTTTAAATGCGaaaaaa > 1:143381/1‑87 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACAATTTTTTTTTCCTCCTCTT‑TTTTTAAATGAGaaaaaa > 2:5206/1‑87 (MQ=255) gTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGAACATTTTTTTTTTTCTCCTCTT‑TTTTTATATGGGAAAAAAGaa > 1:1045793/1‑90 (MQ=255) gTAAATAATAAAAATTTCTCGAGGCGGGGCCAGGAAAAACGCCAAGCACATTTTTTTTTTCCTCCTCTT‑TTTTTAAATGCGGAAAAAGaa > 2:12632/1‑90 (MQ=255) taaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAACg < 1:797346/90‑1 (MQ=255) taaAAATTTCTCGAGGCGGGTACCGGAAAAAAGACAAGAACATTTTTTTTTTCCTCTTCTGTTTTTAAAATG‑GAAACAAGAAATTaaaag > 1:79874/1‑88 (MQ=255) aaaaTTTCTCGAGGCGGGGGCAGAAAAAAAGACAAGCACATTTTTTTTTTCCTCCTCTT‑TTTTTAAATTGGAAAAAAGAAATTTAAAgga > 2:1007607/1‑90 (MQ=255) aaaaTTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCCTCTT‑TTTTTAGATGGGAAACATGAAATTTAACgga > 2:996180/1‑90 (MQ=255) gAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAACGGAGGAAaaataa < 2:1045793/90‑1 (MQ=255) gTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGaa < 1:419146/90‑1 (MQ=255) gTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGaa < 1:592347/90‑1 (MQ=255) | TATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAA > NZ_CP009273/2380251‑2380419 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |