Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP000730 | 408,719 | T→C | 16.6% | L491L (TTA→TTG) | metE ← | 5‑methyltetrahydropteroyltriglutamate‑‑ homocysteine S‑methyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP000730 | 408,719 | 0 | T | C | 16.6% | 63.8 / 5.5 | 30 | L491L (TTA→TTG) | metE | 5‑methyltetrahydropteroyltriglutamate‑‑ homocysteine S‑methyltransferase |
Reads supporting (aligned to +/- strand): ref base T (15/10); new base C (2/3); total (17/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.28e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01 |
CTGTTAAGCTTTGTGCATAAACTGTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCATTTTTTAAGAATGTTTCATA > CP000730/408592‑408856 | cTGTTAAGCTTTGTGCATAAACTGTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTa > 1:187079/1‑141 (MQ=255) tGTTAAGCTTTGTGCATAAACTGTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCCCCCAACCGAATTTAGTTACTAAGAAACCTTGTaa > 2:30855/1‑141 (MQ=255) tAAGCTTTGTGCATAAACTGGTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAAttt < 1:43470/141‑1 (MQ=255) tttGTGCATAAACTGTTTCATCAACAGGTAAAGGCGCTGTCCATTTTACATCACCATAATTGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGACACCTTGTAATTTTactc > 2:191787/1‑141 (MQ=255) gtgCATAAACTGTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCga > 1:44473/1‑141 (MQ=255) gCATAAACTGTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCgaag > 2:225929/1‑141 (MQ=255) gCATAAACTGTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACCAAGAAACCTTGTAATTTTTCTCCgaag > 2:247258/1‑141 (MQ=255) gCATAAACTGTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACCAAGAAACCTTGTAATTTTTCTCCgaag > 2:247245/1‑141 (MQ=255) aTAAACTGTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCgaagaa < 2:280813/141‑1 (MQ=255) aCTGTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCa < 1:286711/141‑1 (MQ=255) tGTTTCATCAACAGTTAAAGGCGCTGTCCCTTTTACATCGCCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAAc < 1:194012/141‑1 (MQ=255) gTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAAcc > 1:38811/1‑141 (MQ=255) aCCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATc > 2:256809/1‑141 (MQ=255) aaTGATTGGTTGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGcc < 1:50072/141‑1 (MQ=255) aaTGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAGCCATGTCATTACGTTCAAATTCCCCGTGAACTAATACATCTACGcc > 1:86434/1‑141 (MQ=255) aTGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAGATTCACCGTGAACTAATACATCTAAGCCa > 1:220927/1‑141 (MQ=255) ggtggtTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTTTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCTTGGACTAATACCTCTACGCCAACGtct > 2:206665/1‑141 (MQ=255) cGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTTAACTAATACATCTAAGCCAATGTCTTCTTGAAttt < 2:246130/141‑1 (MQ=255) ccGTGAACCATATGATTGCACCCAACCCAACTCAGTCACCAACAAGCCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAAt < 1:52747/140‑1 (MQ=255) gTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATcc > 1:101229/1‑141 (MQ=255) gCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAATTTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCAt > 2:102437/1‑141 (MQ=255) cccAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCAtttt > 2:159611/1‑141 (MQ=255) cAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCAtttttt < 1:225929/141‑1 (MQ=255) cAACCGAATTTAGTTACCAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCAtttttt < 1:247245/141‑1 (MQ=255) cAACCGAATTTAGTTACCAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCAtttttt < 1:247258/141‑1 (MQ=255) aaCCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCATTTTTTa < 1:188577/141‑1 (MQ=255) aCCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCATTTTTTaa < 2:38811/141‑1 (MQ=255) aaTTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTCCGTTCAAATTCACCGAGAACTACTACATCTAAGCCAATGTCATCTTGAAATTTAATCCATCGAGCAATATCAtttatgcagact > 1:273818/1‑132 (MQ=255) tttAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCATTTTTTAAGAATGt > 2:123121/1‑141 (MQ=255) tACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCATTTTTTAAGAATGTTTCata < 1:71982/141‑1 (MQ=255) | CTGTTAAGCTTTGTGCATAAACTGTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCATTTTTTAAGAATGTTTCATA > CP000730/408592‑408856 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |