| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | CP000730 | 506,009 | C→A | 14.1% | P137Q (CCA→CAA) | amyA → | alpha, alpha‑phosphotrehalase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 506,009 | 0 | C | A | 14.1% | 62.5 / 2.2 | 31 | P137Q (CCA→CAA) | amyA | alpha, alpha‑phosphotrehalase |
| Reads supporting (aligned to +/- strand): ref base C (16/10); new base A (2/2); total (19/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01 | |||||||||||
AAGACTTGAAAGTGATGTTAGATATTGTCATTAATCATACGTCGACGGAGCATGAATGGTTTAAAGAAGCCCGTAAATCTAAAGATAACCCTTATAGAGATTATTACTTTTTCAGATCATCTGAAGACGGGCCGCCAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGAAGTATGATTCTGAGACAGATGAATATTATTTACATTTATTTGATGTCAGTCAAGCTGATTTAAATTGGGATAATCCGGAAGTACGTCAATCGTTATAT > CP000730/505874‑506145 | aaGGCTTGAATGGGATGTTAGGTATTGTCATTAATCATACGTCGACGGTGCATGAATGGTTTAAAGAAGCCCGTACATCTAAAGATTACCCTTCTAGTAATTCTTTCTTTTTCAGATCTTCTGAAGGCGGGCCGCcaacaa < 2:259488/141‑1 (MQ=255) gAAAGTGATGTTAGATATTGTCATTAATCATACGTCGACGGAGCATGAATGGTTTAAAGAAGCCCGTAAATCTAAAGATAACCCTTATAGAGATTATTACTTTTTCAGATCATCTGAAGACGGGCCGCCAACAAATTGGTa > 1:201946/1‑141 (MQ=255) aatgTGATGTTAGATATTGTCATTAATCATACGTCGACGGAGCATGAATGGTTTAAAGAAGCCCGTAAATCTAAAGATAACCCTTATAGAGATTATTACTTTTTAAGATCATCTGAAGACGGGCCGGGAACAAATTGGTAt > 2:73583/4‑141 (MQ=255) aGTGATGAGAGATATTGTCATTAATCATACGTCGACGGAGCATGAATGGTTTAAAGAAGCCCGTAAATCTAAAGATAACCCTTATAGAGATTATTACTTTTTCAGATCATCTGAAGACGGGCCGCCAACAAATTGGTATTc < 2:218343/141‑1 (MQ=255) tGATGTTAGATATTGTCATTAATCATACGTCGACGGAGCATGAATGGTTTAAAGAAGCCCGTAAATCTAAAGATAACCCTTATAGAGATTATTACTTTTTCAGATCATCTGAAGACGGGCCGCCAACAAATTGGTATTCTa < 1:51490/141‑1 (MQ=255) gTTAGATATTGTCATTAATCATACGTCGACGGAGCATGAATGGTTTAAAGAAGCCCGTAAATCTAAAGATAACCCTTATAGAGATTATTACTTTTTCAGATCATCTGAAGACGGGCCGCCAACAAATTGGTATTCTAAAtt > 1:225860/1‑141 (MQ=255) gTTAGATATTGTCATTAATCATACGTCGACGGAGCATGAATGGTTTAAAGAAGCCCGTAAATCTAAAGATAACCCTTATAGAGATTATTACTTTTTCAGATAATCTGAAGACGGGCCGCCACCAAATCGGTATTCTAACtt > 2:98886/1‑141 (MQ=255) tCATTAATCATACGTCAACGGAGCATGAATGGTTTACAGAAGCCCGTAAATCTAAAGATAACCCTTATAGAGATTATTACTTTTTCAGATCATATGGAGACGCGCGGCCAACAAATTGGTATTCTAAATTCGGTGCTAATg > 2:98916/1‑141 (MQ=255) aTACGTCGACGGAGCATGAATGGTTTAAAGAAGCCCGTAAATCTAAAGATAACCCTTATAGAGATTATTACTTTTTCGGCACTCCTGGAGACGGGCCGCCAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGAAGt > 1:278942/1‑141 (MQ=255) tACGTCGACGGAGCATGAATGGTTTAAAGAAGCCCGTAAATCTAAAGATAACCCTTATAGAGATTATTACTTTTTCAGATCATCTGAAGACGGGCCGCCAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGAAGTa < 2:68128/141‑1 (MQ=255) tCGACGGAGCATGAATGGTTTAAAGAAGCCTGTAAATTTAAAGATACCCCTCATAGAGATTATTACTTTTTCAGATCATCTGAAGAGGGGCCGCCAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGGAGTATGAt < 2:258447/141‑1 (MQ=255) gAGCATGAATGGTTTAAAGAAGCCCGTAAATCTAAAGATAACCCTTATAGAGATTATTACTTTTTCAGATCATCTGAAGACGGGCCGCCAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGAAGGATGATTCTgag > 2:251132/1‑141 (MQ=255) aTGAATGGTTTAAAGAAGCCCGTAAATCTAAAGATAACCCTTATAGAGATTATTACTTTTTCAGATCATCTGAAGACGGGCCGCAAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGAAGTATGATTCTGAGACAg > 1:223299/1‑141 (MQ=255) aTGAATGGTTTAAAGAAGCCCGTAAATCTAAAGATAACCCTTATAGAGATTATTACTTTTTCAGATCACCTGAAGACGGGCCGCAAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGAAGTATGATTCTGAGACAg > 1:223268/1‑141 (MQ=255) cccGTAAATCTAAAGATAACACTTATAGAGATTATTCCTTTTTCAGCTCATCTGAAGACGGCACGCCAACACATTGGTATTCTAAATTAGGTGGTAATGCCTGGAAGTATGATTCTGAGACCGATGACTATAATTTACact > 2:148689/1‑139 (MQ=255) gTAAATCAAAAGCTAACCCTTATAGCTGTTTTTACTTTTTCAGATCATCTGAAGACGGCCCGCCAACCAATTGTTATTCTAAATTCTGTTGTAATGCATGCAAGTATGATTCTGAGACAGATGAATATTATTTACatttat < 2:155839/141‑1 (MQ=255) tAAATCTAAAGATAACCCTTATAGAGATTATTACTTTTTCAGATCATCTGAAGACGGGCCGCCAACAAATTGGTATTCTAAATTCGGTGGTAGTGCATGGAAGTATGTTTCTGAGACGGATGAATATTATTTACttttatt > 1:91493/1‑141 (MQ=255) tAAAGATAACCCTTATAGAGATTATTACTTCTTCAGATCATCTGAAGACGGGCCGCCAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGAAGTATGATTCTGAGACAGATGAATATTATTTACATTTATTTGATgt < 1:94409/141‑1 (MQ=255) aaCCCTTATAGAGATTATTACTTTTTCAGATCATCTGAAGACGGGCCGCCAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGAAGTATGATTCTGAGACAGATGAATATTATTTACATTTATTTTATGTCAGTCaa > 1:237109/1‑141 (MQ=255) tataGAGATTATTACTTTTTCAGATCATCTGAAGCCGGGCCGCCAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGAAGTATGATTCTGAGACAGATGAATATTATTTACATTTATTTGATGTCAGTCAAGCTGAt < 1:6785/141‑1 (MQ=255) gagagATTACTACTTTTTTAGAACATCTGAAGACGGGCCGCCAACAAAGTGGTATTTTAAATTCGGTGGTAATGCATGGAAGTATGATTCTGAGACAGATGAATATTATTTACATTTATTTGATGTCAGTCAAGCTGAttt < 2:7729/140‑1 (MQ=255) gagaTTATTACTTTTTCAGATCATCTGAAGACGGGCCGCCAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGAAGTATGATTCTGAGACAGATGAATATTATTTACATATATTTGATGTCCGCAAAGCTGAtttca > 2:63334/1‑139 (MQ=255) gaTTATTACTTTTTCAGATCATCTGAAGACGGGCCGCCAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGAAGTATGATTCTGAGACAGATGAATATTATTTACATTTATTTGATGTCAGTCACGCTGAGTTGAAt > 1:82868/1‑141 (MQ=255) gaTTATTACTTTTTCAGATCATCTGAAGACGGGCCGCCAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGAAGTATGATTCTGAGACAGATGAATATTATTTACATTTATTTGATGTCAGTCAAGCTGATTTAAAt > 1:162441/1‑141 (MQ=255) tattaCTTTTTCATATCATCTGAAGACGGGCCTCCAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGAAGTATGATTCTGAGACAGATGAATATTATTTACATTTATTTGATGTCAGTCAAGCTGATTTAAATTgg < 2:90964/141‑1 (MQ=255) attaCTTTTTCAGATCATCTGAAGACGGGCCGCCAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGAAGTATGATTCTGAGACAGATGAATATTATTTACATTTATTTGATGTCAGTCAAGCTGATTTAAATTggg > 2:137910/1‑141 (MQ=255) ttCAGATCATCTGAAGACGGGCCGCAAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGAAGTATGATTCTGAGACAGATGAATATTATTTACATTTATTTGATGTCAGTCAAGCTGATTTAAATTGGGATAATCCg < 2:223299/141‑1 (MQ=255) ttCAGAGCAGCTGAAGACGGGCCGCAAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGAAGTATGATTCTGAGACAGATGAATATTATTTACATTTATTTGATGTCAGTCAAGCTGATTTAAATTGGGATAATCCg < 2:223268/141‑1 (MQ=255) tcatcTGAAGACGGGCCGCCAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGAAGTATGATTCTGAGACAGATGAATATTATTTACATTTATTTGGTGTCAGTCAAGCTGATTTAAATTGGGATACTCCGGAAGTa > 2:269895/1‑141 (MQ=255) gAAGACGGGCCGCCAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGAAGTATGATTCTGAGACAGATGAATCTTATTTACATTTATTTGATGTCAGTCAAGCTGATTTAAATTGGGATAATCCGGAAGTACGTCaa > 2:205340/1‑141 (MQ=255) gAAGACGGGCCGCCAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGAAGTATGATTCTGAGACAGATGAATATTATTTACATTTATTTGATGTCAGTCAAGCTGATTTAAATTGGGATAATCCAGCAGTACGTaaa > 2:228421/1‑141 (MQ=255) gccGGGCCGCCAACAAATTGGTACTCTAAACTCGGGGGTAATGCATGGAAGTATGATTCGGAGACAGATGAATATTATTTACATTTATTTGTTGTCAGTCAAGCTGATTTAAATTGGGATAATCCGGAAGTACGTCAATCg < 1:98886/139‑1 (MQ=255) cctccAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGAAGTATGATTCTGAGACAGATGAATATTATTTACATTTATTTGATGTCAGTCAAGCTGATTTAAATTGGGATAATCCGGAAGTACGTCAATCGTTatat < 2:168445/138‑1 (MQ=255) | AAGACTTGAAAGTGATGTTAGATATTGTCATTAATCATACGTCGACGGAGCATGAATGGTTTAAAGAAGCCCGTAAATCTAAAGATAACCCTTATAGAGATTATTACTTTTTCAGATCATCTGAAGACGGGCCGCCAACAAATTGGTATTCTAAATTCGGTGGTAATGCATGGAAGTATGATTCTGAGACAGATGAATATTATTTACATTTATTTGATGTCAGTCAAGCTGATTTAAATTGGGATAATCCGGAAGTACGTCAATCGTTATAT > CP000730/505874‑506145 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |