Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,380,341 | 0 | G | A | 55.0% | ‑3.3 / 22.5 | 20 | intergenic (+31/‑72) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
Reads supporting (aligned to +/- strand): ref base G (4/5); new base A (11/0); total (15/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.13e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACA‑TGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCT > NZ_CP009273/2380273‑2380428 | tGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAGGGAA‑GGGACatttttttttt > 2:326244/1‑89 (MQ=255) ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGAAA‑AGCAAATTTTTTTTTTCCttc > 2:333501/1‑88 (MQ=255) ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCAAAAAAAAGACA‑AGCCCAATTTTTTTTTCCTcc > 2:320954/1‑89 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACA‑TGCACATTTTATTTTCCCTCATc < 2:88746/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACA‑TGCACATTTTATTTTCCCTCATc < 1:119819/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACA‑TGCACATTTTATTTTCCCTCATc < 1:519654/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACA‑TGCACATTTTATTTTCCCTCATc < 1:495638/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACA‑TGCACATTTTATTTTCCCTCATc < 1:43928/90‑1 (MQ=255) ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAAAAA‑AACACAATTTATTTTTTCTCTTCTg > 1:15813/1‑90 (MQ=255) ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGCCCGGGAAAAACGAAA‑TGCACATTTTTTTTTTCCTCATCtt > 2:519625/1‑89 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAAAAAAAAAGACA‑AGCACATTTTTTTTTTCCTCCTCTGtttttaa > 2:40561/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCCGAAAAAAAGACA‑ATCACATTTTATTTTTCCTCATCtttttttaa > 1:478353/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCAAAAAAAAGAAA‑AGCACATTTTTTTTTTTCTCCTCTGtttttaa > 2:493230/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACA‑AGCACATTTTTTTTTTTCTCCTCtttttttaa > 2:343941/1‑87 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAAACA‑TGCAAAATTTTTTTTTTCTCATCTGTTTTTAAATGGGaaaaaa > 1:63899/1‑87 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACA‑TGCAAATTTTTTTTTTTCTCCTTTGTTTTTAGATGTGGAAAATGaaaa > 2:168059/1‑89 (MQ=255) aaaaTTTCTCGAGGCGGGACCAGAAAAAGACACACTGAAACTTTTTTTTTCCC‑CATCTGTTTTTAGAAGCGAAACATTAAAATCAAAgga > 2:398576/1‑90 (MQ=255) gAAAAAGCGACA‑TGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGAtt > 1:169158/1‑90 (MQ=255) aaGCGACA‑TGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAGACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCt > 2:569324/1‑90 (MQ=255) aaGCGACA‑TGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCt > 1:541737/1‑90 (MQ=255) | TGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACA‑TGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCT > NZ_CP009273/2380273‑2380428 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |