| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 1,484,105 | 0 | A | G | 27.3% | 65.8 / 11.9 | 33 | G55G (GGA→GGG) | ydcU | predicted spermidine/putrescine transporter subunit |
| Reads supporting (aligned to +/- strand): ref base A (3/21); new base G (9/0); total (12/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.70e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCG > REL606/1484031‑1484171 | ttACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGAt < 1:355530/76‑1 (MQ=255) tACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGAtt < 1:422609/76‑1 (MQ=255) gctgctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGATTTtat < 2:723153/76‑1 (MQ=255) nttGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGATTTtttatt > 1:1568666/2‑74 (MQ=255) ggCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGATTTTATACTTTTg < 1:329947/76‑1 (MQ=255) gCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGGGGTTTTACtttttt > 2:174465/1‑74 (MQ=255) cccACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGATTTTATACTTTTGAc < 2:220971/76‑1 (MQ=255) ccACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGATTTTATACTTTTGACg < 1:1454465/76‑1 (MQ=255) aCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGGTTTTTTACTTTTTACGAt > 2:1594007/1‑76 (MQ=255) aCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGATTTTATACTTTTGACGAt < 1:939834/76‑1 (MQ=255) tAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGATTTTATACTTTTGACGAttt < 1:1252650/76‑1 (MQ=255) aaTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGGTTTTTTTATTTTTACGAtttt > 1:649427/1‑75 (MQ=255) aaTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGATTTTATACTTTTGACGATTTc < 2:1540819/76‑1 (MQ=255) aTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGGGTTTTTACTTTTGATGATTTct > 2:1107137/1‑75 (MQ=255) aTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGATTTTATTCTTTTGAAGAttttt > 2:1107144/1‑74 (MQ=255) gtggtgTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGATTTTATACTTTTGACGATTTCAcc < 2:469511/76‑1 (MQ=255) tgGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGGTTTTTTATTTTTTACGATTTTACCa > 1:1090739/1‑76 (MQ=255) tgGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGATTTTATACTTTNGACGATTTCACCa < 1:1232485/76‑1 (MQ=255) ggTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGATTTTTTTCTTTTTGCGATTTTAccct > 1:1035269/1‑74 (MQ=255) ttGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGGTTTTTTACTTTTTACGATTTCACCATGTc > 1:5412/1‑76 (MQ=255) gCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGATTTTATACTTTTGACGATTTCACCATGTCGGt < 1:702839/76‑1 (MQ=255) cATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGATTTTATACTTTTGACGATTTCACCATGTCGGTa < 1:1367115/76‑1 (MQ=255) ttGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGGtttttt > 2:1538712/1‑45 (MQ=255) ttGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGGTTTTTTACTTTTTACGATTTCACCATGTCGGTAAc > 2:1120624/1‑76 (MQ=255) ttGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGATTTtat < 1:1538712/47‑1 (MQ=255) ttGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGATTTTATACTTTTGACGATTTCACCATGTCGGTAAc < 2:1098290/76‑1 (MQ=255) tCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGGTTTTTTAATTTTGAAGATTTTAACATTTCCGTAAcccc > 1:1009228/1‑76 (MQ=255) tATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGGTTTTTTTCTTTTTTCGATTTTCCCCTTTCGGGAAcccccg > 2:764718/1‑74 (MQ=255) tATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGGGTTTTTTCTTTTTACGAATTCACCATGTCGGGAAcccccg > 2:1556137/1‑74 (MQ=255) aTTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGGTTTTTTACTTTTTACGATTTTACCATTGCCGGAACGCCggg > 2:884468/1‑75 (MQ=255) ttCGGCTCGCTGCTGACACTGTTATGGCA‑GGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAAc < 2:880040/76‑1 (MQ=255) cGGCTCGCTGCTGACACTGTTATGGCA‑GGGGTTTTTTTCTTTTTACGATTTTACCATGTCGGTAACGCCGGAACTg > 1:1461224/1‑76 (MQ=255) cGGCTCGCTGCTGACACTGTTATGGCA‑GGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTg < 1:197773/76‑1 (MQ=255) ggCTCGCTGCTGACACTGTTATGGCA‑GGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGa < 2:21125/76‑1 (MQ=255) gCTCGCTGCTGACACTGTTATGGCA‑GGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGAc < 2:924403/76‑1 (MQ=255) cTCGCTGCTGACACTGTTATGGCA‑GGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACg < 2:540678/76‑1 (MQ=255) tCGCTGCTGACACTGTTATGGCA‑GGGATTTTATACTTTTGACGATTTCACCATGTCGGGAAACCCGGAAATGAcag > 2:427904/1‑74 (MQ=255) gTTATGGCA‑GGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGt < 1:1722661/76‑1 (MQ=255) aTGGCA‑GGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTgcg > 2:1412762/1‑76 (MQ=255) | TTACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCA‑GGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCG > REL606/1484031‑1484171 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |