| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 2,209,372 | 0 | T | G | 25.6% | 81.4 / 11.1 | 39 | V61V (GTA→GTC) | yeiN | hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base T (3/26); new base G (10/0); total (13/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.50e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GACCCAGTAATTCAATTTCTTCTTTACTTAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGTACTGCGCCCTGTTTACGAATAGTTTCTTCTACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTG > REL606/2209299‑2209446 | gACCCAGTAATTCAATTTCTTCTTTACTTAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGTAc < 2:17374/76‑1 (MQ=255)gACCCAGTAATTCAATTTCTTCTTTACTTAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGTAc < 1:466703/76‑1 (MQ=255) aCCCAGTAATTCAATTTCTTCTTTACTTAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGTACt < 2:1619150/76‑1 (MQ=255) aCCCAGTAATTCAATTTCTTCTTTACTTAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGGACt > 1:596658/1‑76 (MQ=255) aGTAATTCAATTTCTTCTTTACTTAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGGACCgcgc > 1:753429/1‑76 (MQ=255) gTAATTCAATTTCTTCTTTACTTAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGTACTGCGcc < 2:669379/76‑1 (MQ=255) aTTCAATTTCTTCTTTACTTAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGTACTGCGCCCTg < 2:1513838/76‑1 (MQ=255) aTTCAATTTCTTCTTTACTTAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGGAATGCGccccg > 2:484976/1‑74 (MQ=255) ttCAATTTCTTCTTTACTTAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGTACTGCGCCCTGt < 2:204106/76‑1 (MQ=255) ttCAATTTCTTCTTTACTTAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGGAGTGGGCCCCGt > 2:898079/1‑76 (MQ=255) cAATTTCTTCTTTACTTAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGTACTGCGCCCTGttt < 2:1254354/76‑1 (MQ=255) cAATTTCTTCTTTACTTAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGGACTGCCGCCCGttt > 2:578788/1‑76 (MQ=255) aaTTTCTTCTTTACTTAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGTACTGCGCCCTGTTTa < 1:1014134/76‑1 (MQ=255) tcttctTTACTTAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGGACTGCGCCCCGTTTTCGaa > 1:310706/1‑76 (MQ=255) ttcttTACTTAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGT‑GCCGGTACTGCGCCCCGTTTTCGCATAg > 1:1565960/1‑76 (MQ=255) tcttTACTTAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGGACTGCGCCCCGTTTTTGAATAg > 1:1191866/1‑76 (MQ=255) cttTACTTAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGGACTGCGCCCTGTTTACGAATAGt > 1:169772/1‑76 (MQ=255) ttACTTAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGTACTGCGCCCTGTTTAGGAATAGttt < 2:286597/76‑1 (MQ=255) ttAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGGACTGCGCCCCGTTTTCGAATAGGtttttt > 2:1393882/1‑75 (MQ=255) ttAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGGACTGCGCCCCGTTTACGAATAGAttcttt > 2:1006573/1‑75 (MQ=255) aaaCCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGTACTGCGCCCTGTTTACGAATAGTTTCTTCTa < 2:1255779/76‑1 (MQ=255) aaaCCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGGACTGCGCCCCGTTTTCGAATAGttttttttt > 2:1236345/1‑73 (MQ=255) aaaCCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGGACTGCGCCCCGTTTACGAATAGttttttttt > 2:1575060/1‑73 (MQ=255) aaCCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGGACTGCGCCCCGTTTACGAATAGTTtcttcttc > 2:1470516/1‑74 (MQ=255) ccACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGTACTGCGCCCTGTTTACGAATAGTTTCTTCTACTTc < 1:1343583/76‑1 (MQ=255) cACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGTACTGCGCCCTGTTTACGAATAGTTTCTTCTACTTCa < 1:1259766/76‑1 (MQ=255) aCTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGTACTGCGCCCTGTTTACGAATAGTTTCTTCTACTTCaa < 1:1558077/76‑1 (MQ=255) cTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGTACTGCGCCCTGTTTACGAATAGTTTCTTCTACTTCAAt < 2:651960/76‑1 (MQ=255) cTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGTACTGCGCCCTGTTTACGAATAGTTTCTTCTACTTCAAt < 2:1554990/76‑1 (MQ=255) cTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGTACTGCGCCCTGTTTACGAATAGTTTCTTCTACTTCAAt < 2:596658/76‑1 (MQ=255) ttCATCACGCCGCCAATAATGGCGATAGTGGCAGGGACCGCGCCCCGTTTTTTAACAATTTTTTTTTCTTCAATTg > 2:807326/1‑76 (MQ=255) catcaCGCCGCCAATAATGGCGATAGTGGCAGGGACTGCGCCCCGTTTTCGAATAGTTTTTTTTTTTTCAAttttg > 1:967214/1‑73 (MQ=255) caCGCCGCCAATAATGGCGATAGTGGCAGGGACTTCGCCCTGGTTTTGAATAATTTTTTTTTCTTCCATTGCgggc > 2:1337553/1‑74 (MQ=255) cgccAATACTCGCGATAGTGGCAGGTACTGCGCCCTGTTTACGAATAGTTTCTTCTACTTCAATTGCGGTCTGGGc < 2:588121/76‑1 (MQ=255) cgccAATAATGGCGATAGTGGCAGGTACTGCGCCCTGTTTACGAATAGTTTCTTCTACTTCAATTGCGGTCTGGGc < 1:1632639/76‑1 (MQ=255) cgccAATAATGGCGATAGTGGCAGGGACTGCGCCCTGTTTACGAATAGTTTCTTCTTCTTCAATTGCGGGCTGGGc > 1:1183222/1‑76 (MQ=255) gccAATAATGGCGATAGTGGCAGGTACTGCGCCCTGTTTACGAATAGTTTCTTCTACTTCAATTGCGGTCTGGGCa < 1:1615519/76‑1 (MQ=255) gccAATAATGGCGATAGTGGCAGGGACTGCGCCCCGTTTACGAATAGTTTCTTCTTCTTCAATTTCGGTCTGGGcg > 2:572153/1‑75 (MQ=255) cAATAATGGCGATAGTGGCAGGTACTGCGCCCTGTTTACGAATAGTTTCTTCTACTTCAATTGCGGTCTGGGCAtt < 1:323563/76‑1 (MQ=255) aataatGGCGATAGTGGCAGGGACTGCGCCCTGTTTTCGAATAATTTTTTTTTCTTCAATTGCGGTCTGGGGAttt > 1:147203/1‑76 (MQ=255) aataatGGCGATAGTGGCAGGGACTGCGCCCCGTTTTTGAATAGTTTTTTTTTCTTTTATTTCGGTCTGGGGAttt > 1:1248328/1‑76 (MQ=255) ataatGGCGATAGTGGCAGGTACTGCGCCCTATTTACGAATAGTTTCTTCTACTTCAATTGCGGTCTGGGCAtttt < 2:361842/76‑1 (MQ=255) taatGGCGATAGTGGCAGGTACTGCGCCCTGTTTACGAATAGTTTCTTCTACTTCAATTGCGGTCTGGGCATTTtg < 2:190403/76‑1 (MQ=255) gCGATAGTGGCAGGGACTGCGCCCTGTTTACGAATAGTTTTTTTTTCTTTCATTTCGGGCTGGGCCTTTTTTgggg > 2:177346/1‑75 (MQ=255) cGATAGTGGCAGGGACTGCGCCCTGTTTACGAATAGTTTTTTCTACTTCAATTGCGGGCTTGGCATTTTTTgggga > 2:1517126/1‑74 (MQ=255) gATAGTGGCAGGTACTGCGCCCTGTTTACGAATAGTTTCTTCTACTTCAATTGCGGTCTGGGCATTTTGTGGGAAc < 2:1124319/76‑1 (MQ=255) gTGGCAGGTACTGCGCCCTGTTTACGAATAGTTTCTTCTACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCa < 2:206266/76‑1 (MQ=255) cAGGTACTGCGCCCTGTTTACGAATAGTTTCTTCTACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATccc > 1:399444/1‑76 (MQ=255) aGGTACTGCGCCCTGTTTACGAATAGTTTCTTCTACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCg > 2:438237/1‑76 (MQ=255) aGGTACTGCGCCCTGTTTACGAATAGTTTCTTCTACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCg < 2:1323772/76‑1 (MQ=255) gTACTGCGCCCTGTTTACGAATAGTTTCTTCTACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTg < 1:1485298/76‑1 (MQ=255) gTACTGCGCCCTGTTTACGAATAGTTTCTTCTACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTg > 2:1623121/1‑76 (MQ=255) | GACCCAGTAATTCAATTTCTTCTTTACTTAAACCCACTTTCATCACGCCGCCAATAATGGCGATAGTGGCAGGTACTGCGCCCTGTTTACGAATAGTTTCTTCTACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTG > REL606/2209299‑2209446 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |