| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 1,484,105 | 0 | A | G | 34.1% | 56.0 / 25.6 | 44 | G55G (GGA→GGG) | ydcU | predicted spermidine/putrescine transporter subunit |
| Reads supporting (aligned to +/- strand): ref base A (6/23); new base G (15/0); total (21/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.36e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAAT > REL606/1484032‑1484180 | tACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGTTGCTGCCCCTTTTCTGGCGGGGttt > 1:1045681/1‑76 (MQ=255)tACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGtt > 1:1350118/1‑76 (MQ=255) aCTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGAttt < 2:599873/76‑1 (MQ=255) ctgctgCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGtttt > 1:214131/1‑76 (MQ=255) gctgctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTtat < 1:614630/76‑1 (MQ=255) gctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACt < 1:150693/76‑1 (MQ=255) gctTGGCCCACTAATGTGGTTTGGCATTGCCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACt > 2:1381441/1‑76 (MQ=255) ctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACtt < 1:557736/76‑1 (MQ=255) tGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACtttt < 2:1241358/76‑1 (MQ=255) ggCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTg < 2:489077/76‑1 (MQ=255) gCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACtttttt > 1:1589547/1‑74 (MQ=255) gCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGa < 1:134189/76‑1 (MQ=255) ccACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACg < 1:661120/76‑1 (MQ=255) aaTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTc < 2:1043885/76‑1 (MQ=255) gggTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCa < 2:618728/75‑1 (MQ=255) ttGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTc < 2:1443592/76‑1 (MQ=255) ttGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTc < 2:375855/76‑1 (MQ=255) tGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTCCCCCTCTCg > 1:1459399/1‑76 (MQ=255) tGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTCACCATGTCg > 1:574871/1‑76 (MQ=255) tGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCg < 2:1771708/76‑1 (MQ=255) ggCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTTCGATTTTACCATGTCgg > 1:1678257/1‑76 (MQ=255) ggCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTTACCATTTCgg > 1:236867/1‑76 (MQ=255) ggCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTATACTTTTGACGATTTTACCCTGTCgg > 1:1462161/1‑76 (MQ=255) gCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTATTTTTTTCGGTTTCACCCTGTCggg > 1:860130/1‑75 (MQ=255) cATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTATTTTTTTCGATTTTACCATGTCggga > 1:1290807/1‑74 (MQ=255) cATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTa < 2:901967/76‑1 (MQ=255) aTTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTaa < 1:1031238/76‑1 (MQ=255) ttGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAAc < 1:1218482/76‑1 (MQ=255) ttGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAAc < 1:1182841/76‑1 (MQ=255) tGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTTACCATGTCGGGAACg > 1:110183/1‑76 (MQ=255) gTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTATTTTTACGATTTTACCATGTCGGTAAAGc > 1:1621446/1‑76 (MQ=255) gTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTTACCATGTCGGTAACGc > 1:416494/1‑76 (MQ=255) tCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGcc < 1:1687681/76‑1 (MQ=255) tttCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTAAGATTTTACCATGTCGGGAAAGCCggga > 2:495192/1‑74 (MQ=255) tttCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTGACGATTTTACCATGTCGGTAAAGCCggga > 2:1040443/1‑74 (MQ=255) ttCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTCTTTTTATGGTTTCTCCATCTCGGTAAACCCgggaa > 2:354466/1‑73 (MQ=255) gCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGAc > 1:1862948/1‑76 (MQ=255) tCGCTGCTGACACTGTTATGGCAGGGGTTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGGACTGACGc > 1:932563/1‑76 (MQ=255) cGCTGCTGACACTGTTATGGCAGGGGTTTTATACTTTTGACGATTTTACCCTGTCGGTAACGCCCGGACTGCCGcg > 2:628594/1‑75 (MQ=255) ctgctgACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTgg < 2:917745/76‑1 (MQ=255) tgctgACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGc < 2:1118057/76‑1 (MQ=255) gctgACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCg < 1:556007/76‑1 (MQ=255) ctgACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTTGGTANCGCCGGAACTGACGCTGGCGa > 2:296250/1‑76 (MQ=255) acTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATc > 2:159792/1‑76 (MQ=255) ttATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTg < 2:1327487/76‑1 (MQ=255) cAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGt > 2:642731/1‑76 (MQ=255) aGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGtt > 1:1686408/1‑76 (MQ=255) gATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTaa < 2:658994/76‑1 (MQ=255) gATTTTATACTTTTGACCATTTCACCACGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTaa > 2:148537/1‑76 (MQ=255) aTTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAAt < 1:696460/76‑1 (MQ=255) | TACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAAT > REL606/1484032‑1484180 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |