| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,380,362 | 0 | C | T | 55.0% | ‑2.0 / 23.4 | 20 | intergenic (+52/‑51) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
| Reads supporting (aligned to +/- strand): ref base C (3/6); new base T (11/0); total (14/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.17e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-04 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACT > NZ_CP009273/2380275‑2380446 | gATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTccc < 1:577389/90‑1 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCAGAAAAAAAGAAAAGCAC‑ATTTTTTTTTTCCtttc > 2:71098/1‑87 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATc < 1:950891/90‑1 (MQ=255) ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCAGAAAAAACGACATGCAC‑ATTTTTTTTTTCCTCCTCtttt > 2:548255/1‑90 (MQ=255) ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCAGAAAAAAAGAAAAGAAA‑AATTTTTTTTTTCTCATCTGtt > 1:687241/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCAAAAAAAACG‑CAAGAAC‑CCTTTTTTTTTCCCCATCtttttttaa > 2:857756/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCCGAAAAAACGACAAGCAC‑ATTTTTTTTTTCCTCATCCGTTTTTAg > 2:104557/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCCGAAAAAAAGAAAAGCAC‑ATTTTTTTTTTTCTCCTCTGtttttat > 1:283523/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCAGAAAAAACGACAAGGAC‑ATTTTTTTTTTCCTCCTCttttttttg > 2:890676/1‑87 (MQ=255) caacaGTAAATAATAAAAATTTCTCGAGGCGGGG‑GCCGAAAAAAAAACAAGCAC‑AATTTTTTTTTCCTCCTCTGTTTTTAAATGGGaaaa > 1:725209/1‑89 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCCGAAAAAAAGACAAGCACAATTTTTTTTTCCCCCTTTTTTTTTAA‑ATGGGaaaaaa > 2:649050/1‑87 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGG‑CCCGAAAAAACGACATGCAC‑ATTTTTTTTTTCCTCATCTGTTTTTAAATGCGGAAAAAGAAAt > 2:1235508/1‑90 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGG‑CCCGAAAAAACGACAAGCAC‑ATTTTTTTTTTCCTCATCTGTTTTTAAATGCGGAAAAAGaaaa > 2:138038/1‑89 (MQ=255) ataaAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAAc < 1:28095/90‑1 (MQ=255) ttCTCGAGGCGGGG‑CCAGAAAAAACGACATGCAC‑ATTTTTTTTTTCCTCCTCTTTTTTTAAATGGGAAAAAAAAAAATCAAAGGGGGgaa > 2:1302250/1‑90 (MQ=255) aGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAaaataaa < 2:1240825/90‑1 (MQ=255) aaaGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGc < 2:1100269/90‑1 (MQ=255) aaaGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGc < 1:104557/90‑1 (MQ=255) aaGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCt > 1:637002/1‑90 (MQ=255) ttttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAAtt < 1:548255/90‑1 (MQ=255) attttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACt > 2:988684/1‑90 (MQ=255) | GATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACT > NZ_CP009273/2380275‑2380446 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |