| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NZ_CP009273 | 1,560,820 | T→C | E193G (GAA→GGA) | dosC ← | diguanylate cyclase DosC |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,560,820 | 0 | T | C | 90.0% | 47.0 / ‑0.4 | 20 | E193G (GAA→GGA) | dosC | diguanylate cyclase DosC |
| Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (7/11); total (9/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.89e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01 | |||||||||||
TAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCTT > NZ_CP009273/1560734‑1560904 | tAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTcccc < 2:954849/90‑1 (MQ=255) acaGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGaaa < 1:551535/90‑1 (MQ=255) cAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTg < 2:616087/90‑1 (MQ=255) ttAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGtt < 1:185912/90‑1 (MQ=255) aaaGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCt > 1:242439/1‑90 (MQ=255) ggCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGACGTTCtttt < 2:85423/90‑1 (MQ=255) cAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTttc > 2:436649/1‑90 (MQ=255) aGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAg < 2:865429/90‑1 (MQ=255) gAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGt > 1:76375/1‑90 (MQ=255) gAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGt > 2:1462830/1‑90 (MQ=255) aTCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTaa > 2:625427/1‑90 (MQ=255) cagGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGagaa < 1:563454/90‑1 (MQ=255) cagGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGagaa < 1:945316/90‑1 (MQ=255) gATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGtt > 2:292427/1‑90 (MQ=255) gATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGagaa < 2:1292671/78‑1 (MQ=255) gATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGagaa > 1:1292671/1‑78 (MQ=255) atatCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCa > 1:1346528/1‑90 (MQ=255) tCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCt > 2:348684/1‑90 (MQ=255) tCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCt < 1:259169/90‑1 (MQ=255) cTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCtt < 2:619770/90‑1 (MQ=255) cTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCtt < 1:348684/90‑1 (MQ=255) | TAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCTT > NZ_CP009273/1560734‑1560904 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |