| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | CP009273 | 2,111,876 | G→C | 53.9% | L4V (CTT→GTT) | wcaK ← | colanic acid biosynthesis protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP009273 | 2,111,876 | 0 | G | C | 53.9% | 4.8 / 54.3 | 39 | L4V (CTT→GTT) | wcaK | colanic acid biosynthesis protein |
| Reads supporting (aligned to +/- strand): ref base G (13/5); new base C (11/10); total (24/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.23e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01 | |||||||||||
GCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCGTGTGGATTGAGAATGTTGATGGCATCAAGTAAGCCGCGCAGGATGGCGCTGTCGCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCG > CP009273/2111738‑2112018 | gcagcCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCGTGTGGATTGAGAATGTTGATGGCATCAAGTAAGCCGCGCAGGATGGCGCTGTCGCCACGATTGCCGCAAGGGTGGTTGCCCAGGATAACTagt > 1:488930/1‑141 (MQ=255) aagaaCTGTCCGGATCGCGGCCCATCACATGCACTTCGGCGTGTGGATTGAGAATGTTGATGGCATCAAGTAAGCCGCGCAGGATGGCGCTGTCGCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAACTAATTTCATTTTTTCCTCa < 2:238305/149‑1 (MQ=255) aCTGACCGGATAGCGGCTCATCACATCCACTTCGGCGTGTGGATTGAGAATGTTGATGGCATCAAGTAAGCCGCGCAGGATGGCGCTGTCGCCACGGTTGCCGCAAGTGGGGTTGCCCAGAGTAGGTGATTTCATTTTTTCCCCAGaaa > 2:47143/1‑149 (MQ=255) tAGCGGCTCATCACATCCACTTCGGCGTGTGGATTGAGAATGTTGATGGCATCAAGTAAGCCGCGCAGGATGGCGCTGTCGCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAACTAATTTCATTTTTTCCTCATAAATTTGATGCCa < 2:50180/149‑1 (MQ=255) atcaCATCCACTTCGGCGTGTGGATTGAGAATGTTGATGGCATCAAGTAAGCCGCGCAGGATGGCGCTGTCGCCACGATTGCCGCAAGTGTGGTTGCCCAGACTAAGTGATTTCATTTTTTCCTCATAAATTTGATGCCAGGt > 1:59860/1‑143 (MQ=255) atcaCATCCACTTCGGCGTGTGGATTGAGAATGTTGATGGCAGCAAGTAAGCCGCGCAGGATGGCGCTGTCGCCACGATTGCCGGAAGAGTGGTTGCCCAGAATAGGCAATGTCCTTTTTTCCTCATCAATTTGATGCCAGGt > 1:1543/1‑143 (MQ=255) cacgtgggCGTGTGGATTGAGAATGTTGATGGCATCAAGTAAGCCGCGCAGGATGGCGCTGTCGCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAACTAATTTCATTTTTTACTCATTAATTTGATGGCAGGTGAGGCTGc < 1:558031/137‑1 (MQ=255) aCTTCGGCGTGGTGATTGCGAATGTGGATGGCATCAAGTAAGCCGCGCAGGATGGCGCTGTCGCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAACTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCgtttgtt < 2:23677/149‑1 (MQ=255) gcGGCGTGTGGATTGAGAATGTTGATGGCATCAAGTAAGCCGCGCAGGATGGCGCTGTCGCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCgttt < 1:506852/142‑1 (MQ=255) gCGTGTGGATTGAGAATGTTGATGGCATCAAGTAAGCCGCGCAGGATGTCACTGTCGCCGCGATTGCCGCAAGTGTGGTTGCCCAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGGGAGGCCGCgtttgtt > 1:232962/1‑143 (MQ=255) tgtgGATTGAGAATGTTGATGGCATCAAGTAAGCCGCGCAGGATGGCGCTGTCGCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAACTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGtt > 2:51916/1‑149 (MQ=255) gtgGATTGAGAATGTTGATGGCATCAAGTAAGCCGCGCAGGATGGCGCTGTCGCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAAGTAATTTCATTTTTTCCTCATAAAGTTGATGCCAGGTGAGGCTGCGtttgttttgt > 1:428650/1‑143 (MQ=255) tGAGAATGTTGATGGCATCAAGTAAGCCGCGCAGGATGGCGCTGTCGCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAACTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTg > 1:272654/1‑143 (MQ=255) gagaATGTTGATGGCATCAAGTAAGCCGCGCAGGATGGCGCTGTCGCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCttt > 2:15252/1‑149 (MQ=255) gTTGATGGCATCAAGTAAGCCGCGCAGGATGGCGCTGTCGCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCttt < 1:510750/143‑1 (MQ=255) gATGGCATCAAGTAAGCCGCGCAGGATGGCGCTGTCGCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGAt > 2:597692/1‑149 (MQ=255) ctGGCCTCAAGTAAGCCGCGCAGGATGGCGGTTGCGCCACGATTGCCGCAAGTGCGGTTGTCCAGAATAACTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGc < 1:31770/142‑1 (MQ=255) aTGGCATCAAGTAAGCCGCGCAGGATGGCGCTGTCGCCACGATTTGCTCAAGTGTGGTTGCCCAGAATAACTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATg < 2:276887/149‑1 (MQ=255) ggCATCAAGTAAGCCGCGCAGGATGGCGCTGTCGCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCg > 2:245058/1‑149 (MQ=255) tCAAGTAAGCCGCGCAGGATGGCGCTGTCGCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGc < 2:170152/149‑1 (MQ=255) aGTAAGCCGCGCAGGATGGCGCTGTCGCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAACTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGCTGCAACGCTGa > 2:541044/1‑149 (MQ=255) aGCCGCGCAGGATGGCGCTGTCGCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAACTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACCGTGAcccg > 2:545064/1‑149 (MQ=255) ccGCGCAGGATGGCGCTGTCGCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTgc > 1:164292/1‑142 (MQ=255) cTGTCGCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAACTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGCCGCTGACGCGTCTTATCAGGc > 1:50666/1‑143 (MQ=255) tGTCGCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGCTGCGACGCTGACGCGTCTTATCGGGCCTACAAg > 2:65500/1‑149 (MQ=255) tCGCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAACTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCCGGCCTa > 1:170333/1‑143 (MQ=255) gCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCGGACGCGTCTTAGCAGGCCTAGAAGAccc > 2:141914/1‑149 (MQ=255) cGATTGCCGCAAGTGTGGTTGCCCAGAATAACTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGc < 2:111037/149‑1 (MQ=255) gATTGCCGCAAGTGTGGTTGCCCAGAATAACTAATTTCATTTTTGCCGCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGca < 2:234561/149‑1 (MQ=255) aTTGCCGCAAGTGTGGTTGCCCAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGcac > 2:234681/1‑149 (MQ=255) tGCCGCAAGTGTGGTTGCCCAGAATAACTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCGTGCCGACAAGACCCGAGCACAg > 2:14566/1‑149 (MQ=255) gCCGCAAGTGTGGTTGCCCAGAATAACTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAg < 1:365542/143‑1 (MQ=255) cGCAAGTGTGGTTGCCCAGAATAACTAATTTCATTTTTTTTTTTTTTTTTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGca > 1:505728/1‑143 (MQ=37) gtgtGGTTGCCCAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGaa < 1:540092/143‑1 (MQ=255) tgtgGTTGCCCAGAATAACTAATTTCATTTTTTCCTCATAAATTTGATGCCGGATGAGGGTGGGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAgg > 2:282159/1‑149 (MQ=37) tgGTTGCCCAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGAc < 2:335617/149‑1 (MQ=255) gCCCAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGAAGCCGCTCTAGCCGAGCAGGCCTAAACGGCCCAATCACAGAACCGTAGGACGGATa > 2:461504/1‑149 (MQ=255) cAGAATAACTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGACCACAGAACAGTCGCACgg > 1:472415/1‑143 (MQ=38) gAATAACTAATTTCAGTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGATGCCGCTTTTTGCCCCAGTCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCg < 2:307834/149‑1 (MQ=37) | GCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCCACTTCGGCGTGTGGATTGAGAATGTTGATGGCATCAAGTAAGCCGCGCAGGATGGCGCTGTCGCCACGATTGCCGCAAGTGTGGTTGCCCAGAATAAGTAATTTCATTTTTTCCTCATAAATTTGATGCCAGGTGAGGCTGCGTTTGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCG > CP009273/2111738‑2112018 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |