| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 1,484,105 | 0 | A | G | 24.4% | 96.9 / 12.0 | 45 | G55G (GGA→GGG) | ydcU | predicted spermidine/putrescine transporter subunit |
| Reads supporting (aligned to +/- strand): ref base A (7/27); new base G (11/0); total (18/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.14e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAAT > REL606/1484033‑1484180 | aCTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGAttt < 2:1784367/76‑1 (MQ=255) ctgctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTtata < 1:308892/76‑1 (MQ=255) ggCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGGTTTTTATttttt > 1:60274/1‑75 (MQ=255) ggCCCACTAATGTGGTTTGGAATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTg < 1:450117/76‑1 (MQ=255) ccACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGGTTTTTACTTTTTACg > 2:262440/1‑76 (MQ=255) ccACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACg < 1:594536/76‑1 (MQ=255) aCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGAt < 1:1693098/76‑1 (MQ=255) cTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTATTTTTACGAtt > 1:616645/1‑76 (MQ=255) tAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTTTTTTTACGGttt > 2:1496146/1‑76 (MQ=255) tAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTCTTTTTACGGttt > 1:942972/1‑76 (MQ=255) aaTGTGGTTTGGCATTGTCTATTTCGGCTCACTGCTGACACTGTTATGGCAGGGGTTTTATAATTTTTATGATTTc > 1:1365505/1‑76 (MQ=255) aTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCa < 1:415990/76‑1 (MQ=255) tgtgGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCAc < 2:163138/76‑1 (MQ=255) gtgGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGGTTTTTTTTTTTGACGATTTCAcc > 2:625468/1‑76 (MQ=255) tgGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCa < 2:272716/76‑1 (MQ=255) ggTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTATTTTTTACGATTTTACCAt > 2:467770/1‑76 (MQ=255) gTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGGTTTTTTTTTTTTTCGACTTCACCAtt > 1:1316758/1‑75 (MQ=255) gCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGt < 1:1283736/76‑1 (MQ=255) aTTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTATTTTTTCGATTTTACCATGTCGGGaa > 2:845032/1‑76 (MQ=255) ttGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAAc < 1:1865365/76‑1 (MQ=255) ttGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAAc < 2:1620744/76‑1 (MQ=255) gTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTTACCCTGTCGGTAACGc > 1:1423595/1‑76 (MQ=255) gTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTTACCCTGTCGGGAACGc > 1:1296998/1‑76 (MQ=255) gTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGGTTTTTACTTTTTACGATTTTACCATTTCCGTAACGc > 1:1585227/1‑76 (MQ=255) tCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTATACTTTTGACGGTTTTACCATGTCGGTAAAGcc > 2:1417712/1‑76 (MQ=255) tCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGcc < 1:129017/76‑1 (MQ=255) aTTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGa < 1:431016/76‑1 (MQ=255) tttCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGaa < 2:442099/76‑1 (MQ=255) tttCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGaa < 1:929548/76‑1 (MQ=255) cGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTg < 2:435281/76‑1 (MQ=255) ggCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGa < 2:257571/76‑1 (MQ=255) cGCTGCTGACACTGTTATGGCAGGGGTTTTATACTTTTTACGATTTCACCATGTCGGGAAAGCCGGAACTGACGCt > 1:365946/1‑76 (MQ=255) cGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCt < 1:1062310/76‑1 (MQ=255) gctgctGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTg < 1:617901/76‑1 (MQ=255) ctgctgACACTGTTATGGCAGGGGTTTTTTACTTTTGACGATTTCACCATGTCGGGAACGCCGGGACTGACGCTgg > 2:1189572/1‑76 (MQ=255) ctgctgACACTGTTATGGCAGGGGTTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTgg > 1:502251/1‑76 (MQ=255) tgctgACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCCTGTCGGTAACGCCGGAACTGACGCCGGc > 1:255553/1‑76 (MQ=255) gctgACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCg > 2:423679/1‑76 (MQ=255) ctgACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCCTGTCGGTAACGCCGGAACTGACGCTGGCGa > 2:917691/1‑76 (MQ=255) gACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAat < 1:1682916/76‑1 (MQ=255) acacTGTTATGGCAGGGATTTTTTACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAata > 2:508390/1‑76 (MQ=255) cTGTTATGGCAGGGATTTTATACTTTTGACGATTGCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATcc < 2:146498/76‑1 (MQ=255) tGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCg < 2:1528685/76‑1 (MQ=255) gTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGt < 1:736065/76‑1 (MQ=255) tATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTgc > 2:1843737/1‑76 (MQ=255) tATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTgc < 1:665716/76‑1 (MQ=255) aTGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTgcg < 1:400489/76‑1 (MQ=255) cAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGATGCCGAATATCCGTGCGCTGt > 2:1494064/1‑76 (MQ=255) aGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGtt > 2:974480/1‑76 (MQ=255) gggATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGttt < 1:47076/76‑1 (MQ=255) gATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTaa < 2:1473974/76‑1 (MQ=255) aTTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAAt > 2:433012/1‑76 (MQ=255) | ACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAAT > REL606/1484033‑1484180 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |