| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,152,727 | 0 | T | G | 54.0% | 36.1 / 55.2 | 37 | G165G (GGT→GGG) | ycfH | metal‑dependent hydrolase |
| Reads supporting (aligned to +/- strand): ref base T (1/16); new base G (20/0); total (21/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.32e-09 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-08 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTC‑GGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGT > NZ_CP009273/1152640‑1152816 | cGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTaa < 1:205605/90‑1 (MQ=255) tACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGaaaat > 1:297569/1‑88 (MQ=255) cTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTaa < 1:712218/83‑1 (MQ=255) cTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGaa > 2:712218/1‑83 (MQ=255) tGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTgg < 2:416687/90‑1 (MQ=255) aaTTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGAtct < 1:760907/90‑1 (MQ=255) aTTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGAtctc < 1:224100/90‑1 (MQ=255) tGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATTGATCTT‑gggt > 2:838211/1‑88 (MQ=255) tGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATGGGTCTC‑gggt > 1:473298/1‑88 (MQ=255) cgcgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTGGATCTC‑GGGttt > 1:378210/1‑90 (MQ=255) aaGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTC‑GGATTTTACATCT‑CCt < 2:571105/90‑1 (MQ=255) gACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATTGATCTT‑GGGTTTTACATCT‑CCttttt > 1:1636842/1‑89 (MQ=255) cGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTC‑GGATTTTACATCT‑CCTTTTCCg < 1:676145/90‑1 (MQ=255) cGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTC‑GGATTTTACATCT‑CCTTTTCCg < 1:959778/90‑1 (MQ=255) ggATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTC‑GGATGTTACATCT‑CCTTTTCCgg < 2:726700/90‑1 (MQ=255) gTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTC‑GGATTTTACATCT‑CCTTTTCCGGCATTGTGa < 1:594658/90‑1 (MQ=255) ggCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGGAAATATTGGATCTC‑GGGTTTTTAATCT‑CCTTTTTTGGCAGTGTGTcc > 2:1361571/1‑90 (MQ=255) ggCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGATCTC‑GGGTTTTTCATCT‑CCCTTTTCCGCATTGTGGcc > 1:324555/1‑90 (MQ=255) ggCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATTATTGAGATC‑GGGTTTTTTATCT‑CCTTTTCCTGCCAtgtgggc > 2:2249567/1‑87 (MQ=255) gCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTC‑GGATTTTACATCT‑CCTTTTCCGGCATTGTGACCt < 1:122083/90‑1 (MQ=255) tacACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATGGGTCTCTGGGTTTTTAATTT‑CTTTTTCC‑GCGTTGTGACCTTcccga > 2:2200338/1‑87 (MQ=255) acacTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATTGATCTC‑GGGTTTTACATCTCCCTTTTCCGCCATT‑TGACCTTCCCTaa > 2:550781/1‑90 (MQ=255) acTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTGGATCTC‑GGGTTTTATATCT‑CCTTTTTCTGCATTGTTTCCTTTCGTAATg > 2:906199/1‑90 (MQ=255) cTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATGGATCTC‑GGGTTTTTAATTT‑CCCTTTTCGGCATTGTGGCCTTCCCTAATgg > 2:884715/1‑89 (MQ=255) gagGACAGAGAAACGGCGGGGAAAATAATTGGTCTC‑GGGTTTTACATCT‑CCTTTTTCCGCATTGTGGCCTTCCCTAATGCGGGGCAAATg > 1:683222/1‑90 (MQ=255) ggACAGAGAAACGGCGGGGAAATTATTGGGTTTT‑GGGTTTTTAATTT‑CCTTTTTCGGGATTGTGGGCCTCCCTAATGCGGAGGAAATgcg > 2:963471/1‑90 (MQ=255) ggACAGAGAAACGGCGGGGAAAATTCTGGGTCTC‑GGGTTTTTCATCT‑CCTTTTTCCGCATTTTGACCTTCCCTAATTCGGGGGAACTgcg > 2:931810/1‑90 (MQ=255) ggACAGAGAAACGGCGGGGAAAATACTTGGTCTC‑GGGTTTTACATCT‑CCTTTTTCCCCCTTTTGACCCTCCCTAATGCGGAGCAACcgcg > 2:968138/1‑90 (MQ=255) ggACAGAGAAACGGCGGGGAAAATAATTGGGCTC‑GGGTTTTATATTTCCCTTTTCCGGC‑TTTTGGCCTTCCGTTATTCGGGGGAACAgcg > 2:1117641/1‑90 (MQ=255) gACAGAGAAACGGCGGGTAAATTACTGGATCTC‑GGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTgcgc < 2:2247817/90‑1 (MQ=255) gACAGAGAAACGGCGGGTAAATTACTGGATCTC‑GGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTgcgc < 2:2208708/90‑1 (MQ=255) gACAGAGAAACGGCGGGTAAATTACTGGATCTC‑GGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTgcgc < 2:1433142/90‑1 (MQ=255) gagaAACGGCGGGTAAATTATTTGGTCTC‑GGGTTTTAAATCT‑CCTTTTTCCGCATTGTGACCTTCCGTAATTCGGGGGAACTGCGCGAtg > 1:1752677/1‑90 (MQ=255) gagaAACGGCGGGGAAATTAATGGATTTC‑GGAGTTTAAATCT‑CCTTTTTCGGCATTTTGACCTTCCGTAATTCGGAGCAACCGCGCGAtg > 1:2169286/1‑90 (MQ=255) gagaAACGGCGGGGAAAATAATGGGTTTC‑GGGTTTTTTATCT‑CCTTTTTCCGCGTTGTGGCCCTCCCTAATGCGGGGGAACCGCGCGagg > 2:1369475/1‑88 (MQ=255) aaaCGGCGGGTAAATTACTGGATCTC‑GGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGAtgctg < 2:1563740/90‑1 (MQ=255) aaCGGCGGGTAAATTACTGGATCTC‑GGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGAtgctgc < 1:377043/90‑1 (MQ=255) tAAATTACTGGATCTC‑GGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGt < 1:906199/90‑1 (MQ=255) tAAATTACTGGATCTC‑GGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGt < 1:2200338/90‑1 (MQ=255) tAAATTACTGGATCTC‑GGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGt < 1:1468963/90‑1 (MQ=255) | CGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTC‑GGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGT > NZ_CP009273/1152640‑1152816 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |