| 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.8% | ‑1.7 / 14.0 | 13 | intergenic (+31/‑72) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
| Reads supporting (aligned to +/- strand): ref base G (0/6); new base A (7/0); total (7/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.83e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCG‑ACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAA > NZ_CP009273/2380252‑2380419 | aTCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAg < 2:450954/90‑1 (MQ=255) cAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCG‑ACATGCACattt < 2:298993/90‑1 (MQ=255) gATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCG‑ACATGCACATTTTATTTTccc < 2:800053/90‑1 (MQ=255) ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAG‑ACAAGCACAATTTTTTTTTCCTcc > 1:2463/1‑89 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAACGCAAAGGCAAATTTTTTTTTCCCTCTCCtttttttt > 2:786771/1‑88 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACG‑ACATGCACATTTTATTTTTCCTCATCTGtttttaa > 1:42909/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAG‑ACAAAAACATTTTATTTTTTCTCATCTGtttttta > 2:840227/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCACAAAAAACG‑ACATCCACAATTTTTTTTTCCCCCTCtttttttaa > 1:435228/1‑87 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAG‑AAATGCACATTTTTTTTTTCCTCATCTTTTTTTATATGCGAAACaa > 2:375385/1‑89 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACG‑ACATGCACATTTTTTTTTTCCTCATCTGTTTTTAAATGCGGAAAATGaaaa > 1:360898/1‑89 (MQ=255) taaAAATTTCTCGAGGCGGGTCCAGAAAAAGCG‑ACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACg < 2:679906/90‑1 (MQ=255) aaaaaTTTCTCGAGGCGGGTCCAGAAAAAGCG‑ACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACgg < 1:118213/90‑1 (MQ=255) ggCGGGTCCAGAAAAAGCG‑ACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAaaataaat < 1:84331/90‑1 (MQ=255) gTCCAGAAAAAGCG‑ACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGaa < 2:215639/90‑1 (MQ=255) | ATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCG‑ACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAA > NZ_CP009273/2380252‑2380419 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |