Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,380,341 | 0 | G | A | 66.7% | 6.3 / 14.2 | 21 | intergenic (+31/‑72) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
Reads supporting (aligned to +/- strand): ref base G (1/6); new base A (14/0); total (15/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.29e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑A‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGA > NZ_CP009273/2380260‑2380418 | aTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGc < 1:1602187/90‑1 (MQ=255) gCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGcac‑a < 2:1518435/90‑1 (MQ=255) tGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑A‑TTTTATTTTc < 1:1017782/90‑1 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCAA‑A‑ATTTTTTTTTCCTCAt > 2:989389/1‑90 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑A‑TTTTATTTTCCCTCATc < 1:1422190/90‑1 (MQ=255) ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCAC‑A‑TTTTTTTTTTCCTCATCtttt > 2:1246247/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACATACACAA‑TTTTTTTTTCCCCCTTTTTTTTTaa > 2:356477/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCAC‑A‑TTTTTTTTTTCCTCATCtttttttaa > 2:1057100/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAAAAAAAACGACAAGCAC‑A‑TTTTTTTTTTCCTCCTCTGtttttaa > 2:1602514/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCCGAAAAAAAGAAAAGCAC‑A‑ATTTTTTTTTTCTCATCTGttttttt > 1:1517736/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAAGAAAAACAC‑A‑ATTTTTTTTTTCTCCTCtttttttaa > 1:1510369/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGAAAAGCAC‑A‑TTTTTTTTTTTCTCATCtttttttaa > 2:285053/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGCAC‑A‑TTTTTTTTTTCCTCCTCtttttttaa > 2:334653/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGAAC‑A‑TTTTTTTTTTCCTCATCTGtttttaa > 2:1323797/1‑87 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCAAAAAAAAAGACAAGCAC‑A‑TTTTTTTTTTCCCCCTCTTTTTTTAAATGCGaaaaaa > 2:1223496/1‑87 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGAAAAGAAC‑A‑TTTTTTTTTTTCTCCTCTGTTTTTAGAAGGGAAACAAGaaaa > 1:957771/1‑89 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAAACACAGAC‑ATTTTTTTTTTTCCTC‑TTTGTTTTTAAATAGGAAAGATGaaaa > 1:1564917/1‑89 (MQ=255) aaaaTTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCAC‑A‑TTTTTTTTTTCCTCCTCTTTTTTTAAATGGGAAAAAAGAAAATCAACgga > 2:418572/1‑90 (MQ=255) aaaaTTTCTCGAGGCGGGGCCAGAAAAAGCCACATTTCC‑A‑TTTTATTTTTCCTCGTCTGTTTTAATAAGGGGAAAATCAAATTCAACcga > 2:1584276/1‑90 (MQ=255) aGGCGGGTCCAGAAAAAGCGACATGCAC‑A‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAaaataaa < 1:704870/90‑1 (MQ=255) ggTCCAGAAAAAGCGACATGCAC‑A‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGa < 2:685036/90‑1 (MQ=255) | ATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑A‑TTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGA > NZ_CP009273/2380260‑2380418 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |