| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | CP000730 | 1,318,237 | T→A | 43.6% | F124Y (TTT→TAT) | miaB → | 2‑methylthioadenine synthetase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 1,318,237 | 0 | T | A | 43.6% | 2.4 / 23.6 | 23 | F124Y (TTT→TAT) | miaB | 2‑methylthioadenine synthetase |
| Reads supporting (aligned to +/- strand): ref base T (8/5); new base A (5/5); total (13/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.85e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01 | |||||||||||
AATGCACATGACACTGAGGTCATTGCTGGTATACTTGAAGCATTAGGCTAT‑CAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTTTAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAATTTTAAAATCGTATCAAAATGTAGATATGATATTT > CP000730/1318104‑1318373 | aatgCACATGACACTGAGGTCATTGCTGGTATACTTGAAGCATTAGGCTAT‑CAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGaa > 2:81054/1‑141 (MQ=255) cacaTGACACTGAGGTCATTACTGGTATACTTGAAGCATTAGGCTAT‑CAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAAGAAAGGGATTAGtagaagag > 2:159640/1‑134 (MQ=255) acacGACACTGAGGTAAAAGATAGTATGCTTTTAGCATTAGTTTAT‑CTAGCAACGACTGATATTAACACCGCAGAGGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAgg < 1:118923/141‑1 (MQ=255) caTGACACTGAGGTCATTGCTGGTATACTTGAAGCATTAGGCTAT‑CAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTTTAGTGAAATAGGt > 2:140762/1‑141 (MQ=255) gACACTGAGGTCATTGCTGGTATACTTGAAGCATTAGGCTAT‑CAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAAAAAAGTGTTTAGTGAAATAGGTAAt < 1:179732/141‑1 (MQ=255) gggaaTTGCGGGTAGAGTGGAAGCATTAGGCTAT‑CAAGGAAAGACTGATATTAACACTGGAGGTGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAAGAAAGTGTATAGTGAAATAGGTAATTTGAAGCa < 2:256463/137‑1 (MQ=255) tGGTATACTTGAAGCATTAGGCTAT‑CAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGaaaaa < 2:51488/141‑1 (MQ=255) taaattgggaGCATTAGG‑TATACAAGCAAAGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGCTTAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGa < 2:139930/132‑1 (MQ=255) ataCTTGAAGCATTAGGCTAT‑CAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTTTAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGaa < 2:90566/141‑1 (MQ=255) tAGGCTAT‑CAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATAttt > 1:116211/1‑141 (MQ=255) tAGGCTAT‑CAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATAttt > 1:116207/1‑141 (MQ=255) atTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATg < 1:81054/141‑1 (MQ=255) taccacTTCAGATGTTATTTTAATGAACACATGTGCGATAAGACAAAATTCCGAGAACAAAGTGTTTAGTGAAATATGTAATTTAAAAAAAAAAAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTc < 2:268501/138‑1 (MQ=255) aatggaGATGTTATTTTATTGAATATATGTGCTATTAGAGAAAATGCTGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAATTATTTGAAAAAAGAACGACCTTATATTTTAAGCGGTGTTTGTGGTTGTATGTCACaa < 2:55917/136‑1 (MQ=255) cTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTTTAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACaag > 1:227209/1‑141 (MQ=255) aGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTTTAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAagaaga < 1:140762/141‑1 (MQ=255) ttATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTTTAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCag > 2:275624/1‑141 (MQ=255) ttttAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTTTAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCagtag > 1:139127/1‑141 (MQ=255) aCATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTTTAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAAtt > 2:232231/1‑141 (MQ=255) cGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAATTTTAAAATCGTATCAAAATGTAga > 2:16524/1‑141 (MQ=255) agaACAAAGTGTTTAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAATAGTGAATAAAATTTTAAAATCGTATAAAAATGTAgata > 2:256745/1‑141 (MQ=255) aaGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATGGGTGTTTGTGGTTGTATGTCACAAGAAGAGTGAGTCGTGAATAAAATTTTAAAATCGTATCAAAATGTAgctatgatat > 1:108101/1‑141 (MQ=255) gtgtTTAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAATTTTAAAATCGTATCAAAATGTAGATATGATAttt > 1:56955/1‑141 (MQ=255) | AATGCACATGACACTGAGGTCATTGCTGGTATACTTGAAGCATTAGGCTAT‑CAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTTTAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAATTTTAAAATCGTATCAAAATGTAGATATGATATTT > CP000730/1318104‑1318373 |
| 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: ‑ |