Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,380,341 | 0 | G | A | 45.2% | 20.9 / 22.1 | 32 | intergenic (+31/‑72) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
Reads supporting (aligned to +/- strand): ref base G (2/15); new base A (14/0); total (16/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.09e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCA‑C‑ATTTTATTTTCCCTCAT‑CTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGC > NZ_CP009273/2380253‑2380427 | tCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGc < 2:679895/90‑1 (MQ=255) aTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGc < 1:1036477/90‑1 (MQ=255) cAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCA‑C‑attt < 2:1217730/90‑1 (MQ=255) cAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCA‑C‑attt < 1:823553/90‑1 (MQ=255) tGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACC‑attttttttt > 2:1028196/1‑90 (MQ=255) gATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCA‑C‑ATTTTATTTTccc < 2:1381197/90‑1 (MQ=255) gATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCA‑C‑ATTTTATTTTccc < 2:1306150/90‑1 (MQ=255) aTGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAAGAAAAGCA‑C‑atttttttttttct > 2:1476052/1‑86 (MQ=255) gtgtTAT‑ACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCA‑C‑ATTTTATTTTCCCTCAT‑c < 1:24261/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCA‑C‑ATTTTATTTTCCCTCAT‑c < 1:1476052/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCA‑C‑ATTTTATTTTCCCTCAT‑c < 1:447822/90‑1 (MQ=255) ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGAA‑C‑AATTTTTTTTTTCTCct‑ctt > 1:986951/1‑89 (MQ=255) ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGACATGCA‑C‑ATTTTTTTTTTCCTCCT‑Ctttt > 2:1453011/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGAAAAGCA‑C‑AATTTTTTTTTTCTCCT‑CTCtttttat > 2:970335/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACATGCA‑C‑ATTTTTTTTTTCCTCCT‑Ctttttttaa > 2:447822/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGGA‑CAATTTTTTTTTTTCTCCT‑Cttttttta > 1:383232/1‑88 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCAAAAAAACGACAAGGA‑CAATTTTATTTT‑TCTCATCCTGTTTTac > 1:1091422/1‑89 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCATAAAAAAAGGAAAGCA‑A‑AATTTTTTTTCTTTCAT‑CTGTTTTAAg > 1:829155/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAAAAAAGGA‑C‑ATTTTTTTTTTCCTCTT‑CTGtttttaa > 1:1111899/1‑87 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGAAACGCA‑C‑ATTTTATTTTTCCTCGT‑CTGtttttata > 2:1087909/1‑86 (MQ=255) aGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGGAAAAAAGAAAAAGA‑C‑ATTTTTTTTTTTCTTTT‑CCGTTTTTAGATGAGAAACaaaa > 1:831981/1‑87 (MQ=255) ataaAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCA‑C‑ATTTTATTTTTCCTCCT‑CTGTTTTTAAATGCGAAAAAAGAAAATCaaa > 2:87718/1‑89 (MQ=255) aaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCA‑C‑ATTTTATTTTCCCTCAT‑CTGTTTTAAGATGCGAAACATGAAATTCAACgga < 1:412501/90‑1 (MQ=255) aaaaTTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCA‑C‑CTTTTTTTTTTCCTCCT‑CTGTTTTTAAATGGGGAAAAAGAAAATCaaaagg > 2:323055/1‑86 (MQ=255) tcGAGGCGGGTCCAGAAAAAGCGACATGCA‑C‑ATTTTATTTTCCCTCAT‑CTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAaaat < 2:619052/90‑1 (MQ=255) ggggCCCGAAAAAAAGACATGCA‑C‑ATTTTTTTTTTCCTCCT‑CTGTTTTTAGATGGGAAAAAAGAAAATTAAAGGAGGGAAAAAAAATtgg > 1:1292904/1‑90 (MQ=255) gggTCCAGAAAAAGCGACATGCC‑C‑ATTTTATTTTCCCTCAT‑CTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATggg < 2:352407/90‑1 (MQ=255) gTCCAGAAAAAGCGACATGCA‑C‑ATTTTATTTTCCCTCAT‑CTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGaa < 1:323055/90‑1 (MQ=255) tCCAGAAAAAGCGACATGCA‑C‑ATTTTATTTTCCCTCAT‑CTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAt < 2:192994/90‑1 (MQ=255) aaaTCGACATGCA‑C‑ATTTTATTTTCCCTCAT‑CTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGc < 1:67812/90‑1 (MQ=255) aGCGACATGCA‑C‑ATTTTATTTTCCCTCAT‑CTGTTTTAAGATGCGAAACATGAAATTCAACggagga < 2:1015597/66‑1 (MQ=255) aGCGACATGCA‑C‑ATTTTATTTTCCCTCAT‑CTGTTTTAAGATGCGAAACATGAAATTCAACggagga > 1:1015597/1‑66 (MQ=255) | TCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCA‑C‑ATTTTATTTTCCCTCAT‑CTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGC > NZ_CP009273/2380253‑2380427 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |