| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 4,424,386 | G→A | *229* (TAG→TAA) | ulaF → | L‑ribulose 5‑phosphate 4‑epimerase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,424,386 | 0 | G | A | 95.5% | 60.0 / ‑3.6 | 22 | *229* (TAG→TAA) | ulaF | L‑ribulose 5‑phosphate 4‑epimerase |
| Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (9/12); total (10/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.55e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.81e-01 | |||||||||||
AGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTAT > NC_000913/4424260‑4424507 | agaagTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGc > 1:324987/1‑139 (MQ=255) ggATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGga < 1:305460/139‑1 (MQ=255) tGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATa < 2:324987/139‑1 (MQ=255) cccGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc > 1:348817/1‑139 (MQ=255) cccGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc > 1:646086/1‑139 (MQ=255) caACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc < 1:25328/122‑1 (MQ=255) caACTCAATCACATCGACAGCTTCCTGATGAAAAAACACTTCATGCGTAAACACGGTCCGAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc > 2:25328/1‑122 (MQ=255) cTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGccc < 2:646086/139‑1 (MQ=255) cTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCCGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaag < 1:123655/139‑1 (MQ=255) ttCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaa < 2:490906/111‑1 (MQ=255) ttCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaa > 1:490906/1‑111 (MQ=255) tCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaagca < 1:287073/139‑1 (MQ=255) tCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaagca > 2:475742/1‑139 (MQ=255) gCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACt > 1:136915/1‑139 (MQ=255) cTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCg < 2:348817/139‑1 (MQ=255) cTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCg < 1:475742/139‑1 (MQ=255) aaCGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc > 2:344423/1‑62 (MQ=255) aaCGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc < 1:344423/62‑1 (MQ=255) gCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGt < 2:710079/139‑1 (MQ=255) tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCg < 2:352417/47‑1 (MQ=255) tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCg > 1:352417/1‑47 (MQ=255) tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTtat > 1:483983/1‑139 (MQ=255) | AGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTAT > NC_000913/4424260‑4424507 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |