Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 2,124,319 | G→A | P14L (CCG→CTG) | wcaD ← | colanic acid polymerase WcaD |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,124,319 | 0 | G | A | 86.2% | 166.8 / 14.2 | 58 | P14L (CCG→CTG) | wcaD | colanic acid polymerase WcaD |
Reads supporting (aligned to +/- strand): ref base G (4/4); new base A (18/32); total (22/36) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.62e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGT > NZ_CP009273/2124186‑2124454 | ggCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATcagcagcag > 2:1039303/1‑138 (MQ=255) tcatcaATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAg > 1:599504/1‑139 (MQ=255) tcatcaATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAg > 1:1102036/1‑139 (MQ=255) tcatcaATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACAAGCAGATAAATCAGCAGCAGCAACAg > 2:767935/1‑139 (MQ=255) cTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGa < 1:64330/139‑1 (MQ=255) cTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGa > 1:242334/1‑139 (MQ=255) aTGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGa > 1:1202575/1‑136 (MQ=255) tGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGa < 2:622595/139‑1 (MQ=255) tGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGa < 1:767935/139‑1 (MQ=255) cGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAg < 2:951017/139‑1 (MQ=255) cGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAg < 2:937391/139‑1 (MQ=255) ttCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGtt > 2:1277578/1‑139 (MQ=255) tCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGtt > 1:1100791/1‑138 (MQ=255) aaaaCAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACa < 2:1103152/138‑1 (MQ=255) aaaaCAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACa < 2:1049068/138‑1 (MQ=255) aaaaCAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACAt < 2:1202575/139‑1 (MQ=255) aaaCAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATa < 1:751247/139‑1 (MQ=255) aaaCAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATa < 1:1039303/139‑1 (MQ=255) agTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACtcctcc < 1:780368/139‑1 (MQ=255) agTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACtcctcc < 1:764249/139‑1 (MQ=255) cAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCAt < 1:1227742/139‑1 (MQ=255) aGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCt < 1:1050705/139‑1 (MQ=255) aGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCt < 1:674287/139‑1 (MQ=255) ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtgttgt < 1:45298/139‑1 (MQ=255) ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtgttgt < 2:862479/139‑1 (MQ=255) ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtgttgt < 1:504034/139‑1 (MQ=255) gTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCAt > 1:1102109/1‑131 (MQ=255) gTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCAt < 2:1102109/131‑1 (MQ=255) gTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGt < 2:1262330/112‑1 (MQ=255) gTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGt > 1:1262330/1‑112 (MQ=255) gTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTc > 2:438169/1‑139 (MQ=255) gTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTc > 2:1009430/1‑139 (MQ=255) gTAAAAAAGTGAAGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGAAAGAAGTTGACATACTCCTCCAGCATCtgttgt > 1:831148/1‑138 (MQ=255) aaaaaGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtgttgt > 2:841510/1‑135 (MQ=255) aaaaaGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtgttgt < 1:841510/135‑1 (MQ=255) aaaaaGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCa < 1:152221/138‑1 (MQ=255) aaaGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTg < 2:242334/139‑1 (MQ=255) gACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACa > 2:531413/1‑107 (MQ=255) gACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACa < 1:531413/107‑1 (MQ=255) gACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtgtt > 2:67898/1‑126 (MQ=255) gACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtgtt < 1:67898/126‑1 (MQ=255) aCTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGt > 2:473389/1‑134 (MQ=255) aCTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGt < 1:473389/134‑1 (MQ=255) tCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGaa < 2:595063/138‑1 (MQ=255) cGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACt > 1:182291/1‑139 (MQ=255) cGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtgttgt < 2:130309/105‑1 (MQ=255) cGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtgttgt > 1:130309/1‑105 (MQ=255) ggCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATa < 2:413976/76‑1 (MQ=255) ggCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATa > 1:413976/1‑76 (MQ=255) ggCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCg < 2:1093540/138‑1 (MQ=255) ggCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCg < 2:797721/138‑1 (MQ=255) gCGATTTTGACGTTAACCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTg < 1:1009430/139‑1 (MQ=255) aaCCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAg > 2:337204/1‑53 (MQ=255) aaCCAGCAGATAAATCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAg < 1:337204/53‑1 (MQ=255) aaaTCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGTCATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTg > 2:247130/1‑114 (MQ=255) aaaTCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTg < 1:247130/114‑1 (MQ=255) aTCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTAc < 2:1080763/116‑1 (MQ=255) aTCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTAc > 1:1080763/1‑116 (MQ=255) tCAGCAGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCg < 2:599504/139‑1 (MQ=255) agcagcagcAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGaa > 2:251573/1‑96 (MQ=255) agcagcagcAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGaa < 1:251573/96‑1 (MQ=255) agcagcagcAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGt > 1:1166534/1‑139 (MQ=255) | GGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGT > NZ_CP009273/2124186‑2124454 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |