Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,484,105 | 0 | A | G | 32.3% | 53.0 / 13.0 | 31 | G55G (GGA→GGG) | ydcU | predicted spermidine/putrescine transporter subunit |
Reads supporting (aligned to +/- strand): ref base A (3/18); new base G (10/0); total (13/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.45e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAA > REL606/1484033‑1484179 | aCTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGAttt < 2:1054541/76‑1 (MQ=255) tgctgcTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGAttttn < 1:237954/76‑2 (MQ=255) ctgctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTtata < 1:44363/76‑1 (MQ=255) tgctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATAc < 1:1521855/76‑1 (MQ=255) gctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACt > 1:1729660/1‑76 (MQ=255) ctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTAttt > 2:1587997/1‑76 (MQ=255) ttGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACttt < 1:817898/76‑1 (MQ=255) ttGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGTTGACACTGTTATGGCAGGGATTTTATACttt < 1:781269/76‑1 (MQ=255) ttGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGGTTTTTAtttt > 1:1429086/1‑76 (MQ=255) aCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGAt < 2:1393379/76‑1 (MQ=255) cTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGAtt < 1:10876/76‑1 (MQ=255) aaTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTTTTTTTTCGAtttt > 1:689666/1‑75 (MQ=255) aaTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTc < 2:1808795/76‑1 (MQ=255) gTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATg < 2:1158530/76‑1 (MQ=255) tGGCATTGTCTATTTCGGCTCGCTGCTGACCCTGTTATGGCAGGGGCTTTTTATTTTTTTCGGCTTTCCCATGTCg > 1:1609758/1‑76 (MQ=255) ggCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTTCGATTTTACCCTGTCgg > 1:1539606/1‑76 (MQ=255) gCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGt < 1:1330650/76‑1 (MQ=255) aTTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTaa < 1:319605/76‑1 (MQ=255) tCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTATTATTTTTTCGATTTTACCATGTCGGGAACGcc > 1:253072/1‑76 (MQ=255) cTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCg < 1:805669/76‑1 (MQ=255) tttCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTCACCATGTCGGTAACGCCGGaa > 2:1574678/1‑76 (MQ=255) tCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTATTTTTTACGATTTTACCCTGTCGGGAACGCCGGaaac > 1:72985/1‑74 (MQ=255) tCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACt < 1:1504901/76‑1 (MQ=255) cGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTCACCCTGTCGGGAACGCCGGaaacg > 1:698951/1‑73 (MQ=255) cGGCTCGCTGCTGACACTGTTATGGCAGGGGGTTTTTTCTTTTTATGATTTCTCCCTGTCGGGAACGCCCGAAATg > 1:1692748/1‑76 (MQ=255) cTCGCTGCTGACACTGTTATGGCAGGGGTTTTATACTTTTGACGATTTCACCATGTCGGGAACCCCGGAACTGACg > 2:415274/1‑76 (MQ=255) cGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTCACCATGTCGGTAACGCCGGAACGGACGcc > 1:361121/1‑75 (MQ=255) cGCTGCTGACACTGTTATGGCAGGGGTTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAAATTACGCt > 2:1719705/1‑76 (MQ=255) cGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCt < 2:245619/76‑1 (MQ=255) acTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATc < 2:368504/76‑1 (MQ=255) gTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGt < 1:74556/76‑1 (MQ=255) gCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTg > 2:358444/1‑76 (MQ=255) aGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGtt > 1:1315419/1‑76 (MQ=255) gggATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGttt < 1:1865248/76‑1 (MQ=255) gATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTaa > 1:1850054/1‑76 (MQ=255) | ACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAA > REL606/1484033‑1484179 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |