Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 1,401,680 | G→A | R125C (CGC→TGC) | dgcM ← | diguanylate cyclase DgcM |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,401,680 | 0 | G | A | 94.6% | 192.1 / 1.1 | 56 | R125C (CGC→TGC) | dgcM | diguanylate cyclase DgcM |
Reads supporting (aligned to +/- strand): ref base G (1/2); new base A (30/23); total (31/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.81e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01 |
TTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCT > NZ_CP009273/1401553‑1401813 | ttATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTg > 1:392987/1‑139 (MQ=255) aGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACt > 2:93856/1‑139 (MQ=255) aGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACt > 1:769189/1‑139 (MQ=255) aGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACt > 2:1218749/1‑139 (MQ=255) gAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACt < 1:1187046/138‑1 (MQ=255) aTTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCt > 1:321459/1‑139 (MQ=255) aGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCtttt < 2:392987/138‑1 (MQ=255) aGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCtttt > 2:1281996/1‑138 (MQ=255) acgaTTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACt < 1:1189544/110‑1 (MQ=255) acgaTTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACt > 2:1189544/1‑110 (MQ=255) acgaTTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACt > 1:611312/1‑110 (MQ=255) acgaTTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACt < 2:611312/110‑1 (MQ=255) acgaTTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAAAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCtt > 2:310369/1‑139 (MQ=255) tCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATAcc < 2:772452/138‑1 (MQ=255) gTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGa < 2:1286032/135‑1 (MQ=255) gACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAaataa > 1:725694/1‑139 (MQ=255) ggTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCAc < 2:769189/139‑1 (MQ=255) ggTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCAc < 1:930319/139‑1 (MQ=255) gTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCtt < 1:27584/125‑1 (MQ=255) gTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCtt > 2:27584/1‑125 (MQ=255) aatcaaCAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGt > 1:261504/1‑139 (MQ=255) atcaaCAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGtt > 1:328020/1‑139 (MQ=255) tcaaCAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATa > 2:815733/1‑116 (MQ=255) tcaaCAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATa < 1:815733/116‑1 (MQ=255) acaacaTTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGttctgt < 1:1218749/139‑1 (MQ=255) acaTTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGttctgtttc > 1:1257779/1‑139 (MQ=255) tGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGttctgtttctgt > 2:196907/1‑138 (MQ=255) gTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTG‑TTTTTGCACTTGTAATTTGCCGAGCGACTGGCTTTTAAACCCAACGTTGCCGGCGGTTCGATAACTTCGGAAATAATCACCGGGTCTGGTTCTGTCAGc > 1:1279709/1‑139 (MQ=255) agTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCg > 2:177303/1‑139 (MQ=255) agTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCg > 2:603912/1‑139 (MQ=255) aGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGc > 1:544576/1‑56 (MQ=255) aGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGc < 2:544576/56‑1 (MQ=255) aGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCaa > 1:1046133/1‑138 (MQ=255) gTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAaata < 2:364067/105‑1 (MQ=255) gTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAaata > 1:364067/1‑105 (MQ=255) gTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCaaaa > 1:1138442/1‑139 (MQ=255) gTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCaaaa > 1:125214/1‑139 (MQ=255) gTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCaaaa > 1:781196/1‑139 (MQ=255) tCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAAc < 2:1046133/139‑1 (MQ=255) aaaaaaGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCgg < 1:631396/78‑1 (MQ=255) aaaaaaGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCAACTGGCTTTTAAACCCAGCGTTGCCgg > 2:631396/1‑78 (MQ=255) aaaGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGa < 2:328020/139‑1 (MQ=255) aaaGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGa < 1:93856/139‑1 (MQ=255) aaaGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGa < 2:1243586/139‑1 (MQ=255) aaaGCGCGCATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGa < 2:1174286/139‑1 (MQ=255) gcATAAAAACCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTc > 2:567700/1‑138 (MQ=255) cTTGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTc < 1:310369/139‑1 (MQ=255) tGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGgc > 2:408350/1‑57 (MQ=255) tGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGgc < 1:408350/57‑1 (MQ=255) tGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCAc < 1:661434/83‑1 (MQ=255) tGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCAc > 2:661434/1‑83 (MQ=255) tGTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCct > 2:83338/1‑139 (MQ=255) gTTTTTTGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCct < 1:1281996/138‑1 (MQ=255) ttGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCtcttct < 1:226088/138‑1 (MQ=255) tGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATATCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCt > 1:827323/1‑125 (MQ=255) tGCACTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATATCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCt < 2:827323/125‑1 (MQ=255) | TTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCT > NZ_CP009273/1401553‑1401813 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |