| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 880,747 | A→C | 67.7% | intergenic (+66/‑128) | yciM → / → pyrF | conserved hypothetical protein/orotidine‑5'‑phosphate decarboxylase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 880,747 | 0 | A | C | 67.7% | 34.7 / 18.3 | 34 | intergenic (+66/‑128) | yciM/pyrF | conserved hypothetical protein/orotidine‑5'‑phosphate decarboxylase |
| Reads supporting (aligned to +/- strand): ref base A (8/3); new base C (12/11); total (20/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.95e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01 | |||||||||||
ATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCAC > minE/880681‑880807 | aTTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAac < 1:2625722/68‑1 (MQ=255)aTTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAac < 1:1758095/68‑1 (MQ=255) tttAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAacgacg > 1:2505746/1‑69 (MQ=255) tttAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAacgacg > 1:2361791/1‑69 (MQ=255) tttAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAacgacg > 1:2680813/1‑69 (MQ=255) aaaTCCGACTTTAGTTACAACATACTAATTATTGATGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa < 1:1763014/69‑1 (MQ=38) aaaTCCGACTTTAGTTACAACATACTAATTATTGATGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa < 1:2587617/69‑1 (MQ=38) aaaTCCGACTTTAGTTACAACATACTAATTATTGATGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa < 1:2026931/69‑1 (MQ=38) aaaTCCGACTTTAGTTACAACATACTAATTATTGATGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa < 1:2993575/69‑1 (MQ=38) aaaTCCGACTTTAGTTACAACATACTAATTATTGATGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa < 1:3076454/69‑1 (MQ=38) aaaTCCGACTTTAGTTACAACATACTAATTATTGATGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa < 1:1436205/69‑1 (MQ=38) aaaTCCGACTTTAGTTACAACATACTAATTATTGATGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa < 1:3135649/69‑1 (MQ=38) aaaTCCGACTTTAGTTACAACATACTAATTATTGATGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa < 1:3211392/69‑1 (MQ=38) aaaTCCGACTTTAGTTACAACATACTAATTATTGATGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa < 1:3222331/69‑1 (MQ=38) aaaTCCGACTTTAGTTACAACATACTAATTATTGATGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa < 1:3417501/69‑1 (MQ=38) acaacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAacgacg > 1:126230/1‑45 (MQ=255) acaacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAg > 1:1141794/1‑57 (MQ=255) acaacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAg > 1:2740941/1‑57 (MQ=255) acaacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAg > 1:2554712/1‑57 (MQ=255) acaacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAg > 1:695882/1‑57 (MQ=255) acaacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAg > 1:1039924/1‑57 (MQ=255) acaacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAg > 1:2054759/1‑57 (MQ=255) acaacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAg > 1:1279016/1‑57 (MQ=255) acaacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAg > 1:1685367/1‑57 (MQ=255) aTGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa < 1:633800/35‑1 (MQ=25) cTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCAc < 1:2901481/69‑1 (MQ=255) cTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCAc < 1:1801186/69‑1 (MQ=255) cTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCAc < 1:1432217/69‑1 (MQ=255) ggCACCGACGGCGC‑GGAAAAGCCTGTCAGAGGAGAAATCGCAACTGTTAAtt > 1:1201482/1‑52 (MQ=255) ggCACCGACGGCGC‑GGAAAAGCCTGCCAGAGGAGAAATCGCACCTGTTAAtt > 1:2380115/1‑52 (MQ=25) ggCACCGACGGCGC‑GGAAAAGCCTGCCAGAGGAGAAATCGCAACTGtt > 1:3088495/1‑48 (MQ=25) ggCACCGACGGCGC‑GGAAAAGCCTGCCAGAGGAGAAATCGCAACTGTTAAtt > 1:294435/1‑52 (MQ=37) ggCACCGACGGCGC‑GGAAAAGCCTGCCAGAGGAGAAATCGCAACTGTTAAtt > 1:1251027/1‑52 (MQ=37) ggCACCGACGGCGC‑GGAAAAGCCTGCCAGAGGAGAAATCGCAACTGTTAAtt > 1:550660/1‑52 (MQ=37) ggCACCGACGGCGC‑GGAAAAGCCTGCCAGAGGAGAAATCGCAACTGTTAAtt > 1:627587/1‑52 (MQ=37) ggCACCGACGGCGC‑GGAAAAGCCTGCCAGAGGAGAAATCGCAACTGTTAAtt > 1:2167438/1‑52 (MQ=37) ggCACCGACGGCGC‑GGAAAAGCCTGCCAGAGGAGAAATCGCAACTGTTAAtt > 1:674426/1‑52 (MQ=37) ggCACCGACGGCGC‑GGAAAAGCCTGCCAGAGGAGAAATCGCAACTGTTAAtt > 1:2418404/1‑52 (MQ=37) ggCACCGACGGCGC‑GGAAAAGCCTGCCAGAGGAGAAATCGCAACTGTTAAtt > 1:967431/1‑52 (MQ=37) ggCACCGACGGCGC‑GGAAAAGCCTGCCAGAGGAGAAATCGCAAATGTTAAtt > 1:1870687/1‑52 (MQ=25) | ATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCAC > minE/880681‑880807 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |