Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 2,516,091 | C→T | 14.3% | M580I (ATG→ATA) | yfeA ← | putative diguanylate cyclase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,516,091 | 0 | C | T | 14.3% | 62.2 / 3.5 | 28 | M580I (ATG→ATA) | yfeA | putative diguanylate cyclase |
Reads supporting (aligned to +/- strand): ref base C (12/12); new base T (2/2); total (14/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01 |
TGCTCGATGTTGTACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAAGGGGCAGAAACTTATCCGGGGTCATAAT > NC_000913/2515945‑2516237 | tgctCGATGTTGTACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGCCCGCCTCCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGACGCAGCGCGAGCGGCag > 1:7596/1‑148 (MQ=255) gttgtACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGa > 1:441202/1‑149 (MQ=255) ttggACATGCTGCTTTCTGACTTGGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGAc < 1:14917/149‑1 (MQ=255) gtACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTACTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTg > 2:470075/1‑149 (MQ=255) tgctTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCCTTAAATTGACTGAAAAgcgc > 1:714148/1‑149 (MQ=255) tGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGgtg > 1:109738/1‑149 (MQ=255) ggTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTgtca > 2:221289/1‑149 (MQ=255) tGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTgccagc < 1:267394/149‑1 (MQ=255) gATCTCAAGAATGACCGCCGGCGGAGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTTAGCGGCATTAAATTGACTGACAAGCACGTTGCTTTTTGGTCGCAAGGGTGTGTTgccagcc < 2:163857/149‑1 (MQ=255) tGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAg < 2:441202/149‑1 (MQ=255) ggAGATGTGATACCGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCCTTAAATTGACTGAAACGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACtt < 1:593584/149‑1 (MQ=255) gaTGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACCGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCa > 1:319974/1‑149 (MQ=255) ttAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAAcgcgc > 1:309869/1‑149 (MQ=255) tAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAAcgcgcg < 1:248398/149‑1 (MQ=255) cAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTaa < 2:542190/149‑1 (MQ=255) aGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAg > 1:164096/1‑149 (MQ=255) ataatCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTaaa < 1:300657/149‑1 (MQ=255) ccGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAg < 1:470075/149‑1 (MQ=255) cGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGTATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGc > 1:330395/1‑149 (MQ=255) cGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGTATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGc > 1:330489/1‑149 (MQ=255) cGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGc > 1:330478/1‑149 (MQ=255) ttttGCAGCAGCGTGAGCGGTATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAAGGGGCAGaaa < 2:330395/149‑1 (MQ=255) ttttGCAGCAGCGTGAGCGGTATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAAGGGGCAGaaa < 2:330489/149‑1 (MQ=255) ttttGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAAGGGGCAGaaa < 2:330478/149‑1 (MQ=255) gcagcaGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAAGGGGCAGAAACTTa < 2:309869/149‑1 (MQ=255) cagcGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGATGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAAGGGGCAGAAACTTATCCg > 1:499686/1‑149 (MQ=255) gCGGCATTAAATTGACTGAAAAGCCCGGTCCTTTTTTGTCGCAAGCGTTTGTCGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATATGGGGCAGAAACTTATCCGGGGTCATa < 2:321472/149‑1 (MQ=255) ggCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCTAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTGAGGTTAAACTGAGCCATAAGGGGCAGAAAATTATCCGGGGTCCTAAt > 2:736208/1‑149 (MQ=255) | TGCTCGATGTTGTACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTGATAACGCTTAAACAGACGAATAATCCGCCCGGCAATATTCTTTTGCAGCAGCGTGAGCGGCATTAAATTGACTGAAAAGCGCGGTCCTTTTTTGTCGCAAGGGTGTGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAAGGGGCAGAAACTTATCCGGGGTCATAAT > NC_000913/2515945‑2516237 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |