| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 2,886,731 | T→C | I163M (ATA→ATG) | cas3 ← | CRISPR‑associated endonuclease/helicase Cas3 |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,886,731 | 0 | T | C | 90.9% | 88.9 / ‑1.3 | 34 | I163M (ATA→ATG) | cas3 | CRISPR‑associated endonuclease/helicase Cas3 |
| Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (15/15); total (17/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01 | |||||||||||
ATCCTCATTAAACAGAAAGGTATTCGTTGTAGTCCAGGAGCCTAACCAGTCAGCAAGCGAGCAAAAACCTGCTAACAGTGATGAACAATCAGGTGGTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCAC‑GAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAAC‑GCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGG > NC_000913/2886485‑2886971 | aTCCTCATTAAACAGAAAGTTATTCGTTGTAGTCCGAGATCCTAACCAGTCAGCAAGCGAGCAAAAACCTGCTAACAGTGATGAACAATCAGGTGGTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCTCTCCTCTC‑GAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGTAATTCCCCAAC‑GCGACTTATCTTGATCCTGGGAATGTAACATaa < 2:396304/251‑1 (MQ=255) aGCCTAACTAGTCAGCAAGCGAGCAAAAACCTGCTAACAGTGATGAACAATCAGGTGGTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCAC‑GAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAAC‑GCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAAc > 1:579047/1‑251 (MQ=255) aGCCTAACCAGTCAGCAAGCGAGCAAAAACCTGCTAACAGTGATGAACAATCAGGTGGTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCAC‑GAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAAC‑GCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAAc < 1:153488/251‑1 (MQ=255) aaaCCTGCTAACAGTGATGAACAATCAGGTGGTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCAC‑GAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAAC‑GCGACTTCTCTTGATCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGCAACCAGGAATCATAAGGATGAGGAGAGGc > 2:424079/1‑251 (MQ=255) ggatgaACAATCAGGTGTTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCAC‑GAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAAC‑GCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGa < 1:364844/250‑1 (MQ=255) gaACAATCAGGAAGTATATCGTTTATAGATAACCCCGCTGACGTTAAAACTAATGCTTCCAGTAGAGATATCCACTCCTCAC‑GAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCACAAC‑GCGACTTATCTTGCTCATGGGAATGTAACATAAAACCATGATGC‑CCTGTAACAGACACTACCCATGAGAACCAGGACTAATAAGGATGAGGAGCGGAATCAAAAAAACTGaaaa > 2:512438/1‑250 (MQ=255) ggtggtATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCAC‑GAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAAC‑GCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCg > 1:622552/1‑249 (MQ=255) cGTTTATAGAAAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCAC‑GAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAAC‑GCGACTTATCTTGAGCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGctct < 2:520975/251‑1 (MQ=255) aGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCAC‑GAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAAC‑GCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTg < 1:53200/251‑1 (MQ=255) cccGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCAC‑GAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAAC‑GCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTg > 1:460711/1‑251 (MQ=255) cTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTTAC‑GAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAAC‑GCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCATGGtt > 2:158993/1‑250 (MQ=255) ttAAAAATAATGCTTCCAGTACAGATATCCACTCCTCAC‑GAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAAC‑GCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCACGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAAcaa > 2:303502/1‑248 (MQ=255) aaaaaaTAATTCTGCCAGTACATATATCCACCCCT‑ACTGAGTCTGTTTATCGTGCGAAGAATAAGATGCCATAGAAGCTGGCATTTCCCAAC‑GCGACTTATCTTGAGCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCTTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTTAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAAt < 2:603919/250‑1 (MQ=255) aaaTAATGCTTCCAGTACAGATATCCACTCCTCAC‑GAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAAC‑GCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATAc > 1:463497/1‑250 (MQ=255) atGCTTCCAGTAAAGATAGCCCATCCGCAC‑GAGCCTGTTTATCTTGCGCAGCAAAAGATGCCAGAGAAGCTGGCATTTCCCAAC‑GCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCATGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCg < 2:82958/251‑1 (MQ=255) aCTCCTCAC‑GAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAAC‑GCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAAttt > 1:520115/1‑250 (MQ=255) ccTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAAC‑GCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCCTAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTtgt > 1:617236/1‑251 (MQ=255) gttatccttcgcagactaaggttccAGAGAAGGTGGCATTTCCCAACCGGTACTAACCTGGA‑CCCGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGtt < 1:569130/227‑1 (MQ=255) tGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAAC‑GCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATg < 1:303502/247‑1 (MQ=255) agcaTAAGATGCCAGAGAAGCTGGCATTTCCCAAC‑GCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTa > 1:359237/1‑250 (MQ=255) tAAGATGCCAATGAAGCGGGCCTTTCCCATC‑GCGACTTATCGTGATCCTGGGAATGTAACATAAAACCATGA‑GTGCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGCGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACATACCTGCTGCACCAGGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTg < 2:19882/250‑1 (MQ=255) ccAGAGAAGCTGGCATTTCCCAAC‑GCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAACAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGtt > 1:85647/1‑251 (MQ=255) agagAAGCTGGCATTTCCCAAC‑GCGACTTATCTTGATCCTGTGAATGTAATATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCg > 1:309176/1‑189 (MQ=255) agagAAGCTGGCATTTCCCAAC‑GCGACTTATCTTGATCCTGTGAATGTAATATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCg < 2:309176/189‑1 (MQ=255) agaatcgtgCATGTCCCAAC‑GCGACTTACCTTGATCCTGGGAATTCAACATAAAACCATGATGT‑CCTGTAGCGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATGGCTCTGAAAGTGAATCCTGGTTAATCCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTTTTGATTGACCATTAAGTGATGGCGTTGCa < 2:84371/242‑1 (MQ=255) gCTGGCATTGCCCAAC‑GCGACTTATCTTGATCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGAtt < 1:661417/251‑1 (MQ=255) aTTTCCCAAC‑GCTACTTACCTTGATCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTTTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAAttt < 1:64360/251‑1 (MQ=255) ttATCTTGATCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTctgc > 1:163623/1‑251 (MQ=255) aTCTTGATCCTGGGAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTACGTGATGGCGTTGCAGGATTTACTTTCAGCCAACTTTctgctg > 1:157101/1‑251 (MQ=255) gATCCTGGGAATGTAACATAAAACCATGATGT‑CCTGGAACGGCCTCTAACCATGGAAACCAGGACGCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCACCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAAACAAGAAAGACCGGATGCACCATGAGTAAAGTGACGGAACATGTGTGTTGATGGGCCATTAAGTGAGGGCGTTGAAGTATTTAATTCCAGTCAACTTTTTGGTGAtta < 2:96566/250‑2 (MQ=255) aTCCTGGGAATGTAATATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTTATGGACCATTAAGGGATGGCGTTGCAGGGTTTAATTTCAGCCAACttccttctactttat > 2:354382/1‑238 (MQ=255) ccTGGGAATGTAACATAAAACCATTATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTatat < 2:163623/251‑1 (MQ=255) gggAAGGTAAAAGAAACCCATGATGG‑CGTGTAACGGCCTCATCCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGGTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTgg < 2:460711/251‑1 (MQ=255) ggAATGTAACATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTgg > 1:666400/1‑250 (MQ=255) | ATCCTCATTAAACAGAAAGGTATTCGTTGTAGTCCAGGAGCCTAACCAGTCAGCAAGCGAGCAAAAACCTGCTAACAGTGATGAACAATCAGGTGGTATATCGTTTATAGATAACCCCGCTGGCGTTAAAAATAATGCTTCCAGTACAGATATCCACTCCTCAC‑GAGCCTGTTTATCTTGCGCAGCATAAGATGCCAGAGAAGCTGGCATTTCCCAAC‑GCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGT‑CCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGG > NC_000913/2886485‑2886971 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |