| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NZ_CP009273 | 1,401,542 | T→C | M171V (ATG→GTG) | dgcM ← | diguanylate cyclase DgcM |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,401,542 | 0 | T | C | 88.9% | 67.4 / 3.7 | 27 | M171V (ATG→GTG) | dgcM | diguanylate cyclase DgcM |
| Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (10/14); total (12/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.69e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01 | |||||||||||
GACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGC > NZ_CP009273/1401407‑1401679 | gACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGt > 1:811474/1‑139 (MQ=255)gACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACACCGt > 1:516241/1‑139 (MQ=255) gAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGGGCATATTTATTTCCCAGGTATGTTTCTGGCACACCGtctc < 2:375100/139‑1 (MQ=255) gAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTc < 2:811474/139‑1 (MQ=255) ccATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACACCGTTTCATGGTtat < 2:516241/139‑1 (MQ=255) tGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACACCGTTTCATGGTTATAACCATAGAAAtt > 2:733417/1‑137 (MQ=255) aCAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACACCGTTTCATGGTTATAACCATAGAAATTGAgcg < 2:188002/139‑1 (MQ=255) ggTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACACCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAg < 1:90256/139‑1 (MQ=255) tATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACACCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATc > 2:410711/1‑139 (MQ=255) aaTGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACACCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTc < 1:184090/135‑1 (MQ=255) aaTGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACACCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTc > 2:184090/1‑135 (MQ=255) agaTTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACACCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGa > 1:618228/1‑139 (MQ=255) catgcatgATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACACCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGt < 2:858968/138‑1 (MQ=255) catgcatgATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACACCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGt < 1:271293/138‑1 (MQ=255) gATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACACCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATcaac < 2:546694/139‑1 (MQ=255) aTAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACACCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATcaac < 2:602580/138‑1 (MQ=255) aGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACACCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATcaacaa > 2:855832/1‑138 (MQ=255) gCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACACCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACAt < 1:733417/139‑1 (MQ=255) acgcgacgcCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACACCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTg < 1:89266/139‑1 (MQ=255) acgcgacgcCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACACCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTg > 2:40452/1‑139 (MQ=255) acgcCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACACCGTTTCATGGTTATAACCATAGAAATTGAGCgcggc > 2:260130/1‑78 (MQ=255) acgcCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACACCGTTTCATGGTTATAACCATAGAAATTGAGCgcggc < 1:260130/78‑1 (MQ=255) aGGTATGTTTCTGGCACACCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAacaaca < 2:491920/105‑1 (MQ=255) aGGTATGTTTCTGGCACACCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAacaaca > 1:491920/1‑105 (MQ=255) tGTTTCTGGCACACCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGtcaatcaa > 1:613198/1‑95 (MQ=255) tGTTTCTGGCACACCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGtcaatcaa < 2:613198/95‑1 (MQ=255) aTCGTTTCATGGTTATAACCATAGAAATTGAGAGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTgc > 1:740396/1‑139 (MQ=255) | GACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGC > NZ_CP009273/1401407‑1401679 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |