| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,380,357 | 0 | A | T | 48.0% | 76.9 / 47.8 | 77 | intergenic (+47/‑56) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
| Reads supporting (aligned to +/- strand): ref base A (6/34); new base T (37/0); total (43/34) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.97e-16 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACT > NZ_CP009273/2380270‑2380446 | tAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACattttatt < 1:3379445/90‑1 (MQ=255)tAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACattttatt < 2:667522/90‑1 (MQ=255) gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAAGAAAAGGACaatttttttt > 2:3299695/1‑90 (MQ=255) tGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAAAACAAGCACaattttttttt > 2:1403052/1‑89 (MQ=255) gATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTccc < 1:5433/90‑1 (MQ=255) tGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGGAAAAAAGACAAGCACAATTTTTTTTTCCTc > 1:828928/1‑90 (MQ=255) ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCAGAAAAAAAGAAAAGAACAATTTTTTTTTCCTct > 2:2201619/1‑89 (MQ=255) ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACCTTTTTTTTTTCCTcc > 1:191490/1‑89 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCAt < 2:1125757/90‑1 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCAt < 1:24175/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 2:1033323/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 1:76790/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 1:1990512/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 1:868511/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 1:1461846/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 2:1569546/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 2:2179289/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 1:2736629/90‑1 (MQ=255) ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCAAAAAAAAGACAAGCAAAATTTTTTTTTTCTCctctt > 1:3376421/1‑89 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGAAAAGCACAATTTATTTTTTCTCATCTGTTTTTAg > 2:254837/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCCGAAAAAAAGACAAGAAAAATTTTTTTTTTCCCATCTGtttttaa > 2:2138774/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCAGAAAAAAAGAAAAGAAAAATTTTTTTTTCCTCATATGtttttaa > 2:2624656/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCAAAAAAAACGAGACGCGCCTTTTTTTTTTCCTCATCTGtttttta > 1:591431/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACATGCACATTTTTTTTTTCCCCATCtttttttaa > 1:1363496/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACAAGCACATTTTATTTTTTCCCATCtttttttta > 1:629886/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAATCACATTTTATTTTTTCTCCTCtttttttaa > 1:166378/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGGACATTTTTTTTTTCCTCCTCtttttttaa > 1:2234009/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACAATTTTTTTTTCCTCATCTGtttttaa > 1:2320218/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCACAAAAAACCAAAAGAAAAATTTTTTTTTTCTCATCTGTTTTTAg > 2:3205971/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAAGAAAAGAACAATTTTTTTTTCCCCCTCtttttttaa > 2:1976913/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGACAAGGACAATTTATTTTTTCTCATTTGTTTTTAg > 1:2312657/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACATGCACATTTTTTTTTTCCTCATCTGtttttaa > 1:2092111/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCCCATTTTTTTTTTCCTCCTCtttttttaa > 1:1930794/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGGACAATTTTTTTTTTCTCCTCtttttttaa > 2:2893049/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAACAAATTTTTTTTTTTCTCATCTGtttttat > 2:1461846/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCAAAAAAAAAAAAAGCAAAATTTTTTTTTTCTCAACTGtttttaa > 2:737162/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTGtttttaa > 1:1651929/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCATCtttttttaa > 1:1212746/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACAATTTTTTTTTCCTCATCTGtttttaa > 1:586229/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACATTTTTTTTTTCCTCTTCTTTTTTTAg > 1:632664/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACAATTTTTTTTTTCTCTTCtttttttaa > 2:2499095/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAAACAAGCACAATTTTTTTTTTCTCCTCTGTTTTTAg > 2:2736629/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAAAAAAAAAGACATGCACCTTTTTTTTTTCCTCCtttttttttaa > 1:877309/1‑87 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAGGACAAGCACAATTTTTTTTTCCTTATCTGtttttaaa > 1:672990/1‑86 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGa < 1:1926334/90‑1 (MQ=255) ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGa < 2:3375584/90‑1 (MQ=255) ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGa < 2:1579780/90‑1 (MQ=255) ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAAAAAAAAGGAAAAGAACATTTTTTTTTTCCTCCTCTTTTTTTAAATGGGa > 1:1362347/1‑90 (MQ=255) acaacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGaaa < 1:889899/90‑1 (MQ=255) caacaGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGCACATTTTTTTTTTCCTCCTCTTTTTTTAAATGGGaaaa > 1:3384284/1‑89 (MQ=255) caacaGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCCTCTTTTTTTAGATGCGaaaa > 1:54168/1‑89 (MQ=255) caacaGAAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGAAAAAGACATTTTTTTTTTCCTCCTCTGTTTTTATATGTGAAAc > 2:889899/1‑90 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACAt < 1:3299695/90‑1 (MQ=255) aGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGa < 1:2368296/90‑1 (MQ=255) aGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGACAAGCACATTTTTTTTTTCCCCCTCTTTTTTTATATGGGGAAAATGa > 1:1440424/1‑90 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAAAAGGACATTTTTTTTTTTCTCCTCTGTTTTTATATGGGAAAAAAGaaaa > 2:784888/1‑89 (MQ=255) taataaAAATTTCTCGAGGCGGGTCCAAAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTTTTTTTAAATGTGAAAAAAGaaaatta > 1:2647003/1‑86 (MQ=255) taataaAAATTTCTCGAGGCGGGGGCCGAAAAAAAGACATGCACAATTTTTTTTTCCTCATCTTTTTTTTTATGGGGAAAATGaaaatta > 1:1618128/1‑86 (MQ=255) ataaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAAc < 2:3376421/90‑1 (MQ=255) ataaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAAc < 2:191490/90‑1 (MQ=255) ataaAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAAAAGCACAATTTATTTTTTCTCTTTTGTTTTTAGATGCGAAACATGAAAATCAAc > 2:940192/1‑90 (MQ=255) aaaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACgg < 2:1578909/90‑1 (MQ=255) aaaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACgg < 2:450656/90‑1 (MQ=255) aaaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACgg < 1:2239405/90‑1 (MQ=255) aaaaTTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACATTTTTTTTTTCCTCATCTGTTTTTAAATGCGGAAAATGAAAATCaaacgg > 1:980668/1‑86 (MQ=255) aaaaTTTCTCGAGGCGGGGCCAGAAAAAACGAAAAGCACATTTTTTTTTTCCTCATCTGTTTTTAAATGCGGAAAATGAAAATTAAcagg > 1:1095418/1‑87 (MQ=255) aTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACggagga < 1:1108774/90‑1 (MQ=255) ttCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGaaa < 2:1882617/90‑1 (MQ=255) ttCTCGAGGCGGGTCCAGAAAAAAAGACATGCACATTTTTTTTTTCCTCCTCTTTTTTTAGAAGCGAAAAATGAAAATCAACGGGGGaaa > 1:2456598/1‑90 (MQ=255) ttCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACAATTTTTTTTTTCTCCTCTGTTTTTAAATGCGAAAAAAAAAAATTAAAGGAGGgaa > 2:1851396/1‑90 (MQ=255) tctcGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGaaaa < 2:1936251/90‑1 (MQ=255) ctcGAGGCGGATCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGaaaaa < 1:986427/90‑1 (MQ=255) cGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATgg < 2:2969731/90‑1 (MQ=255) gTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGaa < 1:2840214/90‑1 (MQ=255) gTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGaa < 1:2035320/90‑1 (MQ=255) gTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGaa < 1:1379060/90‑1 (MQ=255) tCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAt < 2:1840400/90‑1 (MQ=255) aaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGc < 1:817138/90‑1 (MQ=255) aaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGc < 2:54168/90‑1 (MQ=255) aTGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGc < 1:1976913/90‑1 (MQ=255) caTTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGaaa > 2:2753684/1‑90 (MQ=255) ttttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAAtt > 2:1963959/1‑90 (MQ=255) tattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGAc > 2:2918189/1‑90 (MQ=255) attttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACt < 1:2430605/90‑1 (MQ=255) attttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACt < 1:912143/90‑1 (MQ=255) | TAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACT > NZ_CP009273/2380270‑2380446 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |