Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,380,341 | 0 | G | A | 58.3% | ‑2.1 / 11.6 | 12 | intergenic (+31/‑72) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
Reads supporting (aligned to +/- strand): ref base G (1/4); new base A (7/0); total (8/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.01e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCT‑GTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATT > NZ_CP009273/2380269‑2380424 | gTAGTGGATGGTGTT‑TAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACattttatt < 2:364628/90‑1 (MQ=255) gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGACGAAAAAAAAAACAAGCACattttttttt > 1:877727/1‑90 (MQ=255) tggtataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 1:422380/87‑1 (MQ=255) ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAACACAATCAAATTTTTTTTTTACTCctctt > 2:394858/1‑89 (MQ=255) agttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCT‑GTTTTaa < 1:140357/88‑1 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGACAAGGACAATTTTTTTTTTCTCCTCT‑Gtttttat > 1:537503/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCAGAAAAAAAGACAAGCACATTTTTTTTTTCCTCCTCt‑ttttttaa > 2:875021/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGGAAAAACGACCAGGACATTTTTTTTTTCTTCCTCTGGTTTTaa > 2:665140/1‑90 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCATCT‑GTTTTTAGATGCGAAAcct > 2:140357/1‑88 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACATGCACATTTTTTTTTTCCTCCTCT‑TTTTTTAAATGGGaaaaaa > 2:655004/1‑87 (MQ=255) cagaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCT‑GTTTTAAGATGCGAAACATGAAATTCAACg < 2:140632/87‑1 (MQ=255) gAAAAAGCGACATGCACATTTTATTTTCCCTCATCT‑GTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGAtt > 2:689384/1‑90 (MQ=255) | GTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCT‑GTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATT > NZ_CP009273/2380269‑2380424 |
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: ‑ |