Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,290,243 | 0 | T | G | 58.3% | 29.7 / 43.5 | 37 | L87F (TTA→TTC) | ccmB | heme exporter protein CcmB |
Reads supporting (aligned to +/- strand): ref base T (5/10); new base G (21/0); total (27/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.18e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCTACCAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑CAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAAGCAGC > NZ_CP009273/2290157‑2290330 | gCTACCAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑GAAGGGGAAc > 1:888012/1‑90 (MQ=255) gCTACCAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑GAAGGGGAAc > 1:26831/1‑90 (MQ=255) gCTACCAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑CAAGGGGAAc > 2:990498/1‑90 (MQ=255) cAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑GAAGGGGAACgggca > 1:898750/1‑87 (MQ=255) cAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑CAAGGGTAACGGcaa < 2:1681441/90‑1 (MQ=255) aGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑CAAGGGTAACGGcaac > 2:1573494/1‑90 (MQ=255) gTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑GAAGGGGAAAGGcaaca > 1:1658963/1‑90 (MQ=255) gTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑CAAGGGTAACGGcaaca < 2:1501297/90‑1 (MQ=255) tGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑GAAGGGGGAAGGcaacaa > 1:566165/1‑90 (MQ=255) gCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑CAAGGGGAACGGCAacaaca > 2:751974/1‑90 (MQ=255) cGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑CCAGGGGAACCGCAaaaacaa > 2:1093537/1‑89 (MQ=255) gAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑CAAGGGGAAAGGcaacaaaaac > 1:144135/1‑86 (MQ=255) aaaGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑CAAGGGTAACGGCAACAACATCa < 1:88956/90‑1 (MQ=255) aaaGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑CAAGGGTAACGGCAACAACATCa < 2:849719/90‑1 (MQ=255) aGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑CAAGGGTAACGGCAACAACATCAAt < 1:1634267/90‑1 (MQ=255) gAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑CAAGGGGAACGGCAACAACATCAATTGTTc > 2:355764/1‑90 (MQ=255) aCGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑CAAGGGTAACGGCAACAACATCAATTGTTCAAGAc < 1:1484500/90‑1 (MQ=255) gCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑CAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGc < 1:563221/90‑1 (MQ=255) aGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑GAAGGGGAAAGGGAAAAACATCAATTGTTCAAAACTTccc > 1:873614/1‑89 (MQ=255) ccGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑CAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTcc < 1:1573494/90‑1 (MQ=255) aGTGCGCCATCACCTTCGCCAGCACAACGGCGGGACAAGGG‑GAGGGAAAGAAAAACAATTGTTTAAGACTGCCCTCCCGTAAACCGTCAc > 2:1627836/1‑90 (MQ=255) aGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑CAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCAc > 2:206773/1‑90 (MQ=255) aGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑CAAGGGGAACGGCAACAACATCAATTTTTCAAAACTGCCGGCCTGCAAAACGTCAc > 2:1830620/1‑90 (MQ=255) aGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑CAAGGGGAACGGCAACAACATCAATTGTTTAAGACTGCCGGCCTGCAAATCGTCAc > 1:1029934/1‑90 (MQ=255) aGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑CAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCAc > 2:1866223/1‑90 (MQ=255) gcgcCATCACCTTCGCCAGCACAACGGCGGG‑GAAGGGGAACGGCAACAAAATCAAATTTTTAAGACTGCCGGCCTTCAAATCGTCACGGa > 1:940183/1‑90 (MQ=255) aTCACCTTCGCCAGCACAACGGCGGG‑CAAGGGGAACGGCAACAACATCAATTGTTCAAAACTGCCGTCCCGCAAATCGTCACGGAACAGt > 1:1624413/1‑90 (MQ=255) ttCGCCAGCACAACGGCGGG‑CAAGGGTAACGGcaaca < 2:1522715/37‑1 (MQ=255) ttCGCCAGCACAACGGCGGG‑CAAGGGGAACGGcaaca > 1:1522715/1‑37 (MQ=255) tCGCCAGCACAACGGCGGG‑CAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCa > 2:496040/1‑90 (MQ=255) aGCACAACGGCGGG‑CAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGcc > 2:226020/1‑90 (MQ=255) gCACAACGGCGGG‑CAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCa > 2:671058/1‑90 (MQ=255) aaCGGCGGG‑CAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCaa < 2:873614/90‑1 (MQ=255) ggg‑CAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGa > 1:727418/1‑90 (MQ=255) ggg‑CAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGa > 2:1787572/1‑90 (MQ=255) ggg‑CAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGa > 1:1702414/1‑90 (MQ=255) ggTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAagcagc > 2:1697166/1‑90 (MQ=255) | GCTACCAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGG‑CAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAAGCAGC > NZ_CP009273/2290157‑2290330 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |