| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,959,174 | G→A | 100% | G299G (GGC→GGT) | fbaA ← | fructose‑bisphosphate aldolase, class II |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,959,174 | 0 | G | A | 90.4% | 97.7 / 0.6 | 51 | G299G (GGC→GGT) | fbaA | fructose‑bisphosphate aldolase, class II |
| Reads supporting (aligned to +/- strand): ref base G (4/1); new base A (24/22); total (28/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.62e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
CTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCTTTGTAGTAGTTCAGAACGCCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTC > minE/1959107‑1959213 | cTTCGCCTTTCGGGTTTCCCAGCTGACCCTGCAGATAAGCTTCGTTCGCTTTGTAGTAGTTCAGAACg < 1:276956/68‑1 (MQ=255)cTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCTTTGTAGTAGTTCAGAACg < 1:726915/68‑1 (MQ=255)cTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCTTTGTAGTAGTTCAGAACg < 1:1626444/68‑1 (MQ=255)cTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCTTTGTAGTAGTTCAGAACg < 1:380759/68‑1 (MQ=255) tGACCCTGCAGATAAGCTTCGTTCGCTTTGTAGTAGTTCAGAACg < 1:210232/45‑1 (MQ=255) gCAGGTGAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:7250/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTTTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:2980771/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:200322/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:2233008/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:2029302/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:1042252/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:2261950/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:2278199/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:2529852/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:2787165/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:2929410/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:3434956/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:733216/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:89837/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:973566/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:982606/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:999224/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:1313681/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:1855275/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:1735115/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:1838481/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:1771902/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:1905179/1‑55 (MQ=255) gCAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCa > 1:1818750/1‑55 (MQ=255) cAGGTAAGCTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCAtt < 1:218172/56‑1 (MQ=255) gATAAGCTTCGTTCGCTTTGTAGTAGTTCAGAACGCCTTCCCAGGTTGCCCATTGggtatcggtatc < 1:293455/67‑1 (MQ=255) aGCTTCGTTCGCTTTGTAGTAGTTCAGAACGCCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGtt > 1:1259209/1‑69 (MQ=255) aGCTTCGTTCGCTTTGTAGTAGTTCAGAACGCCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGtt > 1:1822309/1‑69 (MQ=255) aGCTTCGTTCGCTTTGTAGTAGTTCAGAACGCCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGtt > 1:1791243/1‑69 (MQ=255) aGCTTCGTTCGCTTTGTAGTAGTTCAGAACGCCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGtt > 1:289341/1‑69 (MQ=255) cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc < 1:1526647/68‑1 (MQ=255) cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc < 1:1717966/68‑1 (MQ=255) cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc < 1:1714090/68‑1 (MQ=255) cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc < 1:618753/68‑1 (MQ=255) cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc < 1:1594137/68‑1 (MQ=255) cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc < 1:29207/68‑1 (MQ=255) cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc < 1:1397035/68‑1 (MQ=255) cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc < 1:836738/68‑1 (MQ=255) cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc < 1:1230458/68‑1 (MQ=255) cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc < 1:1179747/68‑1 (MQ=255) cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc < 1:3422099/68‑1 (MQ=255) cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc < 1:2792009/68‑1 (MQ=255) cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc < 1:2600841/68‑1 (MQ=255) cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc < 1:2584257/68‑1 (MQ=255) cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc < 1:2445671/68‑1 (MQ=255) cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc < 1:1845109/68‑1 (MQ=255) cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc < 1:18601/68‑1 (MQ=255) cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc < 1:2090839/68‑1 (MQ=255) cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc < 1:2083953/68‑1 (MQ=255) cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTc < 1:1893435/68‑1 (MQ=255) cTTCGTTTTCTTTGTAGTAGTTCAGAACACCTTACCAGGTTGCCCATTGGGTATCGGTATCGATGTTc < 1:3290422/68‑1 (MQ=255) | CTTCGCCTTTCGGGTTACCCAGCTGACCCTGCAGATAAGCTTCGTTCGCTTTGTAGTAGTTCAGAACGCCTTCCCAGGTTGCCCATTGGGTATCGGTATCGATGTTC > minE/1959107‑1959213 |
| 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: ‑ |