Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,380,362 | 0 | C | T | 53.8% | 10.6 / 56.4 | 39 | intergenic (+52/‑51) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
Reads supporting (aligned to +/- strand): ref base C (3/15); new base T (21/0); total (24/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.25e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.81e-09 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAA > NZ_CP009273/2380279‑2380451 | gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCAGAAAAAAAGACAAACACAATTTTTTTTTTCTCAt > 2:1584864/1‑90 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACATGCACATTTTTTTTTTCCTcct > 1:670656/1‑88 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 1:1060488/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 2:37078/90‑1 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc < 2:1507131/90‑1 (MQ=255) ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCTGGAAAAAAGAAATGCACATTTTTTTTTTTCTCATCtt > 2:750746/1‑89 (MQ=255) ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACATGCACCCTTTTTTTTTTCTCctctt > 1:778196/1‑89 (MQ=255) cTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTaa < 1:1584864/90‑1 (MQ=255) cTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAAGAAAAACACAATTTTTTTTTTCTCATCtttttttt > 2:352976/1‑88 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAg < 1:1650453/90‑1 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGGCAAGAACCTTTTTTTTTTCCTCTTCTCttttttg > 1:1313346/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGAAAAAAACAATTTTTTTTTCCTCCTCTGtttttaa > 2:246077/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAAAAAAAAAGAAAAAAACATTTTTTTTTTTCTCCTCTTTTTTTAg > 2:689045/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAAGAAAAACAAATTTTTTTTTTTCTCTTCtttttttaa > 2:593597/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGCACAATTTTTTTTTCCTCTTCtttttttta > 2:1060488/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCCTCTGTTTTTAg > 1:1191226/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAAAAGAACAATTTTTTTTTTCTCATCtttttttaa > 2:1542684/1‑87 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCAAAAAAAAAACAAGCACAATTTTTTTTTCCTCCTCtttttttaaa > 1:662329/1‑86 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAAAAAAAAAGAAAAGAACAATTTTTTTTTTTTTATCTGttttttga > 2:1168071/1‑86 (MQ=255) ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCTGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGa < 1:750746/90‑1 (MQ=255) ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCTCATTTTATTTTCCCTCATCTGTTTTAAGATGCGa < 1:569437/90‑1 (MQ=255) acaacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCTACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGaaa < 1:273856/90‑1 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACATGCACATTTTTTTTTTCCTCATCTGTTTTTAAATGCGaaaaaa > 1:1448437/1‑87 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACATTTTTTTTTTTCTCCTTTTTTTTTAAATGCGaaaaaa > 1:799542/1‑87 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAACACATTTTTTTTTTTCTCCTCCTTTTTTAAATGGGaaaaaa > 2:622352/1‑87 (MQ=255) aGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGa < 1:540878/90‑1 (MQ=255) aGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGa < 1:1020455/90‑1 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAAt < 1:743566/90‑1 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGGCCCAAAAAAAAGAAATGAACAGTTTATTTTTCCTCCTCTTTTTTTAAATGCGAAACATGaaaa > 1:366562/1‑89 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACATGTACAATTTTTTTTTTCTCCTCTTTTTTTAAATGCGGAAAATaaaaa > 1:1718502/1‑89 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAAAAGCACAATTTTTTTTTCCTCCTCTTTTTTTAAAAGGGGAAAATGaaaa > 1:1102002/1‑89 (MQ=255) ataaAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACAATTTTTTTTTCCTCCTTTGTTTTTAGATGGGAAAAAAGAAAATTAAc > 1:40163/1‑90 (MQ=255) taaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACg < 2:1448437/90‑1 (MQ=255) tCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAt < 1:1370657/90‑1 (MQ=255) aGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTa > 1:30264/1‑90 (MQ=255) cGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTAtt < 2:352026/90‑1 (MQ=255) aTGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGc < 2:40163/90‑1 (MQ=255) tCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa > 2:199112/1‑90 (MQ=255) tCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa > 1:368270/1‑90 (MQ=255) cccTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGaa < 1:420853/90‑1 (MQ=255) | GTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAA > NZ_CP009273/2380279‑2380451 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |