| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,839,658 | C→A | 100% | G94G (GGG→GGT) | surE ← | broad specificity 5'(3')‑nucleotidase and polyphosphatase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,839,658 | 0 | C | A | 81.3% | 60.0 / 25.0 | 43 | G94G (GGG→GGT) | surE | broad specificity 5'(3')‑nucleotidase and polyphosphatase |
| Reads supporting (aligned to +/- strand): ref base C (4/4); new base A (35/0); total (39/4) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.67e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.16e-01 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
ACGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGATTCGGCCCGGCGTTAATTCCGGACACAACAATGTCCGGGCGCGGACGCATCAGAGCATTCACGCCAAGATA > minE/1839598‑1839723 | aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:491495/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:2440906/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:263682/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:2790844/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:2815567/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:2861496/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:2898645/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:352423/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:369871/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:385486/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:1015487/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:531685/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:531719/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:667644/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:674761/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:707593/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:732422/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:793995/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:832202/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:1791579/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:1230920/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:1013374/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:1420748/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:1446500/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:1463794/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:1500446/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:1727627/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:1802972/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:1998827/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:207786/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:2139602/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:2166431/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:2319703/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGtatt > 1:2433079/1‑68 (MQ=37)aCGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGGTACGGACCCGCGta > 1:516430/1‑68 (MQ=38) ataaataaCATCATCCCCCAGATTCGGCCCGGCGTTAATTCCg > 1:209232/1‑43 (MQ=255) ataaataaCATCATCCCCCAGATTCGGCCCGGCGTTAATTCCGGACACAACAATGTCCGGGCGCGGACGCa > 1:292564/1‑71 (MQ=255) ataaataaCATCATCCCCCAGATTCGGCCCGGCGTTAATTCCGGACACAACAATGTCCGGGCGCGGACGCa > 1:519465/1‑71 (MQ=255) ataaataaCATCATCCCCCAGATTCGGCCCGGCGTTAATTCCGGACACAACAATGTCCGGGCGCGGACGCa > 1:713543/1‑71 (MQ=255) ttCGGCCCGGCGTTAATTCCGGACACAACAATGTCCGGGCGCGGACGCATCAGAGCATTCACGCCAAGATa < 1:987176/71‑1 (MQ=255) ttCGGCCCGGCGTTAATTCCGGACACAACAATGTCCGGGCGCGGACGCATCAGAGCATTCACGCCAAGATa < 1:131335/71‑1 (MQ=255) ttCGGCCCGGCGTTAATTCCGGACACAACAATGTCCGGGCGCGGACGCATCAGAGCATTCACGCCAAGATa < 1:2511470/71‑1 (MQ=255) tCGGCCCGGCGTTAATTCCGGACACAACAATGTCCGGGCGCGGACGCATCAGAGCATTCACGCCAAGATa < 1:1278385/70‑1 (MQ=255) | ACGGCCTTCCATCGCGGCGGCTACCGTACCGGAATAAATAACATCATCCCCCAGATTCGGCCCGGCGTTAATTCCGGACACAACAATGTCCGGGCGCGGACGCATCAGAGCATTCACGCCAAGATA > minE/1839598‑1839723 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |