| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,298,781 | T→C | 40.5% | K170R (AAA→AGA) | metF ← | 5,10‑methylenetetrahydrofolate reductase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,298,781 | 0 | T | C | 40.5% | 10.3 / 23.3 | 37 | K170R (AAA→AGA) | metF | 5,10‑methylenetetrahydrofolate reductase |
| Reads supporting (aligned to +/- strand): ref base T (15/7); new base C (5/10); total (20/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.98e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TCGACATCGAAGAAGAACTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGGA > minE/2298724‑2298848 | tCGACATCGAAGAAGAACTGAGTAATGGCGCGGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAg > 1:918890/1‑68 (MQ=255)tCGACATCGAAGAAGAACTGAGTAATGGCGCGGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAg > 1:2795630/1‑68 (MQ=255)tCGACATCGAAGAAGAACTGAGTAATGGCGCGGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAg > 1:1415470/1‑68 (MQ=255)tCGACATCGAAGAAGAACTGAGTAATGGCGCGGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAg > 1:1424857/1‑68 (MQ=255)tCGACATCGAAGAAGAACTGAGTAATGGCGCGGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAg > 1:1621915/1‑68 (MQ=255) cTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCt > 1:1081488/1‑69 (MQ=255) cTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCt > 1:868425/1‑69 (MQ=255) cTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCt > 1:2909590/1‑69 (MQ=255) cTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCt > 1:2516599/1‑69 (MQ=255) cTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCt > 1:2348330/1‑69 (MQ=255) cTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCt > 1:1900685/1‑69 (MQ=255) cTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCt < 1:1806242/69‑1 (MQ=255) cTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCt > 1:1047157/1‑69 (MQ=255) cTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCt > 1:1569972/1‑69 (MQ=255) cTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGACTGAGCGCt > 1:1213463/1‑69 (MQ=255) gAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCttt > 1:3226356/1‑69 (MQ=255) gAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCttt > 1:3148501/1‑69 (MQ=255) gAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCttt > 1:2107294/1‑69 (MQ=255) gAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCttt > 1:2767563/1‑69 (MQ=255) ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc < 1:1639675/68‑1 (MQ=255) ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc < 1:3318847/68‑1 (MQ=255) ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc < 1:3104833/68‑1 (MQ=255) ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc < 1:2783955/68‑1 (MQ=255) ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc < 1:177740/68‑1 (MQ=255) ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc < 1:1889088/68‑1 (MQ=255) ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc < 1:1631364/68‑1 (MQ=255) ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc < 1:2202705/68‑1 (MQ=255) ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc < 1:1989293/68‑1 (MQ=255) ggcgcgGTTGGCTCCGGCATCCACTTTTCGTCTCAGGTTGAGTAAATCGGCCTGGGCACTTTTCGCTTc < 1:2086624/68‑1 (MQ=255) ttGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGa < 1:1289613/69‑1 (MQ=255) tGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGa < 1:2375722/68‑1 (MQ=255) gCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGa < 1:439144/69‑1 (MQ=255) gCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGa < 1:2720106/69‑1 (MQ=255) gCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGa < 1:236746/69‑1 (MQ=255) aTCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGa < 1:525237/67‑1 (MQ=255) cACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATAcgc > 1:1201815/1‑69 (MQ=255) cACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATAcgc > 1:2820335/1‑69 (MQ=255) tttCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTcc > 1:2730277/1‑40 (MQ=255) tttCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCa > 1:2231394/1‑65 (MQ=255) tttCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCa > 1:391303/1‑65 (MQ=255) tttCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGga < 1:3355776/69‑1 (MQ=255) tttCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGga < 1:262072/69‑1 (MQ=255) tttCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGga < 1:2472799/69‑1 (MQ=255) tttCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGga < 1:2405833/69‑1 (MQ=255) tttCAGATTAAGCAAATCCGCCTGAGCGCTTTTT‑CTTCCGGGTGAACTTCCGGATAcgccg > 1:1154703/1‑61 (MQ=255) | TCGACATCGAAGAAGAACTGAGTAATCGCGCGGTTGGCTCCGGCATCCACTTTGCGTTTCAGATTAAGCAAATCCGCCTGAGCGCTTTTTGCTTCCGGGTGAACTTCCGGATACGCCGCCACGGA > minE/2298724‑2298848 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |