Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 296,422 | 0 | G | C | 22.2% | 170.3 / 13.3 | 63 | P75R (CCG→CGG) | paoC | aldehyde oxidoreductase molybdenum‑binding subunit PaoC |
Reads supporting (aligned to +/- strand): ref base G (13/36); new base C (14/0); total (27/36) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.37e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGG > NZ_CP009273/296335‑296511 | gTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTGATGAAAGCCAGGGC‑GCCCCgc > 2:373127/1‑90 (MQ=255) gTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAACCCGTAC‑GCCCCgc > 1:63175/1‑90 (MQ=255) gTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGAAAACCCGGAC‑GCCCCgc > 1:34455/1‑90 (MQ=255) ggggCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACACCCCGTAC‑GCCCCgcgc > 1:410783/1‑90 (MQ=255) ggggCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAAGACAGCCCGTAC‑GCCCCgcgc > 2:41554/1‑90 (MQ=255) gccgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCCGCGCttt > 1:142596/1‑90 (MQ=255) gccgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGCCCGGAC‑GCCCCGCGCctt > 1:435246/1‑90 (MQ=255) cgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCCGTAC‑GCCCGGCGCttttt > 1:142431/1‑90 (MQ=255) cgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCttttt > 1:15777/1‑90 (MQ=255) cgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGGATGGACGCCACTACAGCCCGGGGTtttt > 2:231980/1‑90 (MQ=255) cgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGCCACAAC‑CCCCCGGGCttttt > 2:173975/1‑90 (MQ=255) cgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGAAAGCCAGTAC‑GCCCGGCGCttttt > 1:50159/1‑90 (MQ=255) gccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAAGAAAGCCCGTAC‑GCCCCGCGCcttttt > 2:76133/1‑89 (MQ=255) cAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCgcggc < 1:439328/90‑1 (MQ=255) cTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTcc < 1:173975/90‑1 (MQ=255) cTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTcc < 2:15777/90‑1 (MQ=255) cTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTcc < 1:41554/90‑1 (MQ=255) tGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCCGTAC‑GCCCCGCGCTTTTTTTGCGGGGGCCg > 1:351674/1‑90 (MQ=255) gtgtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCa < 1:47938/90‑1 (MQ=255) gtgtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCa < 1:298927/90‑1 (MQ=255) gtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGgc < 2:286277/64‑1 (MQ=255) gtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGgc > 1:286277/1‑64 (MQ=255) tttGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGc < 2:286641/90‑1 (MQ=255) tGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCgg < 2:217654/90‑1 (MQ=255) gacagCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCgg < 2:310148/86‑1 (MQ=255) tCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTg < 2:142431/90‑1 (MQ=255) cGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCCGTAC‑GCCCGGCGCCTTTTTCGCGGCGGCCGGATCAAggggggggg > 1:338562/1‑87 (MQ=255) cTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGc < 2:240628/90‑1 (MQ=255) cTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGc < 2:127358/90‑1 (MQ=255) cTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGc < 2:123618/90‑1 (MQ=255) cTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGc < 1:149457/90‑1 (MQ=255) tttGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCg < 1:91115/90‑1 (MQ=255) gTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCCGCGCTTTTTTTGCGGGGGCCGGATCAAGGGGGGGGGGGGGGCCttttg > 2:149457/1‑88 (MQ=255) gTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGGAC‑CCCCCGCGCTTTTTGCGCGGCGGCCGTATCAAGGGGGGGGGGGGGGCCttttg > 1:326079/1‑88 (MQ=255) gCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCa < 1:140894/90‑1 (MQ=255) gCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCa < 2:149492/90‑1 (MQ=255) cccGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAAt < 1:463312/90‑1 (MQ=255) cccGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAAt < 2:153717/90‑1 (MQ=255) ccGGCGTTACTGGCGGTAAATACAACCAGTAC‑CCCCCGCGCCTTTTTCGCGGCGGCCGCATCAAAGGGGGGGGGGGGGCCTTTTGCAATg > 2:490639/1‑90 (MQ=255) ggCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGc < 2:410783/90‑1 (MQ=255) gCGTTACTGGCGGTAATGACAGCCCGTAC‑GCCCCGCGCCTTTTTCGCGGCGGCCGTATCAAGGGGGGGGGGGCGGCCTTTTGCAATgggg > 2:85430/1‑88 (MQ=255) aCTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAAcc < 1:408942/90‑1 (MQ=255) ggCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcgac < 2:472125/90‑1 (MQ=255) ggCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcgac < 1:76133/90‑1 (MQ=255) ggCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcgac < 2:376003/90‑1 (MQ=255) ggCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcgac < 2:34455/90‑1 (MQ=255) ggCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcgac < 2:351674/90‑1 (MQ=255) cGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAg < 2:147904/52‑1 (MQ=255) cGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAg > 1:147904/1‑52 (MQ=255) tAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGatat < 1:490639/90‑1 (MQ=255) tAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGatat < 1:313742/90‑1 (MQ=255) tAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGatat < 2:332341/90‑1 (MQ=255) aaTGACAGCCAGTAC‑GCCCGCCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGTCGTCCTTTGGCAATGGCGGAACCGACGAtata < 2:204492/90‑1 (MQ=255) aCAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCa < 1:492055/90‑1 (MQ=255) aGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATa > 2:283457/1‑90 (MQ=255) gCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAg > 2:281256/1‑90 (MQ=255) aGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAggcg > 1:335083/1‑90 (MQ=255) gTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAggcgg > 2:495431/1‑90 (MQ=255) gTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAggcgg < 1:73106/90‑1 (MQ=255) aC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGgcggcg < 2:118711/90‑1 (MQ=255) c‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGt > 2:180857/1‑90 (MQ=255) c‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGt < 2:50159/90‑1 (MQ=255) gCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGtt > 2:183656/1‑90 (MQ=255) ggCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcg > 1:15801/1‑63 (MQ=255) ggCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcg < 2:15801/63‑1 (MQ=255) ggCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTgggg > 1:50049/1‑90 (MQ=255) | GTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAC‑GCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGG > NZ_CP009273/296335‑296511 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |