Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,380,362 | 0 | C | T | 44.4% | 14.1 / 10.3 | 18 | intergenic (+52/‑51) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
Reads supporting (aligned to +/- strand): ref base C (2/8); new base T (8/0); total (10/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.05e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACT > NZ_CP009273/2380275‑2380446 | gATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTccc < 2:291955/90‑1 (MQ=255) ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGGCCAGAAAAAAAGAAAAGCACAATTTTTTTTTCCTcc > 1:191881/1‑89 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCAt < 1:1017915/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 2:70421/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 1:789314/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 1:856110/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 2:1021834/90‑1 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCCGGAAAAAG‑GACAAACACACTTTATTTTTTCTCATCTGttttttt > 2:856110/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCCGAAAAAAAGAAAAGAACAATTTTTTTTTTCTCCTCTGTTTTTAg > 2:1057347/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGGCCAGAAAAAAAGACAAGCACATTTTTTTTTTCCTCCTCTGTTTTTAg > 1:1309040/1‑90 (MQ=255) acaacaGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGaaa < 1:580721/90‑1 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGC‑GGGGCCCGAAAAAAAGACAAGCACAATTTTTTTTTCCTCATCTGTTTTTAAATGAGaaaaaa > 2:201032/1‑87 (MQ=255) aGTAAATAATAAAAATTTCTCGAGGC‑GGGGCCCGAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTGTTTTTAAATGGGAAACAAGa > 1:89510/1‑90 (MQ=255) ataaAAATTTCTCGAGGC‑GGGGTCAGGAAAAGAAACAAAGTCAATTTTTTTTTTCCCATCTTTTTTTAGAAGAGAAAAAAGAAATTTaaa > 2:714038/1‑89 (MQ=255) ataaAAATTTCTCGAGGC‑GGGGCCGGAAAAAACGGAATGCACATTTTATTTTTCCTCCTTTATTTTTAAATGTGAAACACAAAATTTAAc > 1:1061231/1‑90 (MQ=255) aaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGc < 2:552760/90‑1 (MQ=255) aaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCt > 1:87290/1‑90 (MQ=255) ttttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAAtt > 2:1454932/1‑90 (MQ=255) attttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACt < 2:698160/90‑1 (MQ=255) | GATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACT > NZ_CP009273/2380275‑2380446 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |